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ABSTRACT 


In this review/tutorial we explore planetary nebulae as a stage in the evolution of 
low-to-intermediate-mass stars, as major contributors to the mass and chemical enrich- 
ment of the interstellar medium, and as astrophysical laboratories. We discuss many 
observed properties of planetary nebulae, placing particular emphasis on element abun- 
dance determinations and comparisons with theoretical predictions. Dust and molecules 
associated with planetary nebulae are considered as well. We then examine distances, 
binarity, and planetary nebula morphology and evolution. We end with mention of 
some of the advances that will be enabled by future observing capabilities. 


Keywords: ISM: abundances, planetary nebulae: general, stars: evolution, galaxies: 
evolution 


1. INTRODUCTION 


In 1764, while observing the sky toward the constellation Vulpecula, the French astronomer Charles 
Messier noted a “nebula without star" that “appears of oval shape” (Messier 1771), adding it as 
number 27 in his catalog of non-stellar objects not to be confused with his main interest, comets; 
thus was discovered the first planetary nebula, known today as the Dumbbell Nebula (or M27 or 
NGC 6853). Eighteen years later, shortly after beginning his tenure as Court Astronomer to King 
George III of England, William Herschel was scanning the heavens for double stars when he, too, 
found a curious non-stellar object, one we now call the Saturn Nebula (or NGC 7009). Herschel 
dubbed objects like these “planetary nebulae,” and he puzzled over their nature for the rest of 
his life (Hoskin 2014). We are fortunate nowadays to know much more, though by no means, all, 
about them: planetary nebulae (PNe) are ephemeral manifestations of the dynamic nature of stellar 
evolution and the ongoing enrichment of the Galaxy's reservoir of star-forming matter. They are the 
penultimate stage in the lives of multitudes of stars, including, possibly, the Sun (Boffin & Jones 2019, 
p. 94). Fig. 1 is an artistic montage of 22 PNe displaying some of the variety of their shapes and sizes. 


Since the time of Messier and Herschel, PNe have become indispensable tools for studying stellar 
evolution, galactic chemical evolution, and the interstellar medium. Studies of PNe also contribute 
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to a variety of other astrophysical topics. PNe progenitor stars contribute the majority of the 
matter that forms the ISM (Dorschner & Henning 1995; Edwards, Cox, & Ziurys 2014). They are 
also excellent tracers of stellar populations in galaxies (e.g., Buzzoni & Arnaboldi 2006; Hartke 
et al. 2020) and serve as valuable test particles in studies of galaxy kinematics (Aniyan, Freeman, 
& Arnaboldi 2018; Aniyan et al. 2020) and galaxy cluster and merger dynamics (Gerhard et al. 
2007). Furthermore, PNe are fertile sites for detecting important molecular species: recently the first 
detection of the cosmologically significant molecule, HeH* via its fundamental rotational transition, 
was reported in NGC 7027 (Güsten, Wiesemeyer, & Neufeld 2019), and soon thereafter, Neufeld et 
al. (2020, 2021) detected additional rovibrational HeH* lines, along with CH* emission seen for the 
first time in an astronomical source. 


A significant fraction of the high luminosity of PNe is concentrated in the Ott emission line at 
5007À (see 85.2.1), allowing them to be observed at extragalactic distances. The first extragalactic 
PNe were detected by Baade (1955), who discovered five in M31. Subsequent searches in M31 using 
increasingly sensitive detectors revealed 311 PNe (Ford & Jacoby 1978). Merrett et al. (2006) found 
more than 2600 PNe; the latest survey in M31 identifies over 4000 (Bhattacharya, Arnaboldi, & 
Hartke 2019). PNe have also been identified in dozens of external galaxies, including members of 
the Virgo and Fornax clusters (Ford, Peng, & Freeman 2002; Spriggs et al. 2021) at d~20 Mpc, and 
even out to the Coma cluster at d~100 Mpc (Gerhard et al. 2005). 


How many PNe are there in the Milky Way Galaxy? We don't know! As of this writing, the 
HASH database (Parker et al. 2016) contains 2667 “true,” 447 “likely,” and 681 “possible” PNe; 
Parker (2020) reviews the status of modern searches. According to Awang Iskandar et al. (2020), 
64% of “likely” and 41% of “possible” PNe are “true” PNe, predicting a total of 3232. The current 
census likely represents only 15 — 30% of the full count (Frew 2016), implying a population between 
roughly 11,000 and 22,000. Jacoby et al. (2010) discuss various estimates of the total number of PNe 
(all much larger than the number detected), and enumerate possible reasons for the shortfall (e.g., 
interstellar dust obscuration, survey biases, poorly searched regions of the sky). As one example, 
Hong et al. (2021) report near-infrared observations of two probable PNe at projected distances 
<20 pc from the Galactic center and behind more than 20 magnitudes of visual extinction; these 
would be the first detected in the nuclear stellar disk, a disk of stars in rotation around Sgr A*. The 
Galactic distribution of the known PNe is shown in Fig. 2. 


There have been several valuable reviews of PNe over the years, such as those by Kaler (1985), 
Kwok (1994), Kwitter et al. (2014), and Zijlstra (2015). Recent IAU Symposia on PNe are #234 
(Barlow & Méndez 2006), #283 (Manchado, Stanghellini, & Schónberner 2012), and #323 (Liu, 
Stanghellini, & Karakas 2016). Books and monographs concentrating on PNe include Aller & Liller 
(1968), Aller (1971), Pottasch (1983), Aller (1987), Gurzadyan (1997), Kwok (2000), and Osterbrock 
& Ferland (2006). Frew (2008) provides an excellent, detailed study of PNe in the solar neighborhood. 


In this review /tutorial we will show how PNe reveal the nucleosynthesis histories of their stars, and 
examine their important impact on the chemical evolution of the Milky Way and other galaxies. In 
§2 we briefly review the evolution of PN progenitor stars (low-to-intermediate-mass stars: LIMS: 0.8 
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— 8 Mo) from the main sequence to the PN stage. Detailed discussion of the evolution and properties 
of PN central stars is not included here; we recommend Weidmann et al. (2015, 2020) for a current 
overview, and Iben (1995) for a historical, physics-heavy, deep dive. Werner (2012) and Ziegler et al. 
(2012) discuss the evolution of H-poor and H-rich central stars! , respectively, and Hajduk et al. (2015) 
explore the temporal variability of nebular fluxes as a way to study young central star evolution. In 83 
we examine in detail the chemical abundances measured in PNe in the Milky Way and in a selection of 
Local Group and more distant galaxies. We also discuss the molecular and dust components of PNe. 
The extensive nebular astrophysics involved in calculating these abundances is not included here, but 
is comprehensively discussed in excellent recent reviews by Peimbert, Peimbert, & Delgado-Inglada 
(2017) and García-Rojas (2020). Then in 84 we discuss radial abundance gradients in galaxy disks, 
and go on to compare observed abundances with theoretical predictions. Distance determinations for 
PNe are explored in 85, and in 86 we consider the important question of binarity in PN central stars. 
Nebular morphology and evolution are discussed in 87. In 88 we summarize and touch on anticipated 
future developments in PN research. Appendix A contains a list of databases and catalogs of PNe 
and central stars, and Appendix B lists references to papers containing model-predicted yields and 
surface abundances of LIMS. 


2. STELLAR EVOLUTION AND PLANETARY NEBULA FORMATION 
2.1. From the Main Sequence to the Asymptotic Giant Branch 


The brief summary that follows is designed to get us to the point where we can discuss how nuclear 
products get from the stellar interior to the envelope and then to the interstellar medium. It is drawn 
from excellent reviews of the evolution of asymptotic giant branch (AGB) stars by Herwig (2005) 
and Karakas & Lattanzio (2014), to which the reader should turn for greater detail. To gain an 
historical perspective on the development of this scenario, the reader should begin with the papers 
by Iben & Truran (1978) and Renzini & Voli (1981). The evolutionary stages described below are 
shown in the context of an HR diagram in Fig. 3 [adapted from Herwig (2005)], where the evolution 
of a 2 Mo star is illustrated. 


When a LIMS exhausts its core H, a H-burning shell surrounding the now-He core ignites and slowly 
moves outward, leaving behind its He product. The star's outer envelope then expands and cools as 
the star moves through the Hertzsprung gap and begins to ascend the red giant branch (RGB) in the 
HR diagram. While the energy needed for this expansion may come from the H-burning shell, there 
is no firm consensus among specialists regarding the actual mechanism by which the star expands 
and becomes a red giant (Miller Bertolami 2021). Once it reaches the upper end of the RGB, He 
ignition occurs in the core, either in a runaway flash under degenerate conditions (if M<2 Mo) 
or more gently (if M22 M). The star now settles into a relatively stable giant-like configuration 
(expanded envelope, dense core), as He is slowly burned to C and O in the core. Following exhaustion 
of He in the core, a He-burning shell surrounding the CO core and interior to the H shell ignites 
and begins burning He into C and O, as the star moves up along the AGB track. The two shells 
are separated by an intershell region, now rich in He, and it is at this point that thermal pulses, 


! This distinction results from the relative timing of the progenitor star's departure from the AGB and its last helium- 
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triggered by degenerate He-burning in the He shell, commence. 


During the evolutionary process just summarized, the He and heavy element products of nuclear 
burning are transported to the outer envelope and from there deposited through mass loss into the 
interstellar medium, thereby seeding the next stellar generation with the enriched materials. Trans- 
porting the nuclear products from the stellar interior to the star's outer envelope is accomplished 
by convective processes involving two or three dredge-up events depending upon the stellar mass. 
Because opacity is high in the outer envelope of a giant star, efficient heat transfer to the surface 
requires convection. Thus, first dredge-up (FDU) occurs as the evolved star ascends the RGB and 
the inner surface of the convective envelope extends downward into the He-rich intershell region. At 
that time hydrogen-burning products such as *He, 13C, 14N are mixed up into the heretofore pristine 
envelope, thereby enriching it with these products of H fusion. Interestingly, model predictions of 
envelope abundances following FDU in the higher luminosity RGB stars markedly disagree with 
observations. For example, the '*C/'’C ratio as well as the abundance of Li are both observed to be 
significantly lower than the models predict. This result suggests the occurrence of a process called 
extra mixing which transfers additional products from the H-burning shell up into the convective 
envelope (see Karakas & Lattanzio 2014, and references therein). 


Following He exhaustion and the subsequent onset of He-shell burning, the star begins to ascend 
the AGB. Stars with masses exceeding 3-4 Mo undergo a second dredge-up (SDU) as the convective 
envelope reaches deep inside the intershell region to move more of the same H-burning products into 
the outer atmosphere. Less massive stars skip the SDU stage. 


At this point the stellar structure, moving outward from the center, includes a carbon-oxygen core, 
a He-burning shell, a He-rich intershell region, a H-burning shell, and a convective envelope now 
polluted by H-burning products (see Fig. 14 in Karakas & Lattanzio 2014). As the star proceeds up 
the AGB, the gravitational effects associated with the increasing mass of the contracting core cause 
the He-shell to become compressed and thermally unstable. A He flash (thermal pulse) ensues, result- 
ing in envelope expansion as well as initiation of convection and homogenization within the intershell. 


Following the He-shell flash, energy produced by the flash causes the intershell region to expand 
and cool, while at the same time, burning in the H-shell ceases. The increased opacity of the 
intershell region resulting from the temperature reduction allows the base of the convective envelope 
to penetrate deeply into the region, causing He-burning products, especially !?C, to be transported 
to the outer envelope during a series of pulsations of the AGB star, a process called third dredge-up 
(TDU). The amount of dredged-up !*C can be significant, and enough to result in C/O>1 in the 
convective envelope, i.e., a carbon star. 


Such thermal pulses continue through dozens of cycles, with the precise number dependent upon 
core mass and metallicity, and each pulse results in further dredge-up. During the interpulse phase, 
the base of the convective envelope lies near or partially within the upper portion of the H-shell 
where temperatures can reach 50 x 10°K in stars of M > 4M. With the envelope now containing an 
increasing amount of ^C from TDU, CNO processing may ensue, converting C to 14N and forcing 
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C/O to values less than unity, i.e., an O-rich envelope. This process is referred to as *hot bottom 
burning? (HBB) and is discussed in more detail in 883 and 4. 


During the post main sequence evolution just described, the star has gone through phases of both 
gentle mass loss as well as the punctuated losses that occur in connection with the thermal pulses. 
'The moment when a star leaves the AGB to begin the journey to the hot, blue side of the HR diagram 
(see Fig. 3) is not obvious. The essence of this milestone is that the thermal pulses described above 
have ceased (as have TDU and HBB), profuse mass loss [up to 107° Mọ/yr; Decin et al. (2019)] 
has reduced the envelope mass to a small fraction of the core mass [~1%; Miller Bertolami (2016, 
hereafter M3B16)], and the envelope has become detached (Lagadec 2016). 


2.2. From AGB to Planetary Nebula 


Observationally, the early stages of the AGB-PN transition are typically marked by stellar spectral 
types of supergiant F or G, if the star can be seen; more often the star is surrounded by its lost 
envelope, obscured by dust of its own making, and optically visible only via scattered light. The 
dust composition will be dominated by either carbon or oxygen (silicate) chemistry, depending on 
the C/O ratio. As just described, this will be mass dependent: below ~ 1.5 Mo” too few thermal 
pulses occur to modify the surface composition toward being C-rich. Up to ~5 Mọ, C/O will be 
enhanced. But above —4.5 Mo HBB sends the surface C/O ratio back below unity and enhances 
surface N (see 84.2.2) (Karakas 2014; Karakas & Lugaro 2016). Both O-rich and C-rich stars will 
have CO in their atmospheres; the less abundant element will be locked into CO, leaving the more 
abundant element able to form dust grains. Dust in PNe is discussed in detail in 83.5. 


The time it takes a star leaving the AGB to move across the HR diagram to the blue depends 
steeply on the core mass. M3B106 defines the “crossing time” as the time it takes for the star to 
evolve from an effective temperature ~ 104 K to its eventual maximum effective temperature ~10° K 
(see his Fig. 8). Crossing times range from ~10* years for initially solar-mass stars to mere decades 
for stars with initial masses ~4 Mọ. Fig. 4 shows the crossing time for solar-metallicity stars of 
varying main sequence masses based on data from his Table 3. 


Depending on the phase of the helium shell flash cycle during which the star departs the AGB, 
a final late or very late thermal pulse can occur (Blócker 2001; Werner & Herwig 2006), looping 
the star back across the HR diagram to cooler temperatures (see Fig. 3), leading to H-deficiency 
and born-again status; see Iben, Kaler, & Truran (1983). Such PNe are rare; in their paper on 
the probable binary central star of born-again PN A 30, Jacoby et al. (2020) list only seven others 
(A 58, A 78, GJJC-1, Sakurai’s Object, WR72, IRAS 15154-5258, and HuBi 1). Weidmann et al. 
(2020) find that 2/3 of PN central stars in their catalog are H-rich and 1/3 are H-poor. Of binary 
central stars, almost 8096 are H-rich, which may indicate that binary evolution can somehow dis- 
rupt the sequence of events that lead to an H-poor central star. We discuss central star binarity in 86. 


? Mass limits for nucleosynthesis processes are metallicity dependent (see $4); values given here apply at solar metallicity 


(Z=0.014). 
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Nomenclature during the blueward HR. diagram crossing of a PN-to-be is murky: post-AGB 
(pAGB) and pre-PN (PPN) are both used, though the former also applies to stars that will not 
produce PNe due to a lack of either sufficiently dense nebular material or ionizing photons — essen- 
tially a mismatch between the expansion timescale and the evolution timescale (i.e., crossing time; 
Renzini 1981). The beginning of the PN stage is fairly easy to define: as the star heats to 220,000 K, 
widespread hydrogen ionization (and thus, emission) begins in the nebula, and at ~30,000 K, oxy- 
gen can be doubly ionized, producing the signature [O III] 5007À line in the spectrum of the new PN. 


3. CHEMICAL ABUNDANCES IN PLANETARY NEBULAE 


The gas forming a PN is the debris ejected by the progenitor star as it nears the end of its lifetime. 
Thus, the abundances of chemical elements in that gas which are not synthesized by the star provide 
information about the composition of the interstellar medium at the time of star formation. At 
the same time, abundances of those elements in the ejected gas which are synthesized within the 
star supply crucial information regarding the nuclear processes that have occurred during the star's 
lifetime. 


Early examples of some of the first work on PN abundances can be found in Aller & Menzel (1945), 
Kaler (1970), Peimbert & Torres-Peimbert (1971), and Peimbert (1978). These important studies 
illustrated what could be done in a more challenging period when dependable values for atomic 
constants as well as estimations of corrections for unseen ionization stages had only just begun to 
become available. While abundance studies of today can take advantage of the significant evolution 
that has occurred in these two areas, each of the earlier studies helped, by their example, to set in 
motion what has become a vibrant industry. 


We begin the discussion by taking an extensive look at important large surveys of PNe in the 
Milky Way Galaxy (MW), initially focusing on elements which are products of hydrogen and helium 
burning, i.e., He, C, N, O, Ne, S, CP, and Ar. Next, we list numerous smaller but important studies 
of MW PNe where more extensive work on individual PNe may be found. Following that, we discuss 
abundance studies of elements on rows 3 and 4 of the periodic table, including Fe-peak and s-process 
elements. Finally, we explore abundances in extragalactic PNe as well as molecules and dust in 
post-AGB objects. 


Several good review papers have been published recently addressing the subjects of nebular physics 
and abundance determinations from recombination and collisionally excited lines*. Readers who 
are unfamiliar with these topics or who desire a refresher course are referred to articles by Kwitter 
& Henry (2011), Pérez-Montero (2017), Peimbert, Peimbert, & Delgado-Inglada (2017); Peimbert 
(2019), and García-Rojas (2020). 


3 Unlike He-burning products such as C, O, Ne, Mg, Si, S, and Ar, Cl is not synthesized directly from alpha reactions. 
?5CL the most abundant Cl isotope, forms when ??S undergoes two neutron captures to produce 34S which in turn 
captures a proton. ?'Cl, a less abundant species, forms by a Bt decay of ?" Ar, itself a product of neutron capture 
onto ?9 Ar (Clayton 2003). 

^A labelled sample spectrum from the near-UV to the near-IR can be found in Dufour et al. (2015); 
see their Fig. 3. See also the templates provided on the Gallery of Planetary Nebula Spectra website: 
https:/ /web.williams.edu/ Astronomy /research/PN /nebulae/legend.php. 
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The reader should keep in mind that persistent problems exist regarding the computation of both 
ionic and elemental abundances, such as the abundance discrepancy factor (ADF) and ionization 
correction factor (ICF). The ADF is addressed most recently by Peimbert (2019) and García-Rojas 
(2020) as well as in an earlier study by Kholtygin (1998) (see also §6.2 here), while the ICF is most 
recently treated by Delgado-Inglada, Morriset, & Stasiriska (2014) and Amayo, Delgado-Inglada & 
García-Rojas (2020, Na, K, Ca). 


Finally, a word of caution. The common use of oxygen as the standard metallicity indicator is 
increasingly controversial these day, because some observations and theoretical models imply that 
this element can be destroyed or synthesized during AGB evolution. If you desire to explore that 
problem before proceeding here, you can find a detailed discussion within the context of galaxy 
abundance gradients in 84.1. Otherwise, if you prefer to wait until you get to that section in the 
natural order, then just keep reading on here. 


3.1. Large Elemental Abundance Surveys of MW Planetary Nebulae 


The goal of this section is to present an extensive listing of individual PN abundance surveys in the 
MW, where each study comprises a relatively large number of objects (usually > 30 for disk objects 
and a bit less than that for bulge PNe). In addition to sample size, two additional selection criteria 
for choosing a survey for analysis here were: 1) the abundances presented in the paper must have been 
computed by the paper’s authors; and 2) line strengths used in the abundance determinations should 
be those observed and measured by the authors themselves, or have been adopted from papers in the 
literature after the authors have critically evaluated them. These two criteria ensure a reasonable 
level of homogeneity within each sample. In addition to the disk and bulge surveys, we also present 
a listing of published abundance studies for each of the 13 documented Galactic halo PNe. Generally 
speaking, the He abundances were derived from recombination lines while the abundances of all other 
elements were determined from collisionally excited lines. A discussion regarding the relevance of 
these results to stellar evolution and nucleosynthesis is postponed until section 4.2. 


3.1.1. MW Disk PN Surveys 


The surveys chosen for abundance analysis and comparison here, as listed in Table 1, are those 
of Barker (1978b), Aller & Czyzak (1983); Aller & Keyes (1987), Kingsburgh & Barlow (1994), 
Perinotto, Morbidelli, & Scatarzi (2004), Stanghellini et al. (2006), Girard, Kóppen, & Acker (2007), 
Pottasch & Bernard-Salas (2010), Maciel et al. (2017), and Kwitter & Henry (2020, hereafter KH20). 
We removed any PNe listed as bulge objects by Chiappini et al. (2009)?, as well as known halo PNe, 
from several of the samples in order to ensure that each comprised only disk PNe. As stipulated 
above, the abundances of all PNe within a sample were determined homogeneously by the authors, 
ie. the same computational method and software were employed consistently across that sample. 
Most authors determined elemental abundances by first computing one or more ion abundances of an 
element directly from spectral data and then applying an ionization correction factor that accounts 
for the unobserved ions of an element. 


5 [n selecting bulge PNe, Chiappini et al. (2009) used the set of criteria proposed by Stasiriska & Tylenda (1994), i.e., 
objects must be located within 10 degrees of the Galactic center, possess an angular diameter smaller than 20 arcsec, 
and have radio fluxes at 5 GHz of less than 100 mJy. 
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The upper portion of Table 1 lists the basic characteristics of each disk survey that we evaluated, 
arranged according to publication date from earliest to latest. The survey name is listed in column 1 
and linked to a specific reference in the footnote. The number of PNe considered in each survey 
is given in column 2, while the relevant wavelength range or ranges of the spectroscopic data upon 
which each survey was based are listed in column 3. The medians of elemental number abundances 
relative to hydrogen are given in columns 4-11, while solar abundances from Asplund et al. (2009) 
are provided in the final row of the table. Survey sizes range from 32 to 230 disk objects. Six studies 
derived their abundances based entirely upon spectral coverage within the optical range, while three 
extended their coverage into the UV and, in one case, the IR. Note that with the exception of the 
surveys by Perinotto, Morbidelli, & Scatarzi (2004) and Stanghellini et al. (2006) all studies employed 
their own observational data to derive their abundances. Perinotto, Morbidelli, & Scatarzi (2004) 
and Stanghellini et al. (2006) adopted previously published line strengths determined by them to be 
of high quality. 


For each survey listed in Table 1, Fig. 5 shows the survey median (filled circles) for each observed 
element, presented in log space?. The solar abundance for each element is also plotted for comparison 
purposes. 


We see in the table that most of the survey values for any particular element are remarkably 
consistent with each other, regardless of the abundance spread found in each survey’. Furthermore, 
as a measure of the extent of agreement of survey median values for a particular element, the value 
of log(X + c) — log(X) ranges from a maximum in the case of S/H (0.13 or 35%) to a minimum in 
the case of He/H (0.02 or 596), where X and o represent the average and dispersion of the survey 
median values for an element, respectively. 

Potential reasons for differences among survey medians are numerous. They include uncertainties 
in: a) measured line strengths; b) reddening corrections; c) atomic data used to convert line strengths 
into ionic abundances; and d) ionization correction factors or photoionization modeling. Also, the 
presence of any radial, azimuthal, or vertical abundance gradients in the disk will likely cause survey 
medians to be affected by the distribution of sample PNe over the disk?. In addition, the surveys are 
arranged chronologically in Table 1 in order to check for systematic changes over time. Although this 
collection of disk surveys extends over approximately four decades in time, reassuringly no noticeable 
effects due to advances in instrumentation (e.g. the advent of CCDs in the early 1980s) or atomic 
data updates, are present. 


We can compare survey results with solar values by scanning down the column of medians for 
individual elements in Table 1. In doing so, we see that this group of surveys collectively seems 
to indicate the following: 1) He, C, N are enhanced by a few tenths of a dex, likely the result of 
nucleosynthetic processes in LIMS (see 82.1); 2) O, Ne, and Ar levels are close to solar; and 3) S and Cl 
tend to be subsolar — S abundances are often found to be depleted in PNe relative to values expected 


6 The standard abundance notation 12 + log X/H is employed, where X/H corresponds to the specific fraction at the 
bottom of a column along the horizontal axis. 

T This assessment does not include the published abundances by Barker (1978b) for N/H, Ne/H, and S/H, which seem 
to be consistently lower than the other samples regarding those same elements. 

8 Radial abundance gradients especially in the case of alpha elements (He burning products) such as O, Ne, S, Cl, and 
Ar are known to exist in the MW and most spiral galaxies. Typically, abundances of these elements decrease by a few 
hundredths of a dex per kiloparsec with galactocentric distance along the disk. See 84.1. 
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based upon their metallicity as measured by O/H (the "sulfur anomaly") for reasons that are not 
currently understood (Henry et al. 2012). In the same manner, Cl abundances are consistently found 
to be sub-solar in surveys of PNe, e.g., see Milingo et al. (2010, Fig. 3), although to our knowledge 
this behavior has not been investigated previously. 


3.1.2. MW Bulge PN Surveys 


Surveys of bulge PNe analyzed here are from Aller & Keyes (1987), Webster (1988), Ratag et al. 
(1997), Stasiriska et al. (1998), Exter et al. (2004), Wang & Liu (2007), Chiappini et al. (2009), and 
Cavichia et al. (2010, 2017). As in the case of the disk, we chose samples based upon homogeneity 
regarding data acquisition and abundance analysis. For all but the Cavichia et al. sample of south- 
ern PNe, objects not listed in the large compilation of northern bulge PNe found in Table A1 of 
Chiappini et al. (2009) were excluded from consideration. For the Cavichia sample, we considered 
only the PNe in their list which satisfied the criteria listed in Stasiriska et al. (1998). 


'The bottom portion of Table 1 shows the results for the eight surveys cited above, where again the 
surveys are ordered time wise beginning from oldest to most recent. To indicate the level of agree- 
ment among the surveys, the dispersion of median values around their average for the bulge surveys 
ranges from 0.03 dex (6.3%) for He/H to 0.19 dex (56%) for Ar/H. (Cl is discounted since only three 
of the surveys reported an abundance for this element.) Thus, for heavy elements there appears to 
be slightly more disagreement regarding their abundances among the bulge surveys compared with 
disk PN studies. This result is likely due to the larger distances and associated reddening effects for 
the bulge objects. And as in the disk, we see no evidence of temporal variations of survey medians 
for any element. 


A comparison of disk and bulge abundance results in Table 1 suggests a similarity between disk and 
bulge PNe in relation to solar abundances of several elements. For example, we see supersolar He/H 
and N/H but roughly solar O/H and Ar/H in both Galactic regions. At the same time, the Ne/H 
abundance is generally higher by about 30% in the disk than in the bulge, while CI/H abundance is 
roughly two times higher in the bulge. As in the disk, S/H in bulge PNe shows significant dispersion, 
with medians usually well below the solar level in both systems. Finally, as we saw for the disk, 
Fig. 6 demonstrates that bulge PN abundances qualitatively follow the solar pattern as expected. 


3.1.3. MW Halo Studies 


There are currently 13 PNe that are generally recognized as belonging to the halo of the MW. Four 
of these objects, BoBn1, DDDM1, H4-1, and K648, have been well-studied by five or more research 
teams each. The relevant studies are listed in Table 2 and stretch over three to four decades in time. 
For each object identified in the left column we list each paper by first author name and year which 
is keyed to a reference provided in the table footnote. Columns 3-10 contain number abundances 
relative to hydrogen for each element. Abundance results for the remaining nine halo PNe, each 
associated with one to four studies, are also featured in Table 2 below the first four. 


For each halo PN and element, Figure 7 shows the median value of all studies for that object. 
Two interesting conclusions can be drawn: 1) O, Ne, S, and Ar, alpha elements which reflect overall 
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metallicity, have median abundances below their respective solar values?; 2) He, C, and N exceed 
solar levels in a few of the halo PNe; and 3) the range of the median abundance values over the 13 
halo PNe is greatest for C (2.6 dex) and least for Ne (1.2 dex). 


We should note that for several PNe listed in Table 2 some uncertainty exists regarding halo mem- 
bership. Based on its radial velocity and location, BoBn1 may be a member of the Sagittarius dwarf 
spheroidal galaxy (Zijlstra et al. 2006). Further, the relatively small (absolute value) radial velocities 
of the following PNe may indicate thick disk membership: NGC2242 (Karachentsev, Karachentseva, 
& Huchtmeier 2004), NGC4361 (Torres-Peimbert, Peimbert, & Pena 1990), PNG243.8-37.1 (Boffin, 
Miszalski, & Jones 2012), PNG006.0-41.9 (Zijlstra et al. 2006); and possible bulge membership for 
M2-29 (Durand, Acker, & Zijlstra 1998). Improved radial velocity and distance determinations will 
be required for definitive classification. 


3.2. Additional MW Elemental Abundance Studies 


Smaller and/or more focused abundance studies published in the last two decades include Tsamis 
et al. (2003, 9 MW PNe), Costa, Uchida, & Maciel (2004, 26 MW anticenter PNe), Krabbe & Copetti 
(2006, 7 MW PNe), Yuan et al. (2011, 3D model of NGC 6153), Dufour et al. (2015, 10 MW PNe), 
Pagomenos et al. (2018, 23 MW anticenter PNe), and Miller et al. (2019, 6 MW disk PNe). There 
is also a series of 11 papers, each treating an individual PN, by Hyung et al. (see the last paper in 
the series, Hyung, Pottasch, & Feibelman (2004), which provides the references to the previous 10 


papers). 


Numerous projects providing abundance measurements for a variety of less studied elements are 
also available in the literature. Neutron capture elements, specifically those formed by the s-process, 
are one important example. Other elements which have received less attention are those residing on 
the third and fourth rows of the periodic table and include F, Na, Mg, AI, Si, P, K, Ca, Fe, Ni, and 
Zn. Below we summarize the current picture of s-process elements and then highlight abundance 
results of four particularly interesting elements: F, Mg, Fe, and Zn. These four elements are featured 
here because AGB stars are the only observationally confirmed production sites of F (Karakas & 
Lattanzio 2014, 85.2.3); as an alpha element, Mg tracks PN progenitor metallicity; and while the Fe 
abundance in PNe is very difficult to measure directly, Zn has been shown to track it. 


3.2.1. PN Abundance Studies of Neutron Capture Elements 


Many PNe are observed to have s-process elements such as Kr and Se present in measurable 
quantities. These post Fe-peak elements with atomic numbers Z > 30 inhabit the fourth row and 
beyond of the periodic table. They are produced when a nucleus (Z,A) first captures a neutron 
and subsequently undergoes a beta decay to yield the nucleus (Z--1,A 4-1). Thus a nucleus of a 
heavier element is produced, rather than a different isotope of the original element. In this instance, 
the process is termed s/ow because the local neutron flux is low enough that beta decay occurs 
before a second neutron is captured. The production of s-process elements in AGB stars is discussed 
in detail by Karakas & Lattanzio (2014, 83.7). In addition, readers can find a detailed review of 
neutron-capture element abundances in PNe in Sterling (2020a). For now we will briefly summarize 


? See the discussion in 84.1.3 regarding the evidence that O may actually be synthesized in low metallicity AGB stars 
following TDU and therefore is not a reliable indicator of the progenitor star's metallicity. 
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the present picture. 


Initial observations of emission lines related to s-process elements in PNe were carried out by 
Péquignot & Baluteau (1994), covering the spectral region of 6540-10460 À!?. They detected forbid- 
den lines of [Kr III] \6826.9 and [Xe IV] A7534.9 in NGC 7027 and proceeded to estimate the total 
abundances of Kr/H and Xe/H of 2.74 x 1078 and 3.36 x 107°, respectively. 


Table 3 provides a listing, from early to recent times, of authors, observed PNe, spectral coverage, 
and observed ions covering the time up to this writing. Note that the atomic number of each element 
is shown in parentheses following the element symbol at the top of the table. Papers listed in the 
first column are linked to their references in the table footnote. All transitions except those of Ba II 
are forbidden. Most of the papers listed report measurements of the [Kr III] and [Se IV] lines, with 
lines of Rb and Xe reported several times as well. Note the significant increase in the number of 
observed ions beginning with the paper by Sharpee, Zhang, & Williams (2007). For the most part, 
this acceleration in s-process work is due to the advent of high resolution optical and (especially) 
infrared instrumentation since that paper (Sterling 2020b). 


Many of the papers listed in Table 3 used their measured emission line strengths along with 
then-available transition energies and probabilities, and collision strengths to compute ionic and 
elemental abundances. Researchers often used ionization correction factors developed within the 
same publications through the use of detailed photoionization models. 


The resulting element abundances are listed in Table 4. Columns 1 and 2 provide the common 
object name and paper designation, respectively, where the latter is defined in a table footnote. 
When available, abundances reported in the papers are listed in the next eight columns. At the 
bottom of the table in a separate section, we list medians and maximum and minimum values of Se/H 
(in 68 objects) and Kr/H (in 39 objects) from a survey of 120 PNe by Sterling et al. (2015). The 
last line in the table provides the solar abundances by Asplund et al. (2009) for comparison purposes. 


Note that in the cases of Se/H, Kr/H, Rb/H, and Xe/H, many of the abundances tend to be above 
their solar values. For instance, recent studies of NGC 7027 by Sterling et al. (2016) and Madonna 
et al. (2018) indicate supersolar abundances of Se/H, Kr/H, and Rb. In addition, Sterling et al. 
(2015)'s large survey found a median value for the Kr abundance that is roughly four times the solar 
value, while the median Se/H level is very close to solar. Evidently, abundances of certain s-process 
elements tend to be enriched in many PNe. 


3.2.2. Fluorine 
The dominant fluorine isotope, F, may be partially produced in four steps summarized by 
MN -LEa-4Hpt— PF +2y+et +v in the He rich intershell of AGB stars. However, the bulk of °F 


in nature is produced in Type II supernova events when a neutrino (there are a lot of them produced 
during core collapse!) impacts a Ne nucleus, forcing out either a proton or neutron, ultimately 


1? Earlier NIR observations by Geballe, Burton, & Isaacman (1991) revealed two unidentified lines later associated with 
[Kr III] 2.199 jum and [Se IV] 2.187 um 
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resulting in the formation of a '?F nucleus (Clayton 2003). Zhang & Liu (2005) employed published 
strengths of the lines [F II] A4789 and [F IV] \4060 to determine F abundances in 13 Galactic PNe. 
Although their median value of F/H = 3.57 x 10^? is very close to Asplund et al. (2009)'s solar value 
of 3.63 x 1075, they nevertheless find a broad abundance range of 5.54 x 107? to 2.99 x 1077. In the 
case of NGC 40, F/Fo = 8.23, the largest in the survey, while the lowest value of that same ratio 
is 0.15, observed in NGC 2440. The authors also suggest that a positive correlation exists in PNe 
between [F/O] and C/O. Otsuka et al. (2008) computed the F abundance in the low-metallicity halo 
PN BoBn-1, using their emission measurements of the [F IV] 43997 and A4060 lines. They found 
that F/H= 3.31 x 1077, roughly an order of magnitude above the solar value and the median value 
observed by Zhang & Liu (2005). Otsuka et al. (2008) invoke a picture similar to that of Zhang & 
Liu (2005) regarding the origin of F enhancement in BoBn-1. Another low metallicity PN, J900, was 
observed by Otsuka & Hyung (2020), who measured a value of F/H=1.45 x 1077 or roughly four 
times the solar F/H value. 


A comprehensive understanding of F in PNe has begun to emerge with the work of Abia et al. (2019), 
who observed 11 AGB carbon stars, five SC-type and six N-type!!. They found that at roughly solar 
metallicity, SC and N-type stars possessed nearly equal F abundances. In the same study, Abia et 
al. (2019) found an inverse correlation between [F/Fe] and [Fe/H] for —2 < [Fe/H] < 0 (see their 
Fig. 3). The authors also showed by plotting theoretical models from Cristallo et al. (2015) that at 
constant [Fe/H], [F/Fe] directly increases with the number of thermal pulses. They conclude from 
their studies, however, that AGB stars are not the main contributors of F in the solar neighborhood 
(core collapse SNe are) in agreement with chemical evolution models extant in the literature, e.g., 
Kobayashi et al. (2011). Further analysis by this same team presented in Vescovi et al. (2021) also 
finds that F abundance in AGB stars is directly related to the average s-process abundance <s> such 
that [F/Fe] and |<s>/Fe] increase together (see their Fig. 2, upper panel). Overall, they also find 
that their stellar models better reproduce the F/H abundance observations when magnetic buoyancy 
*acts as the driving force for the formation of the !?C neutron source" in metal-poor AGB stars 
during TDU. 


3.2.3. Magnesium 


As an alpha element, the Mg found in a PN was likely present in the progenitor star at the time of 
its formation. The dominant Mg isotope is "Mg, which represents about 79% of all Mg in the Sun 
(Asplund et al. 2009). 24Mg is a product of He burning, being produced mainly by the ?Ne(o, ,)?* Mg 
reaction (Clayton 2003). Mg is especially interesting because of its moderate tendency to form dust, 
with a condensation temperature of 1346 K (Lodders 2003, fosterite formation) and a general de- 
pletion factor? that ranges from —0.270 to —1.267 in the ISM, depending upon the line of sight 
(Jenkins 2009). These values imply percent depletions of 4696 to 9596. 


The first observations of Mg emission in a PN were made by Grewing et al (1978), who measured 
fluxes of Mg II AA2795, 2803 and [Mg V] 2783 in NGC 7027 using the IUE. A photoionization model 
analysis of NGC 7027 by Péquignot, Aldrovandi, & Stasitiska (1978), using published line strengths 


11 The spectral evolution of AGB stars as the number of TDU episodes increases is M-MS-S-SC-C(N), where C(N) is 
N-type. 
1? The depletion factor of an element X is defined as [Xgas/H] = log(X/H) gas — log(X/H)o. 
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for a range of elements including Mg from a number of sources, found Mg/H=3.5 x 10^? (—0.06) ^, 
a value very close to Asplund et al. (2009)'s solar value of 3.98 x 10 ?. Followup work on NGC 7027 
by Péquignot & Stasiriska (1980) found a Mg negative radial gradient across the nebula such that 
the gaseous Mg abundance in the inner nebula exceeded that of the outer nebula abundance. The 
authors proposed that Mg depletion into dust was more likely in the cooler outer nebular regions 
where Mg II emission dominated, while Mg precipitated less in the hotter region closer to the central 
star dominated by the [Mg V] emission. Subsequent studies of Mg in PNe included work on IC 418 
by Harrington et al. (1980), Middlemass (1988), and Hyung, Aller, & Feibelman (1994) who found 
Mg/H- 2.0 x 10~°(—0.30), 4.0 x 10^? (+0.002), and 6.91 x 10-9 (—0.76), respectively. Likewise, 
Mg/H in IC 4997 was measured by Middlemass (1988) to be 9.12 x 10-9 (—0.64), while Hyung & 
Aller (1998) found a subsolar value of 2.0 x 10 9 (—1.30) for NGC 2440. Finally, Wang & Liu (2007) 
observed Mg II M4S81 in six MW bulge PNe and found the average Mg/H=4.79 x 10^? (+0.08), close 
to the solar value of 3.98 x 10^? (Asplund et al. 2009). Observationally, then, Mg depletion in the 
PNe appears to range from 096 in NGC 7027 to 9596 in NGC 2440. 


3.2.4. Iron 


According to Jenkins (2009) and Dwek (2016), the Fe depletion factor in the ISM has a range of 
—0.951 to —2.236, corresponding to percentage depletions of 89% to 99%, depending upon the line 
of sight. The Fe observed in PNe was present in the interstellar material at the time of the formation 
of the progenitor star and exists in the nebula in both gas and dust form. Roughly 92% of Fe is °°Fe, 
which results from £* decay of ?9Ni via ?9Co to ??Fe. The original °°Ni is produced by explosive 
Si-burning in single massive stars (Type II supernovae) or C ignition in a C-O white dwarf member 
of a binary stellar system (Type Ia supernovae). 


'The Fe abundance is difficult to measure in emission nebulae because of its weak lines due to heavy 
depletion onto dust. It was first measured in a PN by Shields (1975), who used previously published 
optical line strength data to determine Fe/H in NGC 7027 to be 7.94 x 1077. The corresponding solar 
value is 3.16 x 10^? (Asplund et al. 2009), implying an Fe depletion factor of —1.60 or a percentage 
depletion of 9796. 


In a followup study, again using published data, Shields (1978) repeated measurements of Fe/H 
in NGC 7027 and included five additional PNe, NGC 1535, NGC 2022, NGC 6741, NGC 6886, and 
IC 2165. In the same year Garstang, Robb, & Rountree (1978) used optical line strengths published 
by Kaler (1976) to determine Fe/H in three PNe, NGC 6720, NGC 7009, and NGC 7662. Finally, 
Perinotto et al. (1999) used their own observed optical line strengths to measure Fe/H values and Fe 
depletion factors for NGC 7027, NGC 6543, Hu 2-1, and Cn 3-1. 


'The numerical results for these three studies are presented in Table 5. Note that except for the 
case of NGC 7009, each PN exhibited an Fe depletion of more than an order of magnitude. In the 
case of NGC 7027, which was measured twice by Shields and once by Perinotto, the depletions range 
from —1.51 to —1.90 with an average of —1.70. Marked Fe depletion is seen consistently in these 11 


15 The depletion factors corresponding to the Mg abundances quoted in this paragraph are given in parentheses. 
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PNe. 


More recently, Delgado-Inglada et al. (2009) measured Fe in 33 MW PNe. They found a range in 
Fe/H of 3.09 x 107? to 1.86 x 1078, a median value for Fe/H of 7.08 x 1077, and a median depletion 
value of —1.65 or 98%. 


Fe/H is an important quantity to determine in nebulae, as it can be compared with stellar Fe/H 
measurements in order to contrast metallicities in these two environments'*. Alas, its marked deple- 
tion makes this task difficult. But read on! 


3.2.5. Zinc 


As just discussed, estimating the Fe abundance in PNe is difficult because of its large deple- 
tion factor. Zinc, however, is more volatile and has a depletion factor range of —0.059 to —0.551 
(Jenkins 2009), corresponding to depletions of 13% to 72%, respectively. In addition, Zn/Fe is 
found to maintain very nearly its solar value of 1.15 x 10^? (Asplund et al. 2009) over the range of 
—2 < [Fe/H] € +0.5 (Saito, Honda, & Takeda 2009, Fig. 3). Since stellar metallicities are usually 
represented by the Fe abundance, knowing the Zn/Fe ratio and measuring the Zn abundance in PNe 
allows a direct comparison to be made with the former. 


The stable isotopes of zinc are ^Zn, 99Zn, 9"Zn, and 99Zn, where the first two represent nearly 77% 
of the total element abundance (Asplund et al. 2009). These two isotopes can be formed via alpha 
freezeout following a core collapse supernova event, explosive Si burning, or the s-process in massive 
stars (Clayton 2003). However, 9^9975 are also at the beginning of the chain of s-process elements 
and can be produced in observable quantities in AGB stars (Karakas & Lugaro 2010). 


Dinerstein (2001) first identified [Zn IV] 3.625 um as an emission feature in two Galactic disk PNe, 
NGC 7027 and IC 5117. In the absence of any relevant Zn collision strengths, accurate abundances 
could not be determined. However, the authors concluded that, for reasonable assumptions about 
values of collision strengths based on similar transitions in other elements, the Zn abundance in these 
two objects is not enhanced. 


Smith, Zijlstra, & Dinerstein (2014) measured the strength of the [Zn IV] 3.625 um in nine PNe, five 
of which are located in the Galactic bulge. With a collision strength now available, they computed 
Zn abundances and concluded that Zn/H is subsolar in a majority of their sample objects. Three 
years later Smith et al. (2017) observed six additional bulge PNe and updated their observations of 
the nine objects observed in Smith, Zijlstra, & Dinerstein (2014). Values of Zn/H abundance ratios 
were determined for all 15 objects. The median for both the disk and bulge samples was found to 
be Zn/H = 2.00 x 1078, a bit less than the solar value of 3.63 x 10 *(Asplund et al. 2009). They 
concluded that while Zn is generally subsolar in both the disk and bulge PNe, their ratios of O/Zn 
were above solar and indicated that O is enriched relative to Zn (Smith et al. 2017, Fig. 2). More 
generally, if Zn abundances track those of Fe, the inverse correlation of [O/Zn] vs. [Zn/H] shown in 
Smith et al.’s Fig. 2 is quite consistent with the well established decrease of [O/Fe] with increasing 


14 Usually nebular metallicity is determined using O/H, while in stars, Fe/H serves the same purpose. 
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[Fe/H] that is seen in the solar neighborhood and elsewhere in the Galaxy. Kobayashi, Karakas, & 
Lugaro (2020, Fig. 5) recently compared chemical evolution model results with observations to show 
this behavior between [O/Fe] and [Fe/H], which is likely due to the delayed contribution of Fe-peak 
elements to the Galactic environment by SNIa. The striking similarity between the [|O/Zn] vs. [Zn/H] 
and [O/Fe] vs. [Fe/H] further supports the assumption that Zn and Fe abundances change together 
in lockstep. The next step in research would seem to be to compare Fe abundances in PNe, inferred 
from the constant Zn/Fe ratio, with Fe measurements in central star atmospheres where available. 


3.3. Extragalactic Planetary Nebula Elemental Abundance Studies 


So far we have focused on chemical abundance work pertaining to the MW. Analogous studies 
have also been carried out on at least 19 external galaxies. Just within the Local Group, abundance 
characteristics of four galaxies, M31, M33, and the Large and Small Magellanic Clouds have each 
been studied by numerous research teams. Other Local Group members receiving less attention 
include the spheroidals Fornax, M32, NGC 147, NGC185, NGC 205, and dwarf spheroid Sagittarius 
dSph; a barred Magellanic spiral NGC 3109; barred Magellanic irregulars NGC 4449 and NGC 6822; 
and Magellanic irregulars IC 10, Sextans A, and Sextans B. Finally, abundance work on galaxies 
outside of the Local Group includes the spirals M81 and NGC 300, as well as the lenticular galaxy 
and strong radio emitter NGC 5128 (Centaurus A). 


Tables 6A and 6B provide a summary of the important parameters of each abundance study. For 
each paper listed in the table, the data were both collected and analyzed by the authors themselves. 
The bold-print name of each galaxy appearing in the center column is followed in square brackets by 
the NED classification designation of galaxy type along with the associated T-type index number, 
with morphology extending from elliptical through spiral to irregular galaxies as T values increase. 
Columns 1-5 list the first author name and year in a contracted format (see the table footnote 
for the associated reference), the number of PNe with abundances reported in that paper, the ele- 
ments whose abundances were determined, the spectral range(s) of the observations from which line 
strengths were determined, and the median (23 objects), average (2 objects) or single (1 object) 
value of O/H in the objects observed. T'hus,O/H* is a statistical value derived from O/H values for 
individual PNe within a single sample. Finally, Column 6 provides comments where appropriate. As 
done previously, papers for each galaxy are listed in order of earliest to latest. 


We note that values in Column 5 of Tables 6A and 6B are meant to provide only a rough gauge of 
the metallicity levels of the PN samples observed, using O/H* as the metallicity indicator. Differ- 
ences in O/H* among studies are partially related to the measurement uncertainties associated with 
individual surveys. Additionally, O/H variations among objects within a single survey may also be 
due to the range in progenitor masses of the PNe comprising the survey. Lastly, in the case of galaxy 
disk samples, the details of how the PNe are spread out radially, vertically, and azimuthally in the 
potential presence of gradients or any nonuniformity in element distributions over the disk will also 
impact the value of O/H*. 


Figure 8 is a plot of the logarithmic values of O/H* in Tables 6A and 6B versus galactic mor- 
phological T-type. Individual galaxies listed in the tables have corresponding color coded symbols 
defined in the legend, while each O/H* value listed in the tables for a particular galaxy is shown. For 
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reference purposes, the crosses at T—4 show the log(O/H*) values from Table 1 for the MW disk. We 
note the trend toward lower metallicity with increasing T-type between T—3 (the spiral galaxy M31) 
and T=10 (irregular and dSph galaxies), although there are too few galaxies with T<3 to determine 
if this trend extends to earlier objects. This apparent indirect correlation between metallicity and 
galaxy type is likely related to the complex interactions of total stellar mass, metallicity, and star 
formation history uniquely associated with each galaxy. 


3.4. Molecules in Planetary Nebulae 


Molecules in both the gas and solid phases represent an important component of both PPNe and 
PNe. As an example, in his recent review of molecules in PNe, Zhang (2017) points out that more 
than 80 different molecules have been observed in C-rich AGB envelopes and PPNe. Zhang also 
lists nearly 30 molecules that have been discovered in PNe alone (see his Table 1). These molecules 
range from the simple and abundant two-atom species such as CO, OH, and H» to the complex 
fullerene molecules, Cg) and Cz. In addition, the quantitative measurement of molecular mass in 
the circumstellar envelopes related to PNe likely explains at least some of the difference between the 
combined masses of the central star and ionized nebula and the presumed mass of the main sequence 
progenitor star, the so-called *missing mass" (Kwok 1994; Zhang 2017). 


Because of the significant mass loss that occurs during the AGB phase, a circumstellar envelope 
forms and expands radially around the star. This envelope material is cool enough and of sufficiently 
high density to permit stable chemical bonding to occur between two or more atomic species. As a 
result, molecules presently observed in PNe are primarily the result of the chemical processing that 
occurs beginning with the AGB mass loss phase, continuing through the intermediate PPN stage, 
but diminishing when the central star becomes hot enough to photoionize the circumstellar material 
to form the PN with its hotter environment. The post-AGB chemistry may also be affected by the 
presence of an external ambient UV radiation field. Brief reviews of the physics and chemistry of 
neutral gas in and around PNe are provided by Dinerstein (1991), Hollenbach & Tielens (1997, 84.6), 
Kwok (2000, Ch. 5), Van Winckel (2003), Hasegawa (2003), Bujarrabal (2006), and Zhang (2017). 
And for readers in need of a refresher course in the astrochemistry that is relevant to what follows, 
Kwok (2007) is a good source of information concerning molecular orbitals, polycyclic aromatic 
hydrocarbons, and organic compounds [see 8887.3, 8.5, and 11.3, respectively in Kwok (2007)], while 
Herbst & van Dishoeck (2009) present a detailed discussion of complex organic molecules in the 
interstellar medium. 


In the following discussion, we begin by focusing on the first two molecules to be observed in PNe, 
CO and Hə. After that, we briefly explore many of the more complex molecules, evidence of which 
has been observed more recently. 


3.4.1. CO 


The first detection of a molecule in a PN was made by Mufson, Lyon, & Marionni (1975), who 
detected 2.6mm emission from the J—1-0 transition of CO in NGC 7027, IC 418, and NGC 6543. 
The authors made note of the large volume of gas containing CO and extending beyond the smaller 
ionized region in NGC 7027. From this they speculated that the molecular gas was expelled by the 
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PPN central star as the PN was forming. 


In a set of three papers, the team of P.J. Huggins and A.P. Healy measured the more energetic 
J=2-1 CO transition at 1.3mm in NGC 7293 (Helix) (Huggins & Healy 1986a), NGC 2346 (Hug- 
gins & Healy 1986b; Healy & Huggins 1988) and NGC 6720 (Ring) (Huggins & Healy 1986b) to 
study the distribution and kinematics of CO in these objects. For NGC 7293, the authors con- 
cluded that a significant fraction of the matter lost by the stellar wind during the PPN stage has 
not been ionized by the central star. Thus, roughly 10% of the system’s total mass remains in 
the neutral state. Likewise, relative to its total mass NGC 2346 has an estimated neutral mass of 
greater than 40%, while the estimated neutral mass associated with NGC 6720 is 7% of its total mass. 


Huggins & Healy (1989) carried out a survey of 100 Galactic PNe again at 1.3mm and detected CO 
in 19 of the objects. The investigators noted that these PNe, when compared to the 81 objects with 
no CO detections, appeared to have had high-mass progenitors, based upon each PN’s relatively high 
N/O ratio and strong preference for bipolar morphology. The 19 objects also exhibited a range in 
mass of the molecular gas of 107? to 1 Mo. In a second survey study, Huggins et al. (1996) observed 
91 Galactic PNe, 23 of which were found to harbor CO (three of which they had studied previously). 


Table 3 in Huggins et al. (1996) summarizes the observations of 44 PNe with positive CO detec- 
tions, where most of the objects were part of the two surveys cited above while data for an additional 
few objects were taken from the literature. The table lists, among other characteristics, the radius 
of the ionized nebula, N/O, PN morphological type, total molecular mass M,, (inferred from the 
CO mass), and the ratio of total molecular to ionized nebular mass (M,,/M;). For the 44 PNe, the 
median value of Mm is 0.031 Mo, with minimum and maximum values of 4 x 10 *M (J900) and 
0.79Me (M1-16), respectively. Values for M,,/M; are 0.47 (median) and 5.4 x 107? (NGC 7008; 
minimum), and 490 (CRL618; maximum). Thus, for PNe in which CO is detected, there is a large 
range in Mm as well as M,,/ Mi. 


The most interesting finding in Huggins et al. (1996) is illustrated in their Fig. 11, where a loga- 
rithmic plot of the molecular to ionized mass, M, /Mi, versus the radius of the ionized nebula (PN) 
clearly shows an inverse relation for the objects in which CO is detected. The authors attribute 
this correlation to the growth in size of the ionized nebula as ionizing photons from the central star 
cause the CO to disassociate and atoms to become ionized. In their view, this result supports the 
notion that the matter comprising the PN was initially neutral and molecular in composition but 
became dissociated and ionized by photons as the central star moved horizontally left on the HR 
diagram as its temperature rose. The ionization front moves outward radially at the expense of 
the molecular gas. Thus, the inverse correlation between M,,,/M; and nebular radius represents an 
evolutionary sequence. As for those sample PNe in which CO is undetected, the authors speculate 
that the molecular gas became rapidly photo-dissociated either before or during the formation of the 
visible PN. 


Finally, a few recent examples of CO observations of PNe can be found in Doan et al. (2017), 
Guzman-Ramirez et al. (2018), and Andriantsaralaza, Zijlstra, & Avison (2020). Doan et al. (2017) 
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observed the AGB star c! Gruis, confirming the system's torus-bipolar structure as well as the 
mass, velocity, and momentum of the outflow. Guzman-Ramirez et al. (2018) used ESO's Atacama 
Pathfinder Experiment (APEX) radio telescope facilities to survey 93 Galactic PNe. They detected 
CO emission in 21 of the objects and measured the column density and molecular mass in each of 
them. Andriantsaralaza, Zijlstra, & Avison (2020) used ALMA to observe three isotopic forms of 
CO, ?CO, PCO, and !5CO in the C1 knot of the Helix Nebula, NGC 7293. They determined the 
total molecular mass of the globule to be z 2 x 1074 Mọ (much greater than previously published 
values) and the mass of the progenitor star to be 2 Mo. 


3.4.2. Ho 


Shortly after the initial observation of CO, H, was first observed at 2 um in NGC 7027 by Treffers 
et al. (1976) and followed by Hə observations by Beckwith, Persson, & Gatley (1978) in five addi- 
tional PNe. Webster et al. (1988) later observed Hə in 11 PNe, and concluded that PNe with strong, 
excited Hy emission are Type I PNe (see 84.2.2) which also display a toroidal morphology with faint 
bipolar extensions. Hə was observed in three PNe by Zuckerman & Gatley (1988), whose conclusions 
supported those of Webster et al. (1988) regarding morphology and spectral type. Further solid 
reinforcement of the link among H» emission, bipolar morphology, and the presence of a toroidal disk 
came from the extensive study by Kastner et al. (1996), who searched roughly 60 PNe and found 23 
objects with the 2.122 um emission feature. These same objects also were observed to be located 
at low galactic latitude. The authors concluded that the progenitors of the central stars were likely 
relatively massive. Kastner et al. (1996) also speculated that instead of fluorescence, shocks heated 
the gas to high temperatures and thus colllisional excitation was more likely to be the excitation 
mechanism of the 2.122 um emission. Support for this idea was provided by Guerrero et al. (2000), 
who presented Hə and Bry narrow band imaging in 15 bipolar PNe, finding a strong correlation 
between shock-excitation and bipolar morphology. 


Next, García-Hernández et al. (2002) for the first time extended H, studies to PN precursors, i.e., 
late AGBs, PPNe, and young PNe, finding results similar to those of Kastner et al. (1996) and 
Guerrero et al. (2000) at earlier stages, when Hg is shock excited rather than excited by fluorescence. 
Kelly & Hrivnak (2005) followed up with a study of 51 PPNe, 16 of which were observed to exhibit 
Hə emission. They found that Hz emission is predominantly associated with bipolar nebulae located 
at low galactic latitudes. Rosado & Arias (2003) inferred the total mass of H» in five PNe, NGC 6720, 
NGC 6781, NGC 3132, NGC 2346, and NGC 7048, from the observed mass of shocked H5 gas and 
knowledge of initial density obtained from measured 2.124 m surface brightness and shock velocities 
in each object. This was accomplished in each case by establishing the pre-shock density of the 
molecular torus from the observed 2.12um surface brightness and shock velocity. In this way, Rosado 
& Arias (2003) found a range of 0.1 to 0.8 Mo in total Hə mass. Their study also implied that 
the major source of excitation of the 2.122 and 2.248 um lines is shock heating. Support for this 
conclusion was later provided by the observations by Likkel et al. (2006) of H5 emission in eight 
PNe, where a few of their objects displayed H» line ratios characteristic of shock excitation. 


More recent work on H5 in PNe has helped to refine our view of the location and nature of the 
molecular hydrogen observed in PNe. Examples of studies that have been carried out along this line 
include those by Marquez-Lugo et al. (2013), Manchado et al. (2015), Akras, Gongalves, & Ramos- 
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Larios (2017), and Akras et al. (2020). Using 8um images of Hp emission in PNe, Marquez-Lugo et al. 
(2013) showed that such emission is not produced exclusively in bipolar PNe, as had been previously 
assumed, but can also be associated with ellipsoidal/barrel PNe. Manchado et al. (2015) used highly 
resolved H5 images in the IR of the PN NGC 2346 to show that the H5 gas is composed of knots 
and filaments. They speculated that thermal pressure of the ionized gas resulted in fragmentation 
of the swept-up shell as the system expands. Akras, Gonçalves, & Ramos-Larios (2017) employed 
high resolution images in Hə 2.122 and 2.248um of two PNe, K 4-47 and NGC 7662, to study low- 
ionization knots (internal structures in PNe are discussed in 87.2.1). Their results confirmed the 
presence of Hə gas in both fast and slow-moving knots and suggested that low-ionization structures 
in PNe are characteristically similar to photodissociation regions. In their most recent paper on low- 
ionization structures, Akras et al. (2020) again employed 2um observations of H5 to study NGC 7009 
and NGC 6543. In comparing their own observations and similar ones found in the literature with 
extant shock and UV excitation models of Hz emission, the authors were unable to decide definitively 
which mechanism is most relevant. 


3.4.3. Additional Molecules 


Evolution of the molecular composition of PNe was the subject of a paper by Bachiller et al. (1997), 
who found that molecular abundances in the circumstellar envelope are altered considerably between 
the AGB and PN stages. The same study indicated that as PNe themselves evolve, molecular 
abundances remain relatively constant. This latter point was supported by Edwards, Cox, & Ziurys 
(2014), who measured abundances of CO, CS, and HCO* in five PNe spanning the ages of 900-10,000 
years and found the abundances of each of these molecules tend to be largely constant over time. An 
additional correlation with implications for linking molecular chemistry with central star properties 
is the discovery by Bublitz et al. (2019, their Fig. 15) that the HNC/HCN abundance ratio decreases 
as the UV luminosity of the central star increases, in a power law relation characterized by a slope 
of -0.363 and a correlation coefficient of -0.885. 


In a recent conference review paper, Schmidt & Ziurys (2018) presented results from their decade- 
long work to detect PN emission of numerous molecules besides CO and Hg. Detections of one or 
more simple organic compounds were made in 14 of the 17 objects observed. They concluded that: 
1) polyatomic molecules are observed in PNe of various ages and morphologies; 2) abundances of the 
molecules measured appear not to vary with PN age (see their Fig. 2); and 3) molecules in PNe are 
located in dense clumps where photodissociation is minimal due to shielding by dust. 


Two of the most interesting molecules found in PNe are Ceo and Czo, the largest and most stable 
molecules ever detected in space. The first observation of these species was made by Cami et al. 
(2010), who employed the IRS aboard the Spitzer Space Telescope in a study of the PN Tc 1, where 
the abundance ratio C/O#3.2. The authors observed what appeared to be a complete absence 
of H-containing aliphatic (open chain) molecules such as HCN and C2H2, as well as polycyclic 
aromatic hydrocarbons (PAH), i.e., molecules built around one or more benzene rings. From this 
they suggested that fullerene formation requires an H-free environment, a rather rare situation in PNe. 


This picture was challenged by García-Hernández et al. (2010), who showed that fullerenes are 
observed in PNe with normal H abundances together with other H-rich species like PAHs, and sug- 
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gested the top-down photochemical processing of hydrogenated amorphous (no crystalline structure) 
carbon (HAC) dust grains for fullerene formation. These results were confirmed by IR observations of 
more numerous PNe samples [García- Hernández et al. (2011, 2012a), Otsuka et al. (2014)] and other 
types of objects (see e.g., García-Hernández (2012b) for a review). For example, in García-Hernández 
et al. (2012a) the authors compiled a total of 263 PNe observed by Spitzer/IRS. Out of this sample, 
the 16 PNe found to display fullerene spectral features were analyzed in order to discover the kinds 
of environments which are favorable for the formation of these molecules and from what materials 
do they form. They concluded that fullerene detection in PNe increases with decreasing metallicity 
and that fullerenes (along with PAH and aliphatic species) most likely form from a carbonaceous 
compound that is a mixture of aromatic and aliphatic structures (HAC-like), where the chemical 
processing is driven by UV radiation from the hot central star and/or post-AGB shocks (see also 
Micelotta et al. 2012). This finding is in agreement with predictions based upon laboratory results 
published earlier by Scott, Duley, & Pinho (1997). Interestingly, the idea of HAC-like starting 
materials would be consistent with the efficient formation of non-aromatic carbonaceous molecules 
in circumstellar environments as evidenced by recent top-level experiments (Martínez et al. 2020). 
However, other starting materials for top-down fullerene formation have been proposed in the litera- 
ture; e.g., the photochemical processing of large PAHs (Berné & Tielens 2012) or shock heating and 
ion bombardment of SiC grains (Bernal et al. 2019), among others. Cao and C; fullerenes may be 
just the tip of the iceberg, and many other fullerene derivatives (see e.g., Omont 2016) are potentially 
present in PNe environments. 


Finally, in a more speculative vein, the presence of fullerenes in PNe suggest that molecules sus- 
pected of being precursors of Cgo;9 might also be observable. One such candidate is graphene, 
Cy4. [An example of a proposed mechanism for forming fullerene from graphene was suggested by 
Chuvilin et al. (2010).] García-Hernández et al. (2011, 2012a) observed IR emission features located 
at 6.6, 9.8, and 20um in five PNe in the Magellanic Clouds as well as in the Galactic PN K 3-24. 
These three transitions correspond to the strongest ones associated with C5, that were previously 
predicted from theory by Kuzmin & Duley (2011), implying the presence of graphene in these six 
PNe. Recently, Li, Li, & Jiang (2019) have inferred an upper limit of graphene abundance in the 
ISM of 20 ppm, based upon the absence of a 2755 À graphene absorption feature in the interstellar 
extinction curve. Thus, the presence of graphene in PNe and the ISM remains uncertain. 


So for the final exam, what are the main points to know about molecules in PNe? They are: 
1) PNe with detectable CO and Hz appear to have relatively massive progenitors and tend to possess 
a bipolar morphology; 2) the ratio of molecular to ionized gas mass is inversely proportional to the 
radius of the ionized region (the PN); 3) although the abundances of a variety of molecules relative 
to H5 change with respect to each other during the PPN stage, they tend to remain constant after 
PN formation; and 4) fullerenes in PNe form top-down through a still unknown starting material 
and a rich family of fullerene derivatives is likely present. 


3.5. Dust in Planetary Nebulae 


Early recognition of the possibility that PNe harbor dust came from IR photometric studies by 
Gillette, Low, & Stein (1967); Gillette, Merrill, & Stein (1972), Persson & Frogel (1973), and Willner 
(1972) that showed the presence of IR excess beyond what was estimated to have come from atomic 
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processes. Several years later, Moseley (1980) observed 13 PNe in the far-IR (37-108 um) and 
concluded that all but one object (IC 418) have roughly the same dust masses. Natta & Panagia 
(1981) then employed the IR data from Moseley (1980) to show that the dust/gas ratio, grain size, 
and internal optical depth decreased with increasing values of nebular radius relative to that of 
NGC 7027. The authors speculated that these trends were the result of: 1) an inhomogeneous 
distribution of nebular dust resulting from an increase in the concentration of dust as the AGB 
stellar mass loss progressed; or 2) the partial destruction of grains over time. Finally, an estimation 
of the total PN contribution of interstellar dust to the Galactic ISM of 10 9Myr ! was made by 
Gehrz (1989). Such were the results of the early explorations of dust in PNe. More general discus- 
sions and reviews of dust in PNe may be found in Kwok (2000, Chapter 6), Kwok (2007, Chapters 
10-13), Gail & Sedlmayr (2014), Kwitter et al. (2014, 89), and García-Hernández (2015). For the 
following discussion, we employ the notation for the various dust types used by García- Hernández & 
Górny (2014), namely, Featureless (F), C-rich (CC), O-rich (OC) and mixed dust or dual-dust (DC). 


Dust studies of relatively large samples of PNe have been carried out by Stanghellini et al. (2007, 
2012), Bernard-Salas et al. (2009), García-Hernández & Górny (2014), and García-Hernández (2015). 
Stanghellini et al. (2007) used the Spitzer Space Telescope with its IR spectrograph between 5-40 jum 
to observe 41 PNe, 25 located in the LMC and 16 located in the SMC. Twenty of their sample 
objects (9 in the LMC and 11 in the SMC) possessed broad CC dust features, three displayed 
broad OC features, and the remaining 18 PNe lacked broad dust features altogether. For each 
PN, the authors then compared the inferred dust type with its previously published C/O value 
as determined from collisionally excited lines. They found that C/O>1 for objects with CC dust 
and C/O<1 for objects with OC dust. The latter case was attributed to the effects of HBB, with 
its associated conversion of C to N via the CNO cycle. The authors also found that their CC 
PNe were all morphologically symmetric, while the OC objects appeared highly asymmetric with a 
tendency toward harboring more massive central stars. In a later and much larger dust study of 150 
Galactic bulge and disk PNe, Stanghellini et al. (2012) found the relative population by dust type 
to be 1796 F, 2596 CC, 3096 OC, and 2896 DC. Furthermore, the CC group was found to comprise 
two molecular types, i.e., aromatic and aliphatic, while the OC group was split into two structural 
groups, i.e., crystalline and amorphous.They found that: 1) DC PNe favor locations closer to the 
Galactic center; 2) aliphatic chemistry was favored in CC PNe, while amorphous structure was 
preferred in OC PNe; 3) Nebular radii of CC PNe scaled inversely with both dust temperature and 
IR luminosity; and 4) CC PNe significantly outnumber OC PNe in the Magellanic Clouds in contrast 
to their relative populations in the Galaxy. This last point is likely the result of both the greater 
efficiency of TDU and the heightened effectiveness of increasing the C/O ratio at the low metallicities. 


Bernard-Salas et al. (2009) performed a low resolution spectroscopic study of dust features of 18 
LMC and 7 SMC PNe using the Spitzer Space Telescope. PAH emission features, appearing at 6.2, 
7.7, 8.6, and 11.2 um were seen in 14 of the PNeP, while nine PNe displayed broad SiC and MgS 


15 C-rich PAH molecules in gaseous form are often present in CC environments, and the former are usually taken as 
evidence of the presence of the latter. However, Cohen & Barlow (2005) observed the 7.7 um feature of PAH in 
NGC 6302 and Hu 2-1, two objects in which C/O was found to be 0.31 and 0.48, respectively. In addition, Guzman- 
Ramirez et al. (2011) studied 40 Galactic bulge PNe and found DC dust chemistry and evidence of PAH in 30 objects. 
'They concluded that the DC chemistry is related to hydrocarbon chemistry occurring in the presence of UV radiation 


within dense tori. Therefore, the link between PAH and CC dust chemistry is not absolute. 
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emission at 11 and 25 um, respectively. All of these features are indicative of a CC environment. 
At the same time, the authors detected amorphous silicate emission, usually associated with an OC 
environment, in only two of their sample objects. Again, as found by Stanghellini et al. (2012), the 
CC environment appears to be favored at the low metallicity that is characteristic of the Magellanic 
Clouds. 


García-Hernández & Górny (2014) presented a sample of 131 Galactic bulge and disk PNe for 
which both optical and IR spectroscopic data were available. The authors used Spitzer IR spectral 
data from Gutenkunst et al. (2008), Perea-Calderón et al. (2009), and Stanghellini et al. (2012), 
coupled with optical spectrophotometric data of their own or critically chosen from the literature, 
to determine dust type and gas phase chemical abundances for each PN in order to explore how the 
two properties are connected. 


Fig. 9 is a column plot of data contained in Table 5 of García-Hernández & Górny (2014) showing 
the percentage of the total number of bulge or disk PNe that belong to the indicated dust types. We 
can see that while the bulge and disk exhibit roughly the same proportions of F and OC PNe, the 
disk has a significantly greater proportion of CC PNe and a much smaller proportion of DC PNe. 
In fact, roughly half of the bulge PNe belong in the DC category. For further comparison, we have 
added data from Table 6 of Stanghellini et al. (2012) for the Magellanic Clouds. In this case we can 
see that roughly half of the PN sample belong to the CC group, two-fifths to the F group, and only 
a relatively small number in OC group. 


Values from García-Hernández & Górny (2014) listed here in Table 7 link the median gas phase 
abundances of He, N/O, and Ar/H to the four dust types for both the Galactic bulge and disk PNe, 
where Ar/H is taken by the authors as the metallicity indicator. For comparison, solar abundances 
from Asplund et al. (2009) are provided in the last row of the table. The values in Table 7 indicate 
clearly that in both the bulge and disk: 1) the He/H ratio in PNe of all dust types is about 25% 
greater than the solar value; 2) the ratio of N/O is significantly greater in DC PNe than in the other 
three dust types; and 3) the metallicities of DC PNe are markedly greater than the levels seen in 
both CC and OC PNe. Regarding alpha elements other than Ar, Figs. 8 and 9 in García-Hernández 
& Górny (2014) show that in both the bulge and disk, Ne and S abundances in OC and DC objects 
appear to increase in lockstep with O as is observed in HII regions. But curiously, in the case of 
CC PNe, the S abundances are consistently below what is expected from their O abundances. This 
is perhaps related to the sulfur anomaly described in (Henry et al. 2012), although it is interesting 
that the effect seems to be limited to CC PNe. 


García-Hernández & Górny (2014) also compared their median abundances for each of the several 
dust types with published stellar evolution models by Karakas (2010). Given that DC PNe are 
observed to have solar or supersolar metallicities and elevated N/O values, the authors concluded 
that DC PNe likely had relatively massive progenitor stars (73-5 Mo) that experienced HBB during 
the AGB phase of their evolution. On the other hand, CC and OC PNe with unevolved (amorphous 
and aliphatic) dust tend to have subsolar metallicities and significantly lower N/O values, and are 
likely descended from progenitor stars having masses of less than 3 Mo. However, in a subsequent 
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study of a small subsample of DC PNe for which C/O ratios from faint optical recombination lines! 
could be measured, García-Rojas et al. (2018) showed that updated AGB nucleosythesis predictions 
could explain abundances in some DC objects with progenitor masses of € 1.5 Mo, suggesting that 
DC PNe are not necessarily associated exclusively with massive progenitors. 


Finally, although each element's total relative abundance is conserved throughout the evolution of 
the post-AGB nebula, the form in which the element exists, i.e., ion, atom, or molecule (gas or solid, 
including dust), changes as local conditions such as temperature, density, and the presence of UV 
radiation change. Therefore, when considering elemental abundance measurements that are most 
commonly made by analyzing spectra of the ionized gas in a PN, we need to be aware of how the 
elements are partitioned in the different phases. For example, how much carbon is actually tied up 
in molecules and dust grains instead of the observable ionic and atomic stages? 


Alas, this is a difficult question to answer. According to García-Hernández (2021), a complete 
knowledge of the full inventory of the molecular and dust species in PNe or their immediate AGB 
or post-AGB precursor stages is presently unknown, in spite of having identified a few molecules 
convincingly in the optical, infrared, and radio domains [see e.g., Table 1 in Zhang (2017) Im 
addition, instead of definitive identifications of species from which spectral features arise, researchers 
in this area mostly have candidate carrier species for the several dust features seen in CC PNe (e.g., 
SiC, MgS, HAC, or even mixed aromatic-aliphatic species with O, N, and S impurities) and OC PNe 
(e.g., amorphous/crystalline silicates with different Fe, Si, and Mg content). As an example, Volk, 
Sloan, & Kraemer (2020) have recently reviewed the problems and ambiguities encountered when 
attempting to definitively match carriers to the solid state spectral features located at 30 um and 
21 um in C-rich objects. According to these authors, the 30 wm feature is usually aligned with MgS, 
although they provide reasons to doubt that association, most notably a debate about whether there 
is enough Mg$ to explain the measured intensity. Even less certain is the feature at 21 jum, as its 
detection in PNe is "tentative at best" according to the authors. Thus, for these several reasons, 
the prospect of making a significant step forward in the near future regarding how elements are 
partitioned among the several physical phases is unlikely. 


However, in the absence of positive spectroscopic identifications of molecular and dust species, 
one possible path around this problem is illustrated in the recent analysis by Gómez-Llanos et al. 
(2018) [see also Otsuka et al. (2017)] of the PN IC 418. This team developed a detailed photoion- 
ization model of both the ionized and photodissociated regions of the PN by varying the dust types, 
amounts and distributions relative to the distance to the central star. They were able to match the 
observed IR emission spectrum from 2-200 um by using graphite, amorphous carbon, MgS, and SiC 
as dust ingredients. Their model results showed that refractory elements such as Mg, Si, and S were 
noticeably depleted from the gas phase, while C, N, O, and Ne remained undepleted. Such models 
enable estimates of the amount of depletion of an element, an amount which can then be added 
to the observed gas phase abundance in order to estimate the total elemental abundance in the system. 


16 Delgado-Inglada & Rodríguez (2014) present a comparison of C/O ratios determined from both collisionally excited 
lines and optical recombination lines. They find that the agreement is generally good but caution that in some objects 


the difference is large. 
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5o what have we learned about dust in PNe to date? Substantial progress has been made over 
the past 15 years in identifying specific dust species as well as in understanding the nature of dust 
in PNe, thanks in large part to the availability of quality IR. spectrophotometric data obtained by 
the Spitzer Space Telescope. Observed dust features in the IR coupled with nebular gas abundance 
studies, mostly in the optical and UV, suggest that dust chemistry is largely determined by the value 
of the C/O abundance ratio within the environment in which the dust forms. For example, C-rich 
dust probably forms when C/O>1. Such a rich C environment is likely the result of TDU, in which 
C is transported from the AGB's He-burning shell into the stellar atmosphere and ultimately ejected. 
On the other hand, an O-rich environment can result either when the progenitor star is not massive 
enough for TDU to occur, or in the case of massive progenitors, TDU does occur but HBB converts 
the C to N. Challenges remain before we can claim a thorough understanding of dust in PNe. As 
already pointed out, we cannot yet confidently identify the carriers of many observed dust features. 
In addition, without knowledge of the relative masses of dust and gas in a PN system we cannot 
uncover details of dust production through modeling. Such a condition leaves room for numerous 
PhD theses in the future! 


4. WHAT CAN WE GLEAN FROM PN ABUNDANCE STUDIES? 


It is now time to ask what the outcomes of the numerous PN abundance studies reviewed above 
can teach us about both stellar and galactic chemical evolution. In 84.1 we consider how disk PN 
abundances, when combined with their galactocentric distances, can shed light on disk formation. 
Then in 84.2 we compare elemental abundance information from 83 with stellar evolution model 
predictions in an attempt to learn about stellar nucleosynthesis. 


4.1. Disk Galaxy Abundance Gradients from Planetary Nebula Abundances 


Measurements of chemical abundances in H II regions, PNe, Cepheids, and open clusters in the 
disk of the MW and other galaxies indicate that the metallicity generally decreases with increasing 
radial distance from the center of a galaxy. Recent reviews regarding abundance gradients in galaxies 
include those by García-Rojas (2018) & Maciel & Andrievsky (2019). Chemical evolution models of 
disk galaxies strongly indicate that the slope of an abundance gradient is closely tied to the temporal 
and spatial variation of the infall rate of external gas, the star formation efficiency, and radial gas 
flows [(Fu et al. 2009; Marcon-Uchida, Matteucci, & Costa 2010; Spitoni & Matteucci 2011; Grisoni, 
Spitoni, & Matteucci 2018)]. Even slope differences of the order of a few hundredths of a dex/kpc 
can readily distinguish between models with different input parameter values. 


PNe play an important role in ascertaining the abundance gradient slope in galactic disks, as they 
are visible at large distances and the abundance determination process is relatively straightforward. 
The chemical composition of a PN is determined by the local conditions at the time of progenitor star 
formation as well as the contribution of nuclear products synthesized by the star during its lifetime. 
These two factors must be kept in mind, for example, when using abundance gradients based upon 
PNe to constrain chemical evolution models. 


4.1.1. The Milky Way Gradient 


Examples of PN gradient studies that have been carried out for the MW disk include those by 
Maciel & Kóppen (1994); Maciel & Quireza (1999), Henry et al. (2010), Maciel, Costa, & Cavichia 
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(2015), and Stanghellini & Haywood (2018). What follows is a brief summary of the three most 
recent papers concerning the MW gradient. 


Henry et al. (2010): The PN sample containing both Peimbert Types I and II objects!" employed 
in this study comprised 124 PNe whose galactocentric distances R, were derived from heliocentric 
distances published by Cahn et al. (1992) and Stanghellini et al. (2008, SSV) and ranged from 
0.9-21 kpc, assuming a solar distance of 8.5 kpc. The O/H values came entirely from the authors’ 
own observations and computations. Using the Cahn et al. distances, they computed a gradient for 
12+log(O/H) versus R, by also considering the uncertainties in both distance and abundance and 
found a value of —0.058 + .006 dex/kpc, with a correlation coefficient of r = —0.54. A second regres- 
sion using the SSV distances resulted in a flatter gradient of —0.042 + .004 dex/kpc. In addition, no 
correlation was found between O/H and vertical height from the Galactic plane. 


Maciel, Costa, & Cavichia (2015): These authors combined their own O/H abundance data (Maciel 
& Costa 2013) with like measurements from Stanghellini & Haywood (2010a) to form a PN sample 
of 263 objects. Object distances were determined using the SSV distance scale, assuming a solar 
distance of 8.0 kpc. Objects in the combined sample had a range of « 1 to roughly 15 kpc (estimated 
from their Fig. 2). They obtained a gradient value of —0.025 + .004 dex/kpc and a correlation coeffi- 
cient of r — —0.357. An additional study was made in order to check for a correlation between O/H 
abundance and vertical distance from the Galactic plane. By dividing their sample into two groups 
having either Z < 1000 pc or Z > 1000 pc, the authors computed a linear fit to each group and 
then compared their intercept and slope values. The group closer to the plane had a 12+log(O/H) 
intercept value roughly 0.17 dex (47%) greater and a gradient that was steeper by -0.007 dex/kpc 
than the respective values for the more distant group. 


Stanghellini & Haywood (2018): This recent abundance study of disk PNe in the MW included 
abundances of He, N, O, Ne, and Ar. The sample for PN O/H measurements comprised 248 disk 
objects carefully selected from high quality observational results published by numerous authors. 
Radial distances were determined using the SSV distances and ranged from less than 1 kpc to roughly 
28 kpc. A linear fit that included uncertainty considerations for both distance and O/H yielded a 
gradient of —0.021 + .008 dex/kpc. Two groups were then formed from their complete sample, one 
comprising PNe with progenitor stars younger than 1 Gyr and the other with progenitor stars older 
than 7.5 Gyr, as determined by comparing abundances with stellar model predictions of final surface 
abundances. The younger group was found to have a gradient of -0.027 dex/kpc, while the older 
group's gradient was -0.015 dex/kpc, i.e., shallower than that of the younger group. The authors 
also found that the radial gradient 1) becomes flatter between distances of 10 and 13.5 kpc; and 
2) is steeper near the disk plane compared with its value at greater vertical distances, in qualitative 
agreement with Maciel, Costa, & Cavichia (2015). 


4.1.2. Metallicity Gradients in External Galaxies 


17 Peimbert Types are defined in 84.2 as part of the discussion of Fig. 11. 
18 An attempt to understand the cause of the slope differences among the Henry et al. (2010), Maciel, Costa, & Cavichia 
(2015), and Stanghellini & Haywood (2018) gradients proved unsuccessful. 
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Gradients of O/H, Ne/H, and Ar/H have been measured using PNe in the disks of M31, M33, 
NGC 300, and M81. The gradient in M31 has been recently studied by Kwitter et al. (2012), Balick 
et al. (2013), and Sanders et al. (2012). Magrini et al. (2004); Magrini, Stanghellini, & Villaver 
(2009b); Magrini et al. (2010), Stasiáiska et al. (2005), and Bresolin et al. (2010) have similarly 
studied M33, while Stasiriska et al. (2013) and Stanghellini et al. (2010b) have determined gradients 
for NGC 300 and M81, respectively. 


The results for M31, M33, and NGC 300, along with the MW, have recently been reviewed by 
Pena & Flores-Durán (2019). In their paper, required abundances were collected from the numerous 
extant literature sources, where in most cases the abundances were computed using relevant electron 
temperatures, i.e., the direct method. Separate from the analysis by Pena & Flores-Durán (2019), 
Stanghellini et al. (2010b) used their own optical spectra to determine abundance gradients of O/H, 
Ne/H, and S/H in M81. The results of these two studies are presented in Table 8, where gradient 
values pertaining to combined samples of Type I and non- Type I PNe are provided. Note that each 
gradient has been normalized to the relevant photometric radius, R55?, to allow a direct comparison 
among galaxies possessing different disk radii.Here we can appreciate especially the relative flatness 
of the M31 gradients with respect to those gradients of the remaining four galaxies, while M81 
possesses by far the steepest gradients. Balick et al. (2013) have suggested that M31’s flat gradient is 
the result of a starburst that occurred roughly 2 Gyr ago in a metal-rich ISM in the outskirts of the 
disk following an encounter with M33 approximately 3 Gyr ago. Given the reported uncertainties, 
we also see relatively close agreement among the O/H, Ne/H, and Ar/H gradients of any one galaxy; 
this is likely a reflection of the well known lockstep behavior of alpha elements that is consistently 
seen in PN abundance studies and is discussed in 84.2.1 below. 


4.1.3. Does Oxygen Really Track Metallicity? 


An interesting wrinkle in the abundance gradient picture produced by PNe that has slowly de- 
veloped over the past two decades involves the reliability of using the O/H ratio as a measure of 
metallicity. Péquignot et al. (2000) studied two PNe in the metal-poor Sagittarius dwarf galaxy and 
found that values of Ne/O, S/O, and Ar/O were smaller than the corresponding values found in the 
MW. The authors suggested that He burning followed by TDU in low mass stars (M < 2M.) might 
produce an elevation of both !?C and !9O in the PNe that eventually form. Then Wang & Liu (2008) 
compiled published data for both PNe and H II regions in the MW and several external galaxies and 
reported that in metal-poor galaxies (Z<0.004) both Ne/H and O/H were enhanced in PNe, which 
they ascribed to synthesis of these two elements by LIMS. 


Delgado-Inglada et al. (2015) chose 20 Galactic PNe located in the halo, bulge and disk, all falling 
within +0.25 dex of the solar oxygen abundance and for which both optical and IR data were avail- 
able. They then divided the sample into two groups, i.e., those possessing CC dust and those with OC 
dust, and plotted log values of O/Cl versus O/H. In contrast to the PNe that were associated with 
OC dust (as well as a small sample of H II regions also included in the study) seven of the eight PNe 
with CC dust were offset as a group from the other objects, exhibiting greater O/CI values (see their 
Fig. 3), thereby supporting the conclusion of Péquignot et al. (2000) mentioned above. Since another 


19 Ros is the radial distance at which the blue surface brightness (mg) falls to 25 mag per arcsec?. 
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of their plots indicated strict lockstep behavior between Ar and Cl, the authors suggested that ei- 
ther of these two elements (they preferred Cl) should be used instead of O to gauge metallicity in PNe. 


So what does stellar evolution theory have to say about this matter? As pointed out by García- 
Rojas (2020), AGB stars can synthesize O via He burning and TDU or via HBB and the CNO cycle. 
For example, the stellar evolution models of Pignatari et al. (2016) include convective boundary 
mixing or overshooting to move He shell burning products ^C and '°O into the He-rich intershell 
region and ultimately into the outer atmosphere. They found that the final O mass fraction in the 
outer envelope of stars of mass 2-3 Ms becomes roughly twice what it was in the original material that 
formed the progenitor star. The key to enhanced O here seems to be the inclusion of overshooting in 
the models at the boundary between the He-burning shell and He-rich intershell above it. In addition, 
results of stellar evolution models by Ventura et al. (2017) indicated that when Z«0.004, HBB in 
massive stars destroys oxygen, while TDU in low mass stars causes an increase in oxygen. On the 
other hand, when Z>0.004, O/H in AGB stars shows a tight correlation with metallicity. Thus, at 
low metallicities, O/H may not be a reliable indicator of stellar metallicity. This is a problem which 
is sure to receive much more attention in the future. 


4.2. Measured Elemental Abundances of Planetary Nebula versus Theoretical Predictions 


The principal nucleosynthetic products of LIMS are He, C, N, F, and numerous s-process elements. 
Thus, the abundance measurements of these elements in PNe, as described in 83, serve as valuable 
constraints on LIMS models that make predictions of element levels in the nebular gas. At the same 
time, while O and perhaps Ne may be synthesized in low-metallicity AGB stars during He-burning 
and brought to the surface by TDU (see 84.1.3), generally speaking, the abundances of alpha elements 
O, Ne, S, Cl, and Ar observed in PNe reflect the amounts that were originally present in the star- 
forming gas of the progenitor star. This abundance dichotomy allows PNe to play a second role, that 
of establishing the metallicity of the ISM at the time and location of the formation of PN progenitor 
stars. We shall address these two features in reverse order. 


4.2.1. Planetary Nebula Abundances As Tests of Alpha Element Lockstep Behavior 


In 83.2.4, we briefly described the difficulty of using Fe, the usual indicator of metallicity in stars, 
to establish metallicity in nebulae such as PNe. On the other hand, gauging metallicity in PNe is 
customarily done by measuring abundances of O (usually) or other alpha elements such as Ne, S, Cl, 
or Ar. Since the bulk of these elements is mostly synthesized via He burning in massive stars and 
expelled during Type II supernova events, their interstellar abundances are expected to increase in 
lockstep. 


A direct method for exploring this lockstep behavior using PNe is to form element-versus-element 
plots of alpha element abundances determined for a large survey and then to look for strong linear 
correlations between element pairs. Observed abundance data from KH20 are shown with open 
circles in Figure 10 for Ne/H, S/H, Cl/H, and Ar/H versus O/H, where we have used the standard 
logarithmic form of 12+log(X/H) and X and H are number abundances. The error bars in the 
lower right corner of each panel indicate the average uncertainty expressed logarithmically. We also 
show least-squares fits in each panel, where the solid line represents the fit to the PN data. In the 
cases of Ne, S, and Ar, the dashed line shows the fit to a collection of published abundances for 
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H II regions and blue compact galaxies (H2BCG) compiled by Milingo et al. (2010, Table 6). The 
analogous fit for Cl was inferred from a study of Galactic disk H II regions by Esteban, García-Rojas, 
& Pérez-Mesa (2015). Table 9 displays the regression parameters and uncertainties, i.e., y-intercept, 
slope, correlation coefficient, and number of sample objects considered for each line. The lockstep 
behavior of alpha elements O, Ne, S, Cl, and Ar in the ISM, as probed by the H2BCG, is clearly 
demonstrated by the near-unity slope values and correlation coefficients. The PN results are a bit 
messier (no pun intended!), although with the exception of Ne/H, the slopes are strongly suggestive 
of lockstep behavior. The larger amount of scatter demonstrated by the PNe for all the alpha 
elements may be due in part to the generally higher levels of ionization within the gas, which can 
complicate corrections for unobserved ions when determining total elemental abundances. A different 
explanation may be that one or more of the original alpha element abundances at the time of star 
formation were subsequently altered by nucleosynthetic processes within the star during its lifetime 
as just discussed in 84.1.3. 


Another way of using PNe to confirm the direct abundance relation among alpha elements is to 
compare the slopes of abundance gradients of these elements in galactic disks, where PNe that are 
positioned over a wide range of galactocentric distances are used as abundance probes. Using this 
method, good agreement among gradient values for O/H, Ne/H, and Ar/H has been illustrated 
by Bresolin et al. (2010) for PNe in the disk of M33 and by Stanghellini & Haywood (2018) for 
Galactic disk PNe. Similarities of alpha element gradient slopes in galactic disks are also predicted 
by chemical evolution models such as those published by Chiappini, Matteucci, & Gratton (1997). 


To summarize, with the possible exception of low metallicity regimes (Miller Bertolami 2021), alpha 
element abundances appear to evolve in lockstep in the ISM, a pattern which is then observed in 
PNe that form from that gas. 


4.2.2. Gauging LIMS Nucleosynthesis With Planetary Nebulae 


We now consider elements which are synthesized by LIMS during their lifetime and subsequently 
ejected into the ISM. Stellar models predict both the total contribution of such elements to the local 
environment over the star's lifetime, i.e., the yield, as well as the abundances of these elements in the 
stellar envelope (surface abundances) just prior to ejection and subsequent PN formation. A com- 
parison of PN abundance measurements with predicted surface abundances serves as a constraint on 
the stellar models. A successful model then not only implies the progenitor star's mass and element 
yields, but it also helps to confirm the role of the physical processes that were part of the model 
calculation. Therefore, the back and forth comparison of models and PN abundance observations 
teaches us much about the contributions that LIMS make to the chemical evolution of galaxies. (See 
Appendix B for a list of publications of model-predicted stellar yields and surface abundances of 
LIMS.) 


In this subsection we present and discuss plots of observed values of He/H, C/O, and N/O in MW 
disk PNe. Following that, we compare the observations with updated predictions of the final AGB 
surface abundances made by the same four codes featured and compared in detail in Henry et al. 
(2018), i.e., MONASH [Karakas & Lugaro (2016); Karakas et al. (2018)], LPCODE (Miller Bertolami 
2016), ATON [Ventura et al; (2020) and Marini et al. (2021)], and FRUITY [Cristallo et al. (2011); 
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Piersanti, Cristallo, & Straniero (2013); Cristallo et al. (2015, 2016)]. 


Observed abundance ratios of He/H, C/O, and N/O for the MW disk PNe contained in KH20 are 
displayed with circles in each of the panels of Figs. 11 and 12. In the upper panels of Fig. 11 filled 
circles represent objects whose C abundances were strictly based upon UV data (C III] 1909À and 
CIV 1549À), while open circles indicate those objects where only optical data (C II 4267À) were 
used instead because UV data were unavailable. The sample error bars in the left-most panels show 
median values of the observational uncertainties. In the bottom panels, the upper right quadrant 
marked off by a pair of dashed lines shows the region occupied by Type I PNe. This classification 
was originally introduced by Peimbert (1978), who noted that objects in this group had noticeably 
enhanced abundances of He and N, favored bipolar morphology, and resided in the MW disk??. 
After nearly two decades, and in light of the availability of many additional abundance results in 
the literature, Torres-Peimbert & Peimbert (1997) recommended that MW Type I PNe be defined 
abundance-wise by He/H 70.125 and N/O>0.5. 


'The sample of abundance ratios in the bottom panels of Fig. 11 clearly suggests a direct relation 
between log(N/O) and He/H. A regression analysis of the data produces a correlation coefficient 
of +0.55 between these two variables with a probability of «0.001 of computing such a value from 
an uncorrelated sample. One can see a qualitatively similar correlation in both the Aller & Czyzak 
(1983) and Kingsburgh & Barlow (1994) data (not shown here). In view of the discussion presented 
in 82.1, a log(N/O)-He/H correlation can be reasonably interpreted as the result of the major 
conversion of ^C to ‘N via the HBB process at the base of the H-rich envelope as the He/H ratio 
rises steadily, where both increases are directly related to progenitor mass. This last point is clearly 
illustrated in Fig. 13, where we have plotted values of predicted surface abundances versus mass from 
Karakas & Lugaro (2016); Karakas et al. (2018). In the bottom panel note the abrupt increase in 
log(N/O) beginning at roughly 4 M5, where HBB becomes effective, especially under low metallicity 
conditions. The top panel shows the predicted steady rise of He/H with progenitor mass. These 
trends are also consistent with the notion that Type I PNe result from more massive progenitor stars 
(Torres-Peimbert & Peimbert 1997). 


However, the notion that Type I PNe are exclusively associated with massive progenitors is chal- 
lenged by other empirical results. For example, Stasiriska et al. (2013) observed several Type I 
objects in a sample of 26 PNe in NGC 300, while at the same time estimating the masses of their 
progenitor stars to be within the range of 2-2.5 Mo. In addition, assuming that selection effects are 
not a significant factor, Type I PNe comprise 30% of the entire KH20 sample presented in Fig. 11, 
a proportion significantly larger than predicted by a standard initial mass function if these objects 
are solely the product of massive progenitors. Indeed, Stasiriska et al. (2013), Karakas & Lattanzio 
(2014), and others have suggested that many Type I PNe originate from progenitor stars less massive 
than 3 Mo in which rotationally induced mixing on the main sequence may have resulted in elevated 
N/O following the first dredge-up. Thus, stars with M<3 M; may play a bigger role in influencing 


20 Peimbert (1978) also defined types II, III, and IV for PNe characterized by intermediate population, high velocity, and 
halo location, respectively. However, most conversations today regarding Peimbert types are concerned with whether 
an object is a Type I or non- Type I object, especially when abundances are of prime interest. 
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PN nitrogen abundances than previously thought. 


Regarding the abundances in the upper panels of Fig. 11, note that the absence of a correlation 
between C/O and He/H is confirmed by the presence of a near-zero correlation coefficient for these 
points. Again, we examined the data in the Aller & Czyzak (1983) and Kingsburgh & Barlow (1994) 
papers and found no suggestion of a correlation between C/O and He/H in either reference. 


Continuing on with the same data set, we have also plotted log(C/O) vs. log(N/O) in each of four 
panels of Fig. 12. Symbols in this figure are defined as they are in Fig. 11. Each panel features 
a different unique set of model predictions for our upcoming comparison with the data. Note the 
absence of a correlation between the two abundance ratios, which we can confirm by the presence 
of a near-zero correlation coefficient for these points. Again, we checked the data in the Aller & 
Czyzak (1983) and Kingsburgh & Barlow (1994) papers and found no suggestion of a correlation 
between C/O and N/O in either reference. The absence of a correlation in this case may be related 
to the initial increase followed later by a decrease of C/O as progenitor mass increases, all while N/O 
increases monotonically with mass. 


We now consider the models and what they can tell us about the observed ratios. The surface 
abundance predictions of published models of AGB stars after the final pulse from the MONASH, 
LPCODE, and ATON groups referenced above are shown with X symbols in Fig. 11. Models of the 
same three groups plus those of the FRUITY group are shown in Fig. 12. Each panel in the two 
figures contains three model sets, where all models within a set have the same metallicity and symbol 
color, as defined in the legend. Each model within a metallicity set has a unique progenitor mass 
(not shown). For each metallicity set, Table 10 provides details regarding the source of the models, 
metallicity, progenitor mass range in solar masses, and the estimated onset mass of HBB. This last 
quantity is the progenitor mass at which log(C/O) reaches its highest value before beginning a steady 
descent in a plot (not shown) of log(C/O) versus mass. 


Fig. 11 clearly demonstrates that as set metallicity increases (red to green to blue), set members in 
each panel collectively tend to shift toward lower log(C/O) or log(N/O) values and to higher values 
of He/H. As discussed by Karakas & Lattanzio (2014, 84.1) and references therein, the efficiency of 
'TDU is inversely related to metallicity but directly related to progenitor mass. In the upper panels 
of Fig. 11, the downward shift of the metallicity sets is consistent with progressively less C dredge-up 
(via TDU) as metallicity increases. In addition, as the O abundance rises with metallicity, the direct 
effect of C dredge-up on the value of C/O is likely dampened to a greater extent. The drop in 
dredged-up C with increasing metallicity also limits the amount of that element that can be used to 
synthesize N through HBB and can explain the downward shift in log(N/O) in the lower panels. 


The rightward shift with increasing metallicity tells an equally interesting story. Look at the top 
panel of the MONASH models in Fig. 11 and consider the low metallicity (red) set for a minute. 
With rising He/H, and thus the stellar mass (top panel, Fig. 13 again), we see an initial increase 
in log(C/O) followed by an abrupt drop of 0.78 dex at He/H=0.113, a drop presumably due to the 
onset of HBB. This behavior is repeated in the two higher-metallicity model sets but is more difficult 
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to see, due to the density of the observed data points. In the lower panel of Fig. 11 the abrupt rise in 
log(N/O) in the MONASH models coincides with the sudden drop in log(C/O), likely at the point 
where HBB begins to ramp up. Note again that the models representing the more massive progenitor 
stars appear in the upper right region defined by the Type I PN abundance parameters. A different 
but consistent version of the above scene is played out in Fig. 12. Here we can see how progenitor 
mass splits each metallicity set horizontally into two subgroups. The ATON models also exhibit 
the trends especially well because their output represents three distinct values of He/H instead of a 
continuum. 


To summarize: for constant mass, C/O and N/O appear to correlate indirectly with metallicity but 
directly with mass when metallicity is held constant. 'Thus in Fig. 11, according to the models the 
greater the vertical and horizontal distances of a PN from the origin, the lower its metallicity and 
greater its progenitor mass is likely to be, respectively. According to our current understanding of 
LIMS evolution, as reviewed above, these effects are most likely due to the sensitivity of both TDU 
and HBB to each of these two stellar parameters. 


Fiinally, a general point about the models. The areas occupied by the observations appear to 
be spanned moderately well by the models in both Figs. 11 and 12. No models were available for 
He/H«0.10 and the region where He/H 0.15 is only sparsely populated with models. Models also 
appear to be lacking for log(N/O)<-0.75. But the main core of observed points falling between the 
two He/H extremes appears well covered. This suggests that most if not all of the observed points 
may be successfully modeled simply by making small adjustments in input parameter values and/or 
internal processes. Unfortunately, the problem then shifts to assessing the uniqueness of each model, 
but so far the collective results are encouraging. At least that's our optimistic glass-is-half-full 
interpretation. 


Besides the production of He, C, and N, another significant contribution made by LIMS to the 
chemical composition of the ISM includes the group of s-process elements which were discussed earlier 
in 83.2.1. According to Karakas & Lattanzio (2014, 85.2.5), roughly half of the elements beyond the 
Fe peak are produced by s-process neutron capture, and AGB stars are responsible for most of that 
production. S-process elements are likely produced in the He-burning shell prior to TDU, where 
the source of free neutrons is either the !?C(o,n)!6O or ?Ne(a,n)? Mg reaction, depending upon the 
local temperature. 


One of several good examples of the predicted impact of AGB stars on the chemical evolution of 
s-process elements in the Galaxy is provided by the recent paper by Prantzos et al. (2018). In their 
study, the authors provide results of a detailed one-zone chemical evolution model of the proto-solar 
nebula that predicts the abundances of all 285 stable isotopic species from hydrogen to uranium. 
Their successful model employed the stellar lifetimes and metallicity-dependent yields of Cristallo 
et al. (2015, FRUITY code) for the mass range 1-7 Mo and of Limongi & Chieffi (2018) for stars 
exceeding 7 Mg. The Cristallo yields for LIMS were computed for stars within the metallicity range 
for [Fe/H] of -2.4 to 0.2., while the metallicity range over which massive star yields were computed 
was -3 to 0 for [Fe/H]. In addition, the massive star yields included the effects of stellar rotation on 
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internal mixing. The final model generally reproduced the abundances of the major isotopes between 
12C and Fe to better than 15%. 


The top panel of our Fig. 14, taken from Prantzos et al. (2018), shows the mass fraction of the 
proto-solar abundance (Lodders, Palme, & Gail 2009, for t — -4.5 Gyr) of s-process isotopes that 
are produced by LIMS only (blue open circles), or LIMS and either non-rotating (green crosses) or 
rotating (red filled circles) massive stars as a function of mass number A. The blue vertical error bars 
show the lo observational uncertainties, while the dotted lines show the 10% and 50% deviations 
from the solid line representing the proto-solar abundances. For orientation purposes, the atomic 
mass numbers of the observed s-process elements discussed in 83.2 and listed in our Table 3 range 
from 70-134 (Ge to Xe). 


Since the vertical gap between the blue symbols and the proto-solar line above them represents the 
contribution to the total abundance that is made by LIMS, it is clear that this contribution rapidly 
increases between A=70 and A=95; after that point, the s-process element production appears to be 
completely dominated by LIMS, according to the model predictions. We also note that below about 
A=85, rotating massive stars have a large impact on the abundance of these elements. The middle 
panel illustrates the same results by explicitly showing the fraction of the total s-process production 
of an isotope represented by LIMS. For A795, the fraction is roughly 9596. The bottom panel shows 
the variation of the production factor?! of a given isotope with respect to that value for ?Sm. 


Finally, considering the broader picture that includes species besides s-process elements, what does 
the Prantzos et al. (2018) model suggest about the impact of LIMS on the chemical enrichment of 
the ISM? The authors note that about 1/3 of the total C in their model is contributed by LIMS, 
while the remaining 2/3 originates in massive stars. Similarly, LIMS are responsible for about 4596 
of N production, while the balance comes from massive stars. (The authors also found that the total 
production of N that is necessary to match observations required the inclusion of rotational effects 
in massive stars.) Therefore, the importance of LIMS for creating two basic elements necessary for 
nearly all biological life is obvious. Regarding PF production, Prantzos et al. (2018) point out that 
both ^N and !?C are consumed within the He-burning shell to produce this isotope. In the final 
accounting, rotating massive stars are responsible for about 2/3 of total stellar production of !?F, 
while LIMS make up the remaining 1/3. And if we now add in the contribution of LIMS to s-process 
element abundances in the ISM discussed above, we conclude that stars between 1 and 8 Mo are vital 
sources of important elements. To borrow a quote from Karakas & Lattanzio (2014), “Stellar yields 
are a key ingredient in chemical evolution models. Low- and intermediate-mass stars are an integral 
part of galaxies and help shape their evolution, gas and dust content, as well as their integrated light. 
Even stars as low as 0.9 Mo can, at low metallicity, contribute to the chemical evolution of elements. 
The days of only considering supernovae are over.” 


5. PLANETARY NEBULA DISTANCE DETERMINATIONS 


To convert many observed properties of PNe into astrophysically useful ones requires knowledge of 
the nebular distances. Gaia (GAIA Collaboration et al. 2018) trigonometric distances are becoming 


?! The production factor is the ratio of the yield to the solar abundance. 
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the standard for Milky Way PNe (except in cases where the central star is too faint), but for com- 
pleteness and context we review below the suite of PN distance methods that have been and continue 
to be in use. Straightforward in principle, the determination of accurate distances for PNe has been 
far from simple. Neither the central stars of PNe nor the nebulae themselves present as easily 
characterized standard candles. Carefully chosen subsamples may exhibit some mean parameter 
values, but with significant dispersion, resulting in unsatisfactorily constrained distances. 


Distances for individual Galactic PNe can be derived from central star properties (e.g., trigono- 
metric parallax, companion spectroscopic parallax, gravity distance) and/or nebular properties (e.g., 
reddening, kinematics, expansion parallax). Membership in a star cluster affords additional con- 
straints on stellar and nebular parameters, including distance. For PNe in nearby galaxies, distances 
are presumed to be known, at least within the uncertainty of a galaxy's distance and physical ex- 
tent along the line of sight. In 85.1 we discuss numerous methods commonly used to determine PN 
distances. In 85.2 we will discuss the related matter of the Planetary Nebula Luminosity Function 
(PNLF), whose bright end has been used as an extragalactic standard candle to establish distances 
as part of the cosmological distance ladder. 


5.1. Methods 


Distance methods can be classified as primary or secondary. Primary distance methods (e.g., 
trigonometric parallax) do not rely on generic theoretical models of stars or nebulae and are focused 
on individual PNe; secondary methods (e.g., the “Shklovsky method"), in contrast, are statistical in 
nature and are based on analysis of the correlated behavior of some property/ies of the star or nebula, 
and must be calibrated against PNe whose distances have been obtained independently using primary 
methods. Below we summarize the most important of these methods; a more detailed discussion of 
measurement error and systematic bias in various PN distance techniques can be found in Smith 
(2015). Frew et al. (2016, hereafter F16) also give a detailed discussion and assessment of distance 
measurement methods; their Tables 1-10 list results (theirs and others’) from the techniques they 
describe. 


5.1.1. Primary methods 


Trigonometric Parallax of the PN Central Star 

Trigonometric stellar parallax is the most direct distance method, yet prior to Gaza, it was only 
possible for a small number of nearby PNe: of the 12 objects from the literature listed in Table 1 
of F16, 11 are within 750 pc.” The fractional error for parallax measurements is typically 15-40% 
or more. Harris et al. (2007) carried out a ground-based program at the US Naval Observatory 
that yielded parallaxes for 16 objects (their Table 3). Using HST, Benedict et al. (2009) measured 
parallaxes for four nearby PNe in Harris et al.’s sample, finding good agreement. 


Fortunately, the advent of Gaia has meant that the number of PN central stars with measurable 
parallaxes will continue to grow. Stanghellini et al. (2020) studied over 400 PN central stars from 
DR2 (their 2020 Table 1). Schónberner & Steffen (2019) used DR2 parallaxes for 11 identified central 


22 The remaining PN, K3-35 (d = 3900 pc), is thus far the only PN to have its parallax derived from observations of 
maser emission (Tafoya et al. 2011). 
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stars (their Table 1) to compare with the expansion parallaxes they derived, finding good agreement 
(see section on expansion parallaxes, below). Kimeswenger & Barría (2018) were able to manually 
identify almost 400 PN central stars in DR2. González-Santamaría et al. (2019, hereafter GS19) 
derived distances for 211 PN central stars with good quality Gaia DR2 data. González-Santamaría et 
al. (2021, hereafter GS21) have extended this work with EDR3 data, deriving good quality distances 
to 405 PN central stars; comparing distances in common between the two datasets, good agreement 
is found. Anticipating the burgeoning of both the Gaia archives and PN catalogs, Chornay & Walton 
(2019, 2020, 2021) devised an automated technique for efficient identification and cross-matching of 
PN central stars and rejection of potential contaminants. Applied to the HASH PN Catalog (Parker 
et al. 2016), the method yielded more than 1000 likely Gaia PN central stars, more than twice the 
number from previous manual matching attempts. 


Binary Membership 

Binarity provides additional opportunities to ferret out distances of PNe. If a PN central star has a 
resolved binary companion of a standard spectral type, then the established method of spectroscopic 
parallax can be applied to it, and a distance determined for the system. PN central stars with re- 
solved companions represent a minority of PNe with known binary central stars (see 86 on binaries). 
The best known example is NGC 246. From ground-based photometry Bond & Ciardullo (1999) 
concluded that the central star companion, at a projected separation of 3.8 arcsec, is a G8-K0 main 
sequence star, implying a distance of 495 pc, in good agreement with the Gaia EDR3 value of 530 
pc (GAIA Collaboration et al. 2021). Ciardullo et al. (1999) carried out an HST snapshot survey to 
search for additional resolved PN pairs and found 10 “likely” and six “possible” binary associations. 
As one might expect, calculated distances for these objects range far beyond those for ground-based 
parallax measurements, extending to several kpc, with fractional distance errors averaging around 
2596 (F16, Table 2). A small number of PNe have eclipsing binary nuclei (see 86); astrophysical 
analyses of these systems offer additional distance estimates (F16, Table 3). González-Santamaría 
et al. (2020) found eight wide binaries (including one triple system) among their good-quality Gaia 
DR2 distance sample described in GS19. Most recently, GS21 detected three further possible wide 
binary systems in their EDR3 sample, including the Ring Nebula (NGC 6720). 


Cluster Membership 

A PN located in a star cluster is a rare celestial gift. Reliable methods exist for determining 
cluster distances: multi-color photometry of cluster members allows fitting isochrones of appropriate 
metallicity and age, and distances based on cluster membership typically have errors around 10% or 
better (F16, Table 4). In addition, cluster membership allows the determination of the PN progenitor 
mass. What is not always reliable is the conclusion that a PN is physically associated with a cluster. 
There are examples of such associations being claimed, then deprecated, then possibly re-claimed, on 
either re-analysis or additional evidence (e.g., NGC 2818 in the open cluster of the same designation, 
and Abell 8 in the open cluster Bica 6). The current status of PNe in clusters is summarized by 
Davis et al. (2019). 


Fragkou et al. (2019a) report the association of a PN, BMP J1613-5406, with the young (~90 Myr) 
open cluster NGC 6067. Two Cepheids in the cluster indicate a distance of 1.75+0.10 kpc (Majaess 
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et al. 2013); Gaia DR2 parallaxes for 43 cluster stars yield 1.94+0.07 kpc. Using the PN Ha surface 
brightness-radius relation (see see 85.1.2 below). Fragkou et al. (2019a) derive a distance of 1.71 
+0.29/-0.24 kpc, in good agreement with the other methods. A PN in such a young cluster implies a 
massive progenitor, around 5 Mọ; such PNe are rare and of exceptional interest. Another intriguing 
example is the distant, intermediate-age (70.7 Gyr) open cluster Andrews-Lindsay 1 containing 
the PN PHR 1315-6555, most recently studied by Fragkou et al. (2019b), who derive a distance of 
12.0+0.5 kpc, a 4% error; they are even able to detect the central star. 


As for globular clusters, there are currently four in the Milky Way containing an identified PN: 
Ps1 (whose central star is named K648, often used for the PN as well) in M15 (Pease 1928), GJJC 
1 in M22 (Gillette et al. 1989), JaFu 1 in Palomar 6, and JaFu 2 in NGC 6441 (Jacoby et al. 1997). 
Recent proper motion analysis of these PNe support their cluster membership, except for JaFu 1, 
which still requires more confirmatory data (Bond et al. 2020). Surveying globular clusters in the 
direction of the Galactic bulge, Minniti et al. (2019) identify four "excellent candidates" for PNe, 
one each in four clusters. If confirmed, these would double the known number of globular cluster PNe. 


Bond (2015) carried out a largely unsuccessful HST snapshot survey to detect PNe in the glob- 
ular clusters of Local Group galaxies; he did, however, find a PN with a relatively massive central 
star in an open cluster in M31 (Davis et al. 2019). Larsen (2008) identified a PN in a globular 
cluster in the Fornax dwarf spheroidal galaxy. Searches have been undertaken for PNe in star 
clusters belonging to more distant galaxies, with limited success. Out of 80 clusters observed in 
four galaxies, Larsen & Richtler (2006) detected two PN candidates in young clusters: one in M83 
and one in NGC 3621. Most recently Sun at al. (2019) identified one PN in a globular cluster in M87. 


Self-consistent Nebula + Star Photoionization/Spatiokinematic Modeling 

An approach that produces a self-consistent model description of observed and predicted nebular 
and stellar properties can be a useful exercise, yielding distances with errors around 2096. Using 
photoionization codes like Cloudy (Ferland et al. 2017) and MOCASSIN (Ercolano et al. 2003) 
one tries to construct a self-consistent description of both the star and the nebula. The method is 
data-hungry, requiring spatially resolved narrow-band imaging, spectrophotometry and kinematics, 
and is more likely to be feasible for nearby angularly large and bright PNe. The method has been 
successfully applied to about a dozen PNe; see F16, Table 7 and references therein. 


Gravity Distances 

In principle, if both a star's absolute flux and observed flux are known its distance can be calculated 
via the inverse square law. The gravity distance method uses NLTE model atmospheres to derive 
a PN central star's effective temperature and surface gravity, yielding the absolute flux, which is 
then coupled with the observed magnitude and interstellar reddening toward the star to determine 
the reddening-corrected observed flux, and thus the distance. The method is described in detail by 
Méndez et al. (1988). Despite its apparent simplicity, the application of this method is challenging 
and dependent on model parameters of the stellar atmosphere (see, e.g., Rauch et al. 2007). F16 
list updated gravity distances from the literature for ~60 PNe in their Table 5; the average distance 
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error is about 30%. 


Expansion Parallaxes 

This is another method that is simple in principle, but susceptible to complications. The spectrum 
of a spherical nebula expanding homogeneously will exhibit a Doppler shift due to the expansion. 
Over time, as the nebular radius increases, its angular growth can be measured. Assuming the 
angular rate of change is fueled by the spectroscopically determined radial velocity, the distance can 
be calculated. Unfortunately, PNe are rarely spherical, and more importantly, the observed angular 
motion of a specific nebular feature does not necessarily correspond to the expansion of the nebula 
as a whole, as described by Mellema (2003). As might be expected, PNe with measurable expansion 
parallaxes are bright and within a few kpc. 


Schónberner et al. (2018) used spherically-symmetric 1D radiation hydrodynamical models (see 
87.3) to predict the expansion behavior of various nebular components (shell, rim, shock front), pro- 
viding correction factors to be applied in deriving the distance from the measured radial velocities; 
the factors are greater than 1, implying larger distances. Their Table 7 lists expansion distances for 
15 PNe, with average quoted errors around 2096. Schónberner & Steffen (2019) compared results 
for 11 of these objects that have Gaia DR2 trigonometric parallaxes: the agreement is very good as 
shown in their Table 1 and Figure 2: Schónberner (2021) reports that Gaia EDR3 (GAIA Collabo- 
ration et al. 2021) agreement is even better, with most errors under 10% and less than those for the 
expansion parallaxes, though some outliers remain. 


Distance Mapping 

The distance mapping technique (Akras & Steffen 2012) is a hybrid process combining the expan- 
sion parallax method along with model-derived 3D morpho-kinematics. Gómez-Gordillo et al. (2020) 
used an updated version of this method to derive distances and for several bright PNe, including 
NGC 6302, NGC 6543, and BD+30 3639. Their values agree well with previous determinations from 
the literature, though for Gaza comparisons, they cite the difficulty of accounting for the source- 
dependent zero-point offsets in the Gaia DR2 parallaxes. 


Extinction Distances and Kinematic Distances 

These methods take advantage of the Galactic environment and line-of-sight to the PN. Extinction 
distances are estimated by comparing the extinction of the nebula and/or central star with angularly 
nearby stars whose independently determined distances bracket that of the PN. This is practical 
only for PNe within a few degrees of the Galactic plane; at higher latitudes the dust distribution 
thins and becomes less predictable. Smith (2015) concludes that individual extinction distances are 
generally unreliable. F16 lists 31 extinction distances for PNe in their Table 9, with objects selected 
for individual reported errors <20%. 


To understand the concept of a kinematic distance, consider a PN moving with the general Galactic 
circular rotational velocity at its location, and having only a small peculiar velocity. These constraints 
automatically restrict the method to young, disk PNe, most likely of Type I (see §4.2), originating 
from progenitors of higher mass than bulge or halo PNe. For any given Galactic rotation curve and 
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PN sky position, there will be a relationship between the PN’s distance from the Galactic center (and 
from us), and its observed radial velocity. Distance uncertainties from this method derive primarily 
from uncertainties in the assumed rotation curve and from the size of the velocity dispersion for 
sources of the appropriate age (Phillips 2001). F16 list 20 objects with kinematic distances in their 
Table 8; the quoted errors are around 30%. 


5.1.2. Statistical Methods 


Prior to Gaia, the vast majority of MW PNe were not amenable to direct distance methods 
and astronomers developed statistical techniques that related measurable PN parameters to infer a 
distance. The continuing value of these methods lies in their potential to reveal characteristics of 
different populations of PNe. Statistical techniques have varying degrees of robustness; the main 
difficulty lies in the inherent variation among PNe as well as in evolutionary effects. Methods that 
attempt to account for these tend to be more reliable. Statistical methods must be calibrated using 
PNe that do have direct distances available. Over the last several decades, differences in calibration 
efforts led to scales dubbed “long” (e.g., Zhang 1995) and “short” (e.g., Cahn et al. (1992, hereafter 
CKS) with distance values disagreeing by up to a factor of three. 


The first widely applied statistical method was developed by Shklovsky (1956). Assuming that all 
PNe contain the same mass of ionized material (he assumed 0.2 Mọ), there is a relationship between 
a PN’s (extinction-corrected) H5 surface brightness (SB) and its linear size, so as the nebula expands, 
its electron density (which regulates the surface brightness) decreases. Amassing SB and angular 
size measurements for a set of target PNe, then calibrating with PNe having independent distance 
determinations allows the distances for the target PNe to be calculated. The assumed nebular mass 
actually has a relatively small effect on the resulting distance, which goes only as M?*. Variations on 
the Shklovsky method have included substituting Ha (Abell 1966) or 5 GHz radio fluxes (Daub 1982) 
for HG, as well as eliminating the assumption of a fixed ionized mass, and allowing for consideration 
of both optically thin and optically thick PNe (e.g., CKS; Stanghellini et al. 2008). 


Frew (2008) first presented a statistical scale based on the SB in Ho and the nebular radius. The 
astronomical observables required are the angular size of the nebula and the extinction-corrected 
total Ha flux. The method, refined by F16, includes calibration PNe over a wide range of surface 
brightness along with distances determined by an array of methods; the lack of such broad inclusivity 
had hampered previous statistical methods, leading to the conflicting distance scales, and increasing 
the chances of Malmquist bias affecting the results. They give the best-fit relation between Ha SB 
and nebular radius for their full sample as: 


log SHa = —3.63(+£0.06)log r — 5.34 + 0.05 (1) 


where Sg, is in units of erg cm? s^! sr^! and r is in pc. F16 also give SB-r relations for specific 


nebular subsamples, e.g., known optically thick or optically thin, bulge, and extragalactic PNe. In 
his analysis of multiple statistical methods and the resulting distance scales, Smith (2015) concluded 


23 Analysis of a flux-limited sample may be subject to Malmquist bias (Malmquist 1924), where only brighter objects are 


detected at large distances, leading to erroneously higher average luminosity with increasing distance. A related effect is 
the Lutz-Kelker bias (Lutz & Kelker 1973), where increased distance uncertainties result in statistically underestimated 
distances, thereby affecting derived distance-scale slopes (see, e.g., Bensby & Lundstróm (2001)). 


38 


that the F16 scale is the best currently available, with the only caveat concerning very evolved PNe 
whose distances seem to be systematically underestimated. Kimeswenger & Barría (2018) found that 
the Gaia distance scale agrees with that of F16 to within ~ 8%. 


5.2. The Planetary Nebula Luminosity Function 
5.2.1. Basic Concept 


The varied morphologies and wide range in observed brightness of Milky Way PNe, coupled with 
their difficult-to-determine distances make them unlikely candidates for consideration as a standard 
candle. However, just as Henrietta Leavitt exploited the fact that all Cepheids in the SMC are at a 
(very nearly) identical distance so that apparent differences in brightness reflect absolute differences, 
the same can be done with the ensemble of PNe in an external galaxy. The Planetary Nebula Lumi- 
nosity Function (PNLF) has become a valuable extragalactic distance indicator with the significant 
advantage that PNe originate from low- and intermediate-mass stars found in all types of galaxies, 
so that the PNLF can be applied in active star-forming systems as well as in more sedate, quiescent 
galaxies. 


The brightest PNe in each galaxy exhibit similar absolute fluxes in their strongest emission line, 
[O III] A5007, suggestive of a specific maximum PN luminosity (the “cutoff” ). Models indicate that 
up to ~10% of a PN central star's total luminosity (L,> 10?L5) can be concentrated in this nebular 
line (Dopita et al. 1992; Ciardullo 2010; Schónberner et al. 2010), meaning that [O III] surveys can 
reveal PNe out to large distances; as mentioned in 81, PNe have been detected out to ~100 Mpc. 
These two characteristics (known luminosity and high luminosity) are just what is required in a 
standard candle. Further refinement is possible by using not just the few brightest PNe but all of the 
PNe observed in a given galaxy (Jacoby et al. 1988) and fitting the observed magnitude distribution 
to a theoretical model function to yield the distance. 


5.2.2. The Planetary Nebula Luminosity Function 
The absolute [O III] magnitude of a PN, Msoo7, is defined by Jacoby (1989) as: 


Msooz = —2.5 log(F5007) — 13.74 (2) 


where F5oo; is the absolute flux of the nebula in the 5007 line in ergs cm~*s~! when viewed from a 


distance of 10 pc. Substituting the observed flux, f55o7, yields the apparent [O III] magnitude, m5007. 
The relationship between Moo; and moo; is the same as for any observed vs. absolute magnitude: 
m — M = 51log(d,./10), after correction for foreground reddening. As a point of reference, the 
brightest PNe in M31 have msoo; ^20 (Merrett et al. 2006), while those in the Virgo Cluster have 
msooz ~26 (Jacoby et al. 1990). 


An observed PNLF is constructed as a histogram of the number of PNe per moo; interval (typi- 
cally ~0.2 mag) which can then be fit to a numerical function incorporating the characteristics and 
evolution of an ensemble of PNe and their progenitor stars. Based on their observations of PNe in 
the Small Magellanic Cloud, Henize & Westerlund (1963) derived an approximate PNLF assuming 
a narrow nebular mass range, uniform expansion, and a constant PN formation rate. In a study 
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of PNe in M31 and its companions, Ciardullo et al. (1989b) found that they needed to modify the 
Henize & Westerlund (1963) formulation to reproduce the sharp cutoff observed at the bright end of 
the distribution, leading to the following expression for the empirical PNLF: 


N(M) x gg 0 pim (3) 


where N(M) is the expected number of PNe in a bin at magnitude M, and Mx, the cutoff magnitude, 
is the brightest magnitude observed. The normalization constant will depend on the number of 
PNe in the sample. The value of Mx = —4.54 +0.05 (Ciardullo 2013) has been shown to be invariant 
over stellar populations and galaxy types, with only a slight dependence on metallicity (e.g., Dopita 
et al. 1992; Schonberner et al. 2010). Figure 15 shows the PNLF for M31 out to about 30 kpc from 
the center. Note the sharp cutoff at the bright end, at moo; ~ 20.17. Figure 16 shows the PNLF 
for five galaxies in the Leo I group (Ciardullo et al. 1989a) whose distance is ~10 Mpc. Despite 
the spread in Hubble type from E0 to SBb — and therefore the stellar populations sampled — the 
bright-end behavior of the PNLF is remarkably consistent. This observed consistency across galaxy 
types extends down to —2 magnitudes below the bright end cutoff; at fainter magnitudes the star 
formation history of the population comes into play, causing deviations above or below the empirical 
formulation (Ciardullo 2010). In fact, the detailed shape of a galaxy's (or galactic subpopulation's) 
PNLF carries information about the stellar mass and age of the parent population, and can, for 
example, constrain models of merger scenarios (Bhattacharya et al. 2021). 


The PNLF works, but is not completely understood. It was surprising that models using relatively 
simple assumptions matched the early PNLF data well (Jacoby 1989; Ciardullo 2010), though the 
required distribution of central star masses (see below) was narrow and higher than the mean of 
known PN central stars or of white dwarfs. It continues to be true that decades after its formulation 
and after application to PNe in dozens of galaxies of different types, the PNLF works better than it 
should, given that it still lacks a firm theoretical underpinning. Ciardullo (2012) presents extensive 
discussion on this point, arguing that it should not work in the first place. 


Derivations by several investigators of the central star masses required to populate the bright 
end of the PNLF have sat stubbornly at or above ~ 0.6 Mo (Jacoby 1989; Marigo et al. 2004; 
Schónberner et al. 2007). The accelerated evolutionary tracks of M3B16 used in the PNLF models 
of Gesicki et al. (2018) help by allowing lower-mass stars to produce bright PNe. In their models 
the bright end of the PNLF is populated by PNe with progenitor masses from 1.1 — ~2 Mọ with 
lifetimes from 1-~7 Gyr. Recently, Galera-Rosillo et al. (2021) obtained deep spectrophotometric 
observations of a sample of bright PNe in M31’s extended disk. They found that models of the four 
Mx PNe in their sample indicated progenitor masses clustered at the “knee” of the 1.5 Mo M3B16 
track; these stars, according to Gesicki et al. (2018), can spend ~ 1000 years at the bright end of 
the PNLF. Their average remnant mass is 0.574 Mo, which, interestingly, agrees with the 0.58 Mo 
^maximum final mass" derived from PNLF simulations by Valenzuela et al. (2019). However, Davis 
et al. (2018) analyzed a sample of bulge PNe in M31 and concluded that after accounting for the 
dimming effects of circumstellar dust, the brightest objects at the observed PNLF cutoff are in 


24 For a discussion of a generalized form of the PNLF see Longobardi et al. (2013). 
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fact up to 1 mag brighter than the canonical Mx value (implying even more massive central stars 
20.66 Mo). These authors speculate that the constancy of the bright-end cutoff in the PNLF might 
be more indicative of higher circumstellar extinction with increasing central star luminosity than 
of a constant luminosity upper limit. Apparently, a comprehensive explanation for the observed 
consistency of the PNLF still eludes us. 


Though we still do not entirely understand its foundations, the bottom line is that the PNLF 
promises to remain a valuable tool. Studies like those by Gesicki et al. (2018), Valenzuela et al. 
(2019), and Bhattacharya et al. (2021) demonstrate the significant potential of the PNLF to study 
star formation history and stellar populations. As a distance indicator, the PNLF is benefitting from 
new instrumental capabilities and analysis techniques that provide increasingly accurate and deeper 
measurements of PNe in distant galaxies. For example, Roth et al. (2021) used data from the MUSE 
integral field spectrograph on the VLT to develop an improved method of measuring PN brightnesses 
that can extend the reach of the PNLF to ^40 Mpc, a distance where it would become an additional 
cosmological tool. 


6. CENTRAL STAR BINARITY 


Multiple lines of evidence point to binary systems as the birthplaces of many — if not most, 
or even all — PNe. Excellent and detailed discussions of the role of binaries in PN formation, 
morphology and evolution are found, for example, in De Marco (2009); Jones & Boffin (2017); De 
Marco & Izzard (2017); Frank et al. (2018); and Boffin & Jones (2019). The current tally?? of 
~60 confirmed short-period, ~10 long-period binary central stars, plus a few dozen systems with 
IR photometric excesses indicative of low-temperature companions, will surely grow as observational 
sensitivity to binary indications sharpens. For example, a search for new binary CSPN (Central 
Star[s] of a PN) candidates using data from TESS (Transiting Exoplanet Survey Satellite; Ricker et 
al. 2015), yielded seven out of eight targets showing some kind of photometric variation (Aller et 
al. 2020), though the authors note that three of these could result from contamination by a nearby 
variable source. In addition, the intrinsic variability of CSPN winds can be mistaken for evidence 
of binarity, so observed variations should be strictly periodic before binarity is claimed. Using Gaia 
epoch photometry Chornay et al. (2021) recovered several previously known binaries, and found 
several more objects consistent with binarity, demonstrating promise for future confirmations from 
Gaia data.?6 


6.1. Observational and Theoretical Predictions 


PNe are predicted as a result of common envelope (CE) evolution (De Marco 2009). Stars in close 
binaries can experience a short-lived period during which the more evolved of the pair expands and 
overflows its Roche lobe. The outflowing material increases in angular momentum at the expense of 
the binary orbit, which shrinks in response; the secondary star spirals in toward the primary before 
the primary’s envelope is ejected and eventually ionized by the exposed hot core. A related scenario, 
“grazing envelope evolution," may occur in systems where, prior to the secondary being engulfed, an 


?5 D. Jones’ website: http:/ /www.drdjones.net/bcspn/ 


?6 Despite the observed non-spherical morphologies and complex structures in PPNe, searches for close binaries have 
failed (Hrivnak et al. 2011, 2017), and the resulting limits on orbital parameters imply separations too large to have 
caused the observed features and compatible only with periods on the order of years. For example, Manick et al. 
(2021) recently discovered the longest-known period in a PPN binary: (IRAS 08005-2356, P~7.3 yr). Hillwig (2018) 
concludes that the population of currently-observed PPNe are not the precursors of PNe with close binary central 


stars. 
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accretion disk generates jets that efficiently remove mass from the envelope (Soker 2015), preventing, 
or at least delaying, the onset of a CE phase. In wider binaries (i.e., in non-contact binaries), hydro- 
dynamic simulations of AGB mass loss via “wind Roche-lobe overflow” (Mohamed & Podsiadlowski 
2012), show that part of the accelerated stellar wind slingshots around the accreting star toward the 
L2 Lagrange point, producing spiral-shaped density enhancements in the orbital plane resulting in 
decidedly non-spherical ejecta. Another possible mechanism for producing a spiral pattern lies in the 
orbital motion of the binary (Mastrodemos & Morris 1999; Mauron & Huggins 2006). Such spiral 
patterns have been observed in, e.g., the carbon star AFGL 3068 (Mauron & Huggins 2006) and the 
AGB star R Sculptoris (Maercker et al. 2012). 


Binaries may provide a natural explanation for the morphological and kinematical features seen 
in many PNe. Only ~20% of PNe can be described as having the basically spherical morphology 
expected for an envelope ejected and formed by interacting winds from a single star (Parker et al. 
2006) (see 87 on morphology). The markedly non-spherical morphology of the remaining 8096 is hard 
to explain as arising from a star that has never had a companion. For example, García-Segura et 
al. (2014) concluded that a single star cannot achieve the rotation rate necessary to form a bipolar 
PN, nor can it muster the magnetic field strength required to break spherical symmetry, whereas 
models of common envelope (CE) ejection naturally form bipolar nebulae. Furthermore, theory 
predicts that the symmetry axis of the resulting post-CE nebula should be oriented perpendicular 
to the binary orbital plane, which is precisely what has been observed in every close-binary PN for 
which the data are available (Hillwig et al. 2016; Munday et al. 2020). Kinematically, the outflow 
features in many PNe with binary central stars appear older than the central regions, which would 
be possible if these structures originated from an accretion disk prior to the general envelope ejection 
forming the PN proper (e.g., Corradi et al. 2011). Observed physical features such as point sym- 
metry (morphologically and kinematically similar structures on opposite sides of the central star), 
low-ionization microstructures, and collimated fast outflows (jets), likewise a challenge to single-star 
evolution, have been observed among a significant proportion of the OGLE sample of binary PN 
central stars (Miszalski et al. 2009b). However, the CE ejection scenario is not without flaws: it is 
unable to account for the very small masses (~ 0.01 Mo) measured for the nebulae, which would 
have been ejected over a short timescale during the post-CE phase, long after the main evolved 
stellar envelope has had time to diffuse away (Corradi et al. 2015b). Clearly there is more to be done 
to understand the relationship between CE binaries and PNe. 


Globular cluster PNe (Bond 2015), which are rare, (4 known, with another 4 candidates; see 85.1.1), 
may provide another argument for binary origin. The four known PNe are located in old globular 
clusters and have central star masses too high to have originated from a single star (Jacoby et al. 
2017). A PN ejected from a single current post-AGB remnant with a typical mass ~0.50 Mo in a 
globular cluster will have expanded and dissipated before the star has heated sufficiently to ionize 
the nebular gas. Alves et al. (2000) concluded that K648 (= Ps 1) in M15, the first PN discovered 
in a globular cluster, looks like a merger remnant with a mass ~0.6 Mo, higher than the typical 
remnant of a single AGB star. This conclusion was based on evolutionary tracks by Vassiliadis & 
Wood (1994) and Schónberner (1983). The M3B16 tracks, incorporating updated microphysics and 
improved modeling of TDU events on the thermally-pulsing AGB, lead to significantly shortened 
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timescales by as much as a factor of 10. These faster evolutionary rates make single stars once again 
tenable as globular cluster PN progenitors, leading to higher predicted numbers, but they appear 
inconsistent with the shorter kinematic ages of the known globular cluster PNe (Jacoby et al. 2017). 
The question is far from closed: globular clusters in other Local Group galaxies appear comparably 
deficient in PNe (Bond 2015), so whatever may be the eventual explanation for the low observed 
numbers in Milky Way globular clusters, it is likely applicable more broadly. 


6.2. Close Binaries: Chemical Abundances and the ADF 


Abundance discrepancies between analyses of permitted vs. collisionally-excited emission lines, 
where the former imply higher chemical abundances than the latter, are observed in many PNe. 
The “abundance discrepancy factor" (ADF) is the abundance ratio (typically O^* /H*) derived from 
recombination lines to those from collisionally-excited lines (see 83). Disagreements are found to be 
more extreme in PNe with close binary central stars, sometimes exceeding a value of 100 (Corradi et 
al. 2015a). A compilation of H II regions and PNe with measured ADFs is available on R. Wesson's 
website?" Wesson et al. (2018) posit that the observed PN ADF distribution in PNe (see Fig. 5.4 
in Boffin & Jones 2019) consists of two components: one a “baseline” effect (seen in H II regions as 
well) of a factor of a few, and a second that arises only (but not always) in PNe with close binary 
central stars. One proposed explanation for these extreme ADFs is the existence of spatially and 
chemically distinct gas phases in the nebula (Liu et al. 2000): ^normal" gas (emitting collisionally- 
excited lines) and cooler, H-deficient gas (emitting mostly recombination lines), possibly resulting 
from discrete ejection events (Corradi et al. 2014). There appears to be a threshold orbital period 
around 1.2 days, under which ADFs tend to be larger (75; Wesson et al. 2018), which presumably 
relates to the history and timing of ejection events leading to the extreme ADF values. In addition, 
several authors have shown that the emitting regions are also kinematically distinct, in many cases 
with the H-deficient gas not following the general outward acceleration seen in the ^normal" plasma, 
e.g., Richer et al. (2013); Pena et al. (2017). 


6.3. Evidence of Mass Transfer 


Since the first identification of lines from Kr and Se (Dinerstein 2001), many s-process nuclei such 
as Ge, Se, Br, Rb, Cd, Te and Xe, have been detected in PN gas. Their origin is straightforwardly 
understood, as described in 83.2. As for detection in PN central stars, there is an intriguing class 
called “barium stars,” a subset of which are giants, whose spectra show strong lines of Ba and other 
s-process nuclei plus carbon, even though the stars have not evolved sufficiently to produce these 
elements. It is now believed that barium stars are members of binaries, and that their enriched 
material was accreted via wind Roche Lobe overflow as previously described, when a now-faint 
companion white dwarf was an AGB star (Boffin & Jorissen 1988; Jorissen et al. 2019). We know of 
several PN central stars with cool subgiant or giant barium components in wide binaries (‘Tyndall et 
al. 2013; Boffin & Jones 2019, Table 4.2), providing compelling evidence for a mass-transfer episode 
prior to the PN phase. For close binaries, there was no such smoking gun for mass transfer until the 
discovery of a carbon-rich dwarf secondary star in the Necklace Nebula, a post- CE system with a 
binary period of a little over one day (Miszalski, Boffin, & Corradi 2013). So we see firm evidence 


27 https:/ /www.nebulousresearch.org/adfs/ 
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among both short- and long-period PN central star binaries for earlier episodes of mass transfer with 
the signature of TDU. 


6.4. Demographics of Binary Companions 


The first confirmed PN binary central star was UU Sge, a previously known variable identified by 
Bond et al. (1978) as the central star of Abell 63, a bipolar PN; this system has become an archetype 
for close-binary CSPN. UU Sge is an eclipsing binary with a period of 0.465 days (Bond et al. 1978), 
(more recently measured to be 0.456 days by Bell et al. (1994)); this is near the peak of the period 
distribution of known close binary CSPN (Jones & Boffin 2017, Fig. 4). With periods clustered 
around a few hours to a few days, the distribution exhibits a very sparsely populated longer-period 
tail. There are presently only four known long-period spectroscopic binaries with measured periods, 
all greater than 3 years (Boffin & Jones 2019). Discovery is biased toward shorter periods which are 
easier to detect, both photometrically and spectroscopically; De Marco et al. (2004) argue that a 
significant fraction of CSPN may be members of intermediate-period spectroscopic binaries. Binarity 
may not be the only relevant configuration; Soker (2016) and Bear & Soker (2017) posit that 1/6 to 
1/8 of PNe are “too messy” to have been formed in a binary, but require interacting triple stars (not 
wide triples like NGC 246; Adam & Mugrauer 2014) to explain their shapes (Soker 2021). 


Of the 60 or so known close-binary central stars, only a dozen or so have meaningful constraints 
on their components (Table 3.2 in Boffin & Jones 2019). Most are found to have late main sequence 
companions; two of the three double-degenerate (WD/WD) systems have total masses above the 
Chandrasekhar limit (albeit with large uncertainty), along with a merger time shorter than the age 
of the universe. These might be candidates for eventual Type Ia supernova events (see 87.5) . 


The total binary fraction of all PN central stars — a zealously sought grail with implications for 
both stellar and PN evolution — is unknown. Recent estimates of the PN post-CE close-binary 
fraction (periods less than about a day) from surveys based on data from OGLE and Kepler, cluster 
around ~10-20% (e.g., Miszalski et al. 2009a; De Marco et al. 2015; Jacoby 2019). But estimates 
of the total binary fraction calculated by various investigators do not appear to be converging on 
any kind of consensus at present (Boffin & Jones 2019). Interpretation and comparison of ob- 
served and predicted PN central star binary fractions is muddied by the fact that some merger 
remnants may not now be identifiable as having had a binary origin. And if most PNe do come 
from binaries, then the numbers detected in the Milky Way and other galaxies agree reasonably well 
with predictions from evolution of stellar populations in galaxies. Based on evolutionary models 
available at the time, Moe & De Marco (2006) predicted several times the observed number of 
PNe if single stars are their presumed progenitors; such an analysis needs to be redone with the 
new accelerated tracks of M3B16. Weidmann et al. (2020) find that roughly half (88 out of 175) 
of the central stars whose luminosities, temperatures and gravities they were able to derive imply 
masses and ages that are consistent with models of single star evolution using M3B16 tracks. To 
account for the number of PNe observed, Jacoby (2006; during discussion at IAU Symp. #234) 
posited a multi-channel origin whereby the predominantly non-spherical, bright PNe are the result of 
binary interactions, while single stars produce the less numerous and fainter, extended spherical PNe. 
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7. MORPHOLOGY AND EVOLUTION OF PNE 


There has been an abiding historical question of how an originally spherically-symmetric AGB star 
metamorphoses into a distinctly non-spherical PPN and ultimately, PN. However, that appears to 
have been the wrong question, since evidence has been accumulating that asymmetries may begin 
while the star is still on the AGB. MHD (magnetohydrodynamic) models by Pascoli (2020) indicate 
that a spherical-appearing AGB star conceals anisotropic structures that will break out once the 
fast stellar wind commences. Sahai & Trauger (1998) carried out an HST Ha imaging survey of 
young PNe and found that all of their objects exhibited “highly aspherical” morphology, leading 
them to conjecture that high-speed jets or collimated outflows (presumably from a close companion) 
during the late AGB or early post-AGB phase are responsible for shaping PNe?*. Guélin et al. 
(2018) investigated the three-dimensional envelope structure of the nearby carbon star IRC+10 216 
(~130 pc; Menten et al. 2012) using VLA millimeter-wave interferometry, concluding that mass loss 
modulated by the orbit of a low-mass binary companion might explain the regularly-spaced shells 
also observed in the optical (Leào et al. 2006), as hypothesized by Soker (1998). Decin et al. (2020) 
carried out ALMA observations of oxygen-rich AGB stars that also show evidence of asymmetrical 
features such as shells, disks, and spirals. 


Advances in imaging and spectroscopic capabilities have spurred progress in observing and mod- 
eling the formation and development of nebular structures. Some examples: Schmid et al. (2017) 
used the SPHERE adaptive optics system and ZIMPOL polarimeter at the VLT to image the binary 
star R Aqr, (consisting of a Mira variable and a hot companion), and obtained detailed images of its 
strong and complex jet. Using the NACO infrared adaptive optics system at the VLT, Kervella et 
al. (2014) discovered a dust disk surrounding Lə Puppis, a nearby (d=64 pc) AGB star. Kervella et 
al. (2015) then went on to use the SPHERE/ZIMPOL instrument at the VLT to image the disk in 
scattered optical light. A color composite of their polarimetric visible-light images (Fig. 17) shows 
the nearly edge-on disk, from which bipolar plumes emerge, along with indications of the primary and 
companion star. Using these data, Chen et al. (2016) performed three-dimensional hydrodynamical 
simulations and showed that these outflows could be due to a thermal pulse, by periastron passage of 
the companion, or even by planet ingestion. Subsequent modeling by Zou et al. (2020) demonstrated 
that it is possible to form young bipolar PPNe with a variety of shapes. Fig. 1 of Balick et al. 
(2019) shows scattered-light HST images of a dozen PPNe with “candle-shaped” lobes presumed to 
have been carved by collimated jets plowing into ambient slower AGB wind. And in a recent series 
of papers, Balick and colleagues developed hydrodynamic models for the mass-ejection histories of 
PPNe, including winds, lobes, jets, and other structures (see Balick et al. 2020). In that paper they 
extend their analysis to explore effects of toroidal magnetic fields; see their Fig. 1 for images of very 
compact PPNe whose knots and protrusions along the axis of symmetry appear to have been shaped 
by such magnetic fields. 


Our understanding of the menagerie of PN shapes has advanced both theoretically, with the 
development of evolution and hydrodynamical codes, and observationally, with the advent of mul- 


28 The morphological classification systems for PPNe (Sahai et al. 2007) and for young PNe (Sahai et al. 2011) cite lobes, 
waists, and haloes, presaging that for PNe in general (e.g., Balick & Frank 2002; Manchado et al. 2000), as will be 


seen below. 


45 


tiwavelength milli-arsecond resolution plus deep, wide imaging. Combining morphology with other 
characteristics can help us better understand nebular history and evolutionary pathways. Below we 
describe nebular classification, and then go on to outline ideas about nebular shaping and evolution. 


7.1. Morphology and Classification Systems 


Given the twin notions of a spherical AGB star expelling its outer atmosphere via winds originating 
near its surface, and its hot stellar remnant creating a canonical Strómgren sphere with isotropic 
thermal pressure, it is easy to understand why the archetypal mental model of a PN is round. 
However, a quick inspection of any compilation of modern PN images, e.g., Balick's catalog of HST 
images?, quickly dispels this misconception. PNe display morphologies encompassing elliptical, 
bipolar, multipolar, point-symmetric, multiple-shell, and yes, round shapes in a fraction (~ 20%) of 
cases. The origin of this variety is a vigorous area of current research, involving stellar evolution, 
binarity, gas dynamics and nebular astrophysics. Reviews of PN morphology can be found in Balick 
(1987), Manchado (2003), and Shaw (2012). 


Curtis (1918) was the first to classify PNe into morphological categories, motivated by the belief 
that *a collection of illustrations showing the forms assumed by the planetary nebulae would be 
of considerable value as contributing to theories of the structure and life histories of these bodies." 
He was not wrong. And while his categories include familiar descriptors like halo, annular, and 
ellipsoidal, he could only dream of the detailed images and complex features whose interpretations 
challenge us today. 


Investigation of PN morphology specifically as a consequence of shaping by current and previous 
stellar winds (see 87.2), was first undertaken by Balick (1987). Based on the new availability of 
high-dynamic-range CCD images, he devised a two-dimensional classification sequence for global 
morphology (as shown in his Fig. 15) along the dimensions of shape (round, elliptical, butterfly) and 
evolutionary state (early, middle, late). Balick & Frank (2002) describe the third shape class slightly 
differently, dividing what they call bipolar into butterfly and bilobed, (see their Fig. 1). In the former, 
the extended lobes converge into the nebular center, while in the latter they connect to a central 
elliptical or round component. Other classifications focus on symmetry properties of the PN (e.g., 
Manchado et al. 1996a,b) ranging from reflection symmetries across one or both axes or point-wise 
across the nebular center. Manchado et al. (2000) classified 255 northern PNe, settling on three 
main classes: round (2096), elliptical (6196) and bipolar (1396) — see their Fig. 2. Other classification 
systems have been described in Corradi & Schwarz (1995); Parker et al. (2006); Manchado (2004) 
and Sahai et al. (2011). Some examples of PN morphological classifications are shown in Figs. 18-20. 


Study of the interconnection between nebular kinematics and morphology has benefitted from 
the availability of high-resolution, spatially-resolved spectra. The kinematic catalog compiled by 
López et al. (2012b) has enabled such detailed analysis; e.g., López et al. (2012a) constructed three- 
dimensional models of the bipolar PNe Hb5 and K3-17, both of which exhibit evidence of more 
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recent additional lobe formation, possibly leading to future poly-polar structure. 


It is important to bear in mind that the apparent shape and extent of a PN is dependent on 
numerous factors, including how it is observed. Observed shapes are a two-dimensional projection, 
and apparently different morphologies might be due to orientation effects (e.g., Dalick & Frank 
2002); for example, a bipolar PN viewed nearly pole-on will look like a ring. Morphology is also 
wavelength-dependent; e.g., Kwok (2018) cites bipolar structure in several PNe seen in the infrared 
by Zhang & Kwok (2009) but not in the optical. Kinematics and expansion velocities observed in 
PNe also depend on the emission line used (Chu et al. 1984). 


Dynamic range limitations and field of view are further issues affecting how a PN is seen. As an 
obvious example, multiple outer shells and faint extended halos have been discovered surrounding 
many PNe, (e.g., Jewitt et al. 1986; Chu et al. 1987; Balick et al. 1992; Ramos-Larios et al. 2016; 
Guerrero et al. 2020b; Kastner et al. 2020), presumably representing earlier episodes of mass ejection. 
An outstanding example is NGC 6543, the Cat's Eye Nebula, seen in Fig 21 (Balick et al. 2001). The 
central nebula is surrounded by no fewer than nine ~evenly-spaced, concentric shells (projected as 
rings), extending out to about twice the main nebular radius. Encompassing all of this, extending out 
roughly five times farther still, with an angular radius of 2.5 arcminutes, is a patchy, but massive, halo 
(Fig. 22) comprised of material ejected earlier during the star's AGB evolution (Villaver, Manchado, 
& García-Segura 2004). Such extensive shell systems are seen in PPNe as well (Sahai et al. 1998; 
Balick et al. 2011). 


7.2. Nebular Shaping: The Interacting Winds Model 


The shape of a planetary nebula is the result of multiple factors, beginning with often-collimated 
momentum flux in the outflowing matter from the stellar atmosphere which alters as the star evolves. 
The physical environment of the evolving star also plays a significant role in terms of companion 
objects, as well as the ambient ISM. 


The kernel of current theories describing PN morphology is the interacting stellar winds (ISW) 
model, as originated by Kwok et al. (1978) and further detailed by Kwok et al. (1978), Volk & 
Kwok (1985) and many others. In this scenario, the faster, less dense stellar wind from the evolving 
central star catches up with and plows into slower gas ejected previously during the AGB phase, 
creating a compressed shell of high density observed as the PN. The essential anatomy of such an 
interstellar bubble was described by Castor, McCray, & Weaver (1975) and further elucidated by 
Weaver et al. (1977). Though their focus was on nebulae formed by winds from massive early-type 
stars blowing bubbles into the ambient ISM, the basic mechanisms and structure are similar. More 
recent PN-specific hydrodynamic simulations are discussed below in 87.3. 


A cartoon schematic of a typical established PN bubble is shown in Fig. 23. The volume nearest the 
star is filled with free-streaming stellar wind (gray) with velocity ~1000 km s^!, and mass-loss rate 
7-107? Mo yr! (e.g., Krticka, Kubát, & Krticková 2020). This is surrounded by shocked stellar wind 
(blue) created by an inward-moving (reverse) shock formed when the wind first reached the AGB 
material. This shocked gas can reach temperatures ~ 10’ K, capable of producing detectable x-ray 
emission (Kastner et al. 2012; Montez et al. 2015), which has been observed in several PNe, including 


47 


NGC 6543 (Chu et al. 2001), as seen in Fig. 21, as well as the young PN BD+30° 3639 (Kreysing et 
al. 1992; Arnaud, Borkowski, & Harrington 1996), NGC 7009 (Fig. 27), and others listed in Table 1 
of Kastner et al. (2008). One thing these PNe always have in common is an intact, continuous inner 
rim that presumably confines the shocked wind. Radiation-hydrodynamic simulations of these hot 
bubbles and their x-ray emission have been carried out by several groups. Steffen, Schonberner, & 
Warmuth (2008) and Sandin et al. (2016) demonstrated the importance of thermal conduction at 
the interface between the shocked wind and the nebular gas. Toalá & Arthur (2014, 2016, 2018) 
performed high-resolution two-dimensional calculations following the evolution of these hot bubbles. 


The next shell out is composed of the ionized former AGB slow wind (red) that we see as the 
PN, with a temperature ~10* K. Between the hot shocked wind and the cooler nebular gas, a 
conduction layer (mentioned above) at ~10° K, is expected to produce collisionally-excited O VI 
emission, which has been observed by Ruiz et al. (2013). If the PN is radiation-bounded rather than 
matter-bounded, neutral swept-up material (tan) will lie outside an ionization front, and beyond 
that, the outward-moving leading shock demarcates the as yet unaffected remnant AGB wind (white). 


Since the nebula expands much faster than the sound speed in the displaced neutral gas, the 
original ISW formulation assumed momentum conservation, where the kinetic energy in the fast 
wind is radiated away and does not contribute significantly to the acceleration and compression of 
the swept-up material. It was able to explain some, but not all observed PN expansion velocities 
and masses (Kwok 2000; Balick & Frank 2002). Subsequently, Kwok (1982, 1983); Kahn (1983) and 
Volk & Kwok (1985) examined the energy-conserving case, and estimated that typically ~30% of 
the fast wind's energy is dumped into the expanding shell, which is then able to adequately explain 
PN expansion velocities. Detailed analysis of the dynamics and energetics of such systems can be 
found in the radiation-hydrodynamic models referred to in 87.3 below. 


In the original ISW, the general geometry is taken to be spherically symmetric, but the abundance of 
aspherical PN shapes requires that this assumption be relaxed. In the generalized ISW scheme (GISW; 
Balick 1987), the fast stellar wind expands into a pre-existing asymmetrical density environment, 
usually framed as a dense equatorial torus, resulting in fast lobes directed along the axis of the 
torus. The origin of this “waist” is not yet completely understood. Rotation in a single star is an 
unlikely source (Frank et al. 2018). It may result from extra momentum supplied by a companion 
in a post-CE binary system (Lagadec 2018), or perhaps from asymmetric AGB mass loss in a single 
star, mediated by the stellar magnetic field (García-Segura et al. 2005; Balick et al. 2020). Recently, 
Kastner et al. (2020) constructed UV—-to-IR. panchromatic HST overlay images of the young PNe 
NGC 7027 (Fig. 24) and NGC 6302 (Fig. 25) to study and illustrate PN shaping processes. The new 
NGC 7027 image reveals the intricacy of the dust filaments seen against the bright nebular shell, as 
well as the wispy, concentric dust rings. For the bipolar NGC 6302, the most conspicuous features 
of the new image are the red point-symmetric arcs seen in [Fe II] emission, which is an indicator of 
fast shocks. 


7.2.1. Internal Structures. 


About half of PNe, across almost all morphologies, display internal small-scale structures described 
by various terms including filaments, knots, jets, caps, or ansae (Balick 1987; Corradi et al. 1996; 
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Gongalves et al. 2001). These structures are best visualized in the light of low-ionization species 
like [O I], [N II], [O II], and [S II], hence they are called “LIS,” for low-ionization structures. LIS 
can manifest in various configurations and display velocities in a wide range from a few to hundreds 
of km s~!, indicating different formation and excitation mechanisms (Goncalves et al. 2001; Akras, 
Gongalves, & Ramos-Larios 2017; Balick et al. 2020). 


NGC 7293 (the Helix Nebula) exhibits a multitude of cometary knots (Fig. 26): dense, neutral 
clumps of material in the nebular interior, ionized on their star-facing side (O’Dell & Handron 1996). 
Such knots are also seen in other nearby PNe, e.g., NGC 6853 and NGC 6720 (O’Dell, Balick & 
Hajian 2002). O’Dell, Henney, & Ferland (2007) found that for NGC 7293 the central star’s radiation 
is sufficient to power the near-IR H5 emission observed from the knots, eliminating the need to appeal 
to shock excitation. In addition, Matsuura et al. (2009) found that the H, associated with the knots 
is likely a relic of molecular gas formed while the star was on the AGB, and not formed in situ (see 
83.5 on the molecular content of PNe). 


Balick (1987) noted the presence of low-ionization “inclusions” emitting prominently in the [N H] 
A6584 line in ~ 20% of the bright PNe he studied. Balick et al. (1994) dubbed these regions 
“FLIERs” for fast, low-ionization emission regions due to their supersonic relative velocities (7450 
km s 1). They tend to occur in pairs along the nebula's symmetry axis on opposite sides from the 
center. Examples of PNe with FLIERs are NGC 6826 and NGC 7009 (both in Fig. 27), and also 
NGC 3242 and NGC 6543 (Fig. 21). The origin of FLIERs has been difficult to determine. Balick 
et al. (2020) and García-Segura, Taam, & Ricker (2020) conclude that they can be explained as 
the remnants of slow, dense, collimated jets shaped by a toroidal magnetic field during the PPN 
phase (akin to the bullets seen in CRL 618; Balick et al. 2013), manifesting as shocked gas at the 
tips of the jets. Guerrero et al. (2020a) performed a statistical analysis of jets (inclusive of all the 
different names for these features) in 58 PNe. They found a bimodal space velocity distribution, with 
means of 66 km s^! (70%; low-speed) and 180 km s^! (3096; high-speed). The low-speed jets are a 
challenge to explain, and are perhaps related to a lighter outflow more vulnerable to deceleration via 
interaction with the slow AGB wind. 


'The data regarding magnetic fields in PNe are, unfortunately, very sparse. Observations of dust 
polarization have revealed toroidal fields in only three PNe: NGC 7027, NGC 6537, and NGC 6302 
(Sabin, Zijlstra, & Greaves 2007). Zeeman splitting in molecular maser lines is also a probe of 
magnetic fields (Hou & Gou 2020); thus far only four PNe exhibit distinct or possible pairs of 
Zeeman lines: IRAS 1633-4807, K3-35 (see footnote 22), IRAS 17393-2727, and JaSt 23 (Qiao et al. 
2016). 


1.2.2. Binary Effects. 


Binary companions can have profound effects on the nebular morphology: post-CE PNe are in- 
variably bipolar, and have their orbital plane oriented perpendicular to the nebular symmetry axis 
(see 86.1). Soker et al. (1998) presented arguments to explain the correlation between bipolarity 
and higher-mass progenitors. Frank et al. (2018) modeled the interaction of fast stellar wind with 
previous CE ejecta and were able to recreate bipolar morphology. Simulations of post-CE systems 
by García-Segura, Ricker, & Taam (2018) likewise reproduced bipolar morphology for young PNe, 
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with later transitions to barrel or elliptical shapes, sometimes including nested lobes. Numerical 
simulations of AGB winds in interacting binary systems (Bermüdez-Bustamante et al. 2020) can 
serve as the starting point for hydrodynamic simulations of CE evolution and eventual PN formation. 


The PN KjPn8 (López, Vázquez, & Rodríguez 1995; López et al. 2000) exhibits precessing 
“BRETs,” bipolar rotating episodic jets that have carved out the most extensive bipolar lobes seen 
in PNe, more than 4 pc end to end. Kwok (2018) suggests that such bipolar/multipolar lobes 
represent low-density cavities formed after outflows breach the dense equatorial torus. López et al. 
(2002) describe KjPn8 as a “double planetary nebula,” that could be the result of a binary system of 
two relatively massive stars evolving closely in time, experiencing two separate PN-creating events, 
each with its own symmetry axis. Another BRET nebula, Fg 1 (López, Meaburn, & Palmer 1993), 
displaying spectacular S-shaped point-symmetric jets (Fig. 28), was found by Jones et al. (2014) 
to contain a post-CE double-degenerate nucleus with a period of ~1.2 days, strengthening the 
association of these outflows with binary evolution. Additional evidence is provided by Guerrero et 
al. (2021), who have for the first time directly imaged the jets in NGC 2392, S-shaped like those in 
Fg 1, and which they deduce are being actively launched now, rather than representing remnants of 
earlier ejections. Also, based on the observed hard x-ray stellar emission and high level of nebular 
excitation, they conclude that the companion to NGC 2392's central star is a white dwarf, adding to 
the short list of double-degenerate systems in PNe. Longer-period systems with periods between 100 
and 1000 days, including “barium star" PNe (86.3) like WeBo 1 (Bond, Pollacco, & Webbink 2003) 
and A70 (Miszalski et al. 2012), exhibit similar morphologies, with equatorially-enhanced densities 
and bipolar extensions. Outflow speeds in bipolar PNe appear to be "tapered," i.e., moving slowest 
at low nebular latitudes, and fastest in the polar direction (e.g., Corradi 2004). The gas impacted by 
these outflows is often observed to exhibit homologous expansion, moving at velocities that increase 
with distance from the star (e.g., Corradi et al. 2001). 


Planetary-size bodies can also influence nebular morphology by enhancing a star's mass-loss rate 
on the RGB and AGB, and if a planet spirals into the AGB envelope, it can spin up the envelope 
causing asymmetrical mass loss (Soker 2018). In their modeling of the evolution of several exoplanet 
systems, Hegazi, Bear, & Soker (2020) concluded that a planet engulfed into an AGB envelope could 
lead to an elliptical PN, as posited by De Marco & Soker (2011). 


7.3. Expansion and Modeling of PN Evolution 


PNe have been known to be expanding at least since Zanstra (1932) noted as evidence the Doppler 
doubling of nebular emission lines.? Catalogs of nebular expansion velocities include Sabbadin 
(1984); Weinberger (1989), and more recently, the PN kinematic catalog of López et al. (2012b). The 
GISW model provides the phenomenological underpinning describing the response of ejected and 
surrounding material to the changing winds and radiation field from the star. The details of nebular 
expansion and evolution have been explored in depth by Schónberner and colleagues (Schónberner, 
Jacob, & Steffen 2005; Schónberner et al. 2005b, 2007, 2014, 2018), who used spherically-symmetric 


30 Earlier papers exploring the nature of doubled lines and at least flirting with the idea of expansion are those by 
Campbell & Moore (1918), Perrine (1929), and Zanstra himself (Zanstra 1931). Thanks to A. Zijlstra for these 
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one-dimensional radiation-hydrodynamics simulations to examine nebular behavior. 


Schónberner et al. (2014) detail their high-resolution echelle observations and models of various 
nebular components (see their Fig. 1) over time as the shock fronts produced by ionization and stellar 
wind advance. Their generated plots of variables including density, surface brightness, velocity and 
line profiles as a function of time allow comparison with observations. Among their conclusions 
is that it is the propagation of the outer shell's shock (outer edge of the tan ring in Fig. 23) that 
represents the PN's actual expansion velocity. This finding implies that the typical measured Doppler 
velocity (traditionally taken to be ~20 km s~', measured as half the line splitting in [O III] A5007, or 
half the half-width if unresolved)?! must be augmented by a factor ~1.3-1.5 (Jacob, Schónberner, & 
Steffen 2013; Schónberner et al. 2018) to more accurately capture the actual expansion rate. Jacob, 
Schónberner, & Steffen (2013) calculate a mean true expansion velocity of 42 + 10 km s^ for their 
sample of 59 PNe within a distance of about 2 kpc. GS19 find a similar value, 38 + 16 km s^, 
for their sample of 45 objects in Gaia DR2; GS21 find 32 + 13 km s^! for their EDR3 sample of 
65 PNe, all three samples agreeing within the uncertainties. As a consequence of larger expansion 
velocities, PN distances determined by the expansion parallax method will likewise be pushed to 
larger values (see 85.1.1). Combined with ongoing high-quality imaging and spectroscopic studies of 
PNe, meticulous investigations such as those described promise increased understanding of nebular 
evolution. 


7.4. Ages and Sizes of PNe 


Together with their distances, observed PN expansion velocities allow us to estimate kinematic 
ages (sometimes called *dynamical ages"). GS19's analysis of PNe in DR2 also included age es- 
timates for the sample of 45 PNe central stars for which they have nebular expansion velocities. 
They found that most of their objects were under 15,000 years old, though several appear more 
than 50,000 years old (see their Fig. 10). The average age was found to be ~24,000 yr; how- 
ever, the authors acknowledge a likely bias toward smaller ages, given the easier detectability of 
younger PNe. The larger EDR3 sample of GS21 likewise exhibits a majority of young PNe («10,000 
yr); the average age is 17,800 + 3000 yr. Strongly bipolar PNe appear to tend younger: from their 
3D models of a sample of bipolar PNe, Gesicki et al. (2016) derived ages of only a few thousand years. 


The physical radii of PNe are straightforwardly calculable from their distances and angular radii. 
Fig. 1 is a montage of 22 PNe with well determined distances, arranged according to relative physical 
size; the scale bar indicates 5 light years (~1.5 pc). The radii of the PNe shown range from about 
0.03 — 0.8 pc. Fig. 8 in GS19 shows the distribution of radii for 206 PNe, calculated with DR2 
distances and angular sizes (10% Ha isophotes) from the HASH database (Parker et al. 2016). More 
than half are under 0.2 pc in radius, though the distribution has a long tail beyond 1 pc with small 
numbers of larger PNe up to a maximum radius of almost 2.2 pc (the PN EGB 6; Ellis, Grayson, 
& Bond 1984). In the GS21 sample, 54% of the objects have nebular radii «0.3 pc, while 14% are 
greater than 1 pc in radius (see their Fig. 12). The mean radius of their sample is 0.59 pc. 


7.5. Late Evolution of PNe 


31 This is a simplistic view, at best, since it is well known that velocities depend on the ion observed, as well as the 
position in the nebula (e.g., Chu et al. 1984; Gesicki & Zijlstra 2000). 
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As a PN and its central star evolve, the nebular density decreases and the surface brightness 
declines toward invisibility. Morphologically, lobes become washed out and the nebula as a whole 
tends toward becoming more spherical. Pereyra, Richer, & López (2014) studied the kinematics of 
evolved PNe, finding, not unexpectedly, that the most evolved objects expand most slowly. 


Some PNe have shapes that reveal expansion into an inhomogeneous ISM, denser in some directions 
than others.? Studies, both observational (e.g., Tweedy & Kwitter 1996; Xilouris, Papamastorakis, 
& Paleologou 1996; Kerber et al. 2000) and theoretical (e.g., Soker, Borkowski, & Sarazin 1991) 
have been carried out, some including interaction with previously-ejected AGB material as well (e.g., 
Villaver, García-Segura, & Manchado 2003; Villaver, Manchado, & García-Segura 2012; Villaver 
et al. 2014). In their Atlas of Ancient PNe, Tweedy & Kwitter (1996) found 21 out of 27 objects 
exhibiting evidence of ISM interaction. A striking instance of PN-ISM interaction shaping is seen 
in Sh2-188 (Fig. 29); the central star is clearly offset from the nebular center toward the brighter 
region in the southeast. This PN was analyzed in detail by Wareing et al. (2006), whose model, 
incorporating stellar winds plus the effect of the nebula's motion through the ISM, was able to 
reproduce the nebula’s shape. They found that the nebula's peculiar velocity (with respect to the 
local ISM) is ~125 km s~}, directed toward the brightest part of the limb in the southeast, consistent 
with the proper motion of the central star. Wareing (2010) describes the effect as “rebrightening” 
when a PN shell moving through the ISM interacts with a pre-existing bow shock formed earlier in 
the star's evolution. Abell 21 offers another example, as shown in Fig. 30. Like Sh2-188, it appears 
much brighter on one side and the central star is also clearly offset from the nebular center. One more 
object worthy of mention here is Sh2-216, the closest PN (d=129 pc; Harris et al. 2007). It is also 
one of the largest currently known, angularly as well as physically: 1.6°=3.8 pc in diameter. Tweedy, 
Martos, & Noriega-Crespo (1995) constructed deep mosaics in Ho and [N II] that demonstrate that 
the bright eastern rim, where ISM interaction appears to be taking place, is part of a much fainter, 
nearly round, nebula (see their Figs. 1-4). The central star, a hot DAO white dwarf (Tweedy & 
Napiwotzki 1992), is displaced from the nebula center, and has a proper motion consistent with such 
an interpretation (Cudworth & Reynolds 1986). 


Intriguingly, there is an example of a classical nova occurring inside a PN. When V458 Vul erupted 
in 2007, the outburst flash-ionized a pre-existing bipolar nebulosity that Wesson et al. (2008) con- 
clude is a ~14,000 yr-old PN produced by CE evolution of the binary central star (see 86.1). Nova 
Persei 1901 (GK Per) has been suggested as another example of a nova inside an ancient PN (Bode 
et al. 2004). 


Finally, the fascinating possibility of a Type Ia supernova (thermonuclear white dwarf cataclysm) 
exploding inside a pre-existing PN has been explored by Tsebrenko & Soker (2015a,b). They estimate 
this to occur in at least ~20% of Type Ia events, those resulting from a binary merger of a white 
dwarf with an AGB core or with another white dwarf. Their three-dimensional hydrodynamical 
simulations can explain the observed shapes of some supernova remnants (SNRs), in particular those 
resembling elliptical PNe with “ears,” small, incomplete lobes on opposite sides of the center, like Ke- 


3? Note that the ISM surrounding the PN is not expanding. 
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pler's SNR (Tsebrenko & Soker 2013) and G1.9+0.3 (Tsebrenko & Soker 2015b). A prime candidate 
for such a Type Ia merger event is the halo PN TS 01 (PN G135.94-55.9), whose double-degenerate 
binary central star has a total mass very close to the Chandrasekhar Limit (Tovmassian et al. 2010). 
Another potential candidate might be V458 Vul, mentioned above, which is thought to contain two 
compact stars with a total mass above the Chandrasekhar Limit (Wesson et al. 2008; Rodríguez-Gil 
et al. 2010; Corradi 2012). 


8. SUMMARY AND OUTLOOK 


We examined PNe in the context of stellar evolution in 82, as the penultimate stage in the lives 
of low- and intermediate-mass stars. In 83 we laid out the current landscape of chemical abundance 
determinations, covering PN abundances in the Milky Way, including disk, bulge and halo samples. 
We discussed neutron-capture elements, the less-abundant elements fluorine, magnesium, iron and 
zinc, as well as the molecular and dust content of PNe. We then described abundance studies of 
PNe in more than a dozen external galaxies in and beyond the Local Group. In 84, we presented an 
analysis of radial abundance gradients in the Milky Way and external galaxies. Next, the observed 
elemental abundances in PNe were compared with those predicted from theoretical models. In 
85, we described techniques for determining distances to PNe, and discussed the Planetary Nebula 
Luminosity Function. Central star binarity and its consequences were considered in 86. Nebular 
morphology, shaping, and evolution were explored in 87. 


Our understanding of PNe has advanced significantly, but numerous open questions remain. Areas 
ripe for future investigations have been pointed out by, e.g., Kwitter et al. (2014), including: ex- 
ploiting non-optical wavelengths for detection of new PN candidates; better characterization of AGB 
mass loss, especially as it relates to metallicity; improved atomic data for abundance determinations, 
especially in the case of s-process elements; continued efforts at resolving the ADF; the formation 
and fate of dust and molecules around PNe; and improved understanding of the effects of binarity. 
Existing observing facilities (e.g., ALMA, Gemini, the VLT and GTC) continue to provide vital 
data to address such issues. Future telescopes and instrumentation promise to extend these capa- 
bilities. Guerrero (2020) discusses future x-ray missions, including eROSITA, which will produce 
an all-sky survey with a sensitivity capable of detecting the hot gas in the interior of hundreds of 
PNe, thereby improving the statistics of x-ray emitting PNe. In the near- to mid-IR JWST will 
enable detection of lower-mass companions in PN central star binaries and dusty disks in PNe, along 
with high-resolution mapping of UV-irradiated knots to probe the evolution of molecules and dust 
grains in the PN environment (Sahai 2020). Lagadec (2020) anticipates contributions from the ELT's 
first instruments, imager/spectrometers and IFUs, that will map dust and molecule distributions, 
radio continuum surveys like ASKAP*? will enable determination of accurate angular diameters for 
PNe, including compact objects as well as those suffering heavy optical extinction. Finally, the 
capabilities of the ngVLA will allow determinations of light-element abundances, measurement of 
precise nebular proper motions, and studies of close binary interactions in PN central stars (Kastner 
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et al. 2018), as well as investigations of molecular lines to study the AGB-PN transition. 


In this review/tutorial we have seen how the observed properties of planetary nebulae disclose the 
workings of stellar evolution, nucleosynthesis, and time. It is fitting to end with a quote (Aller 1971) 
from the late Lawrence H. Aller, renowned scholar of planetary nebulae and mentor to one of us 
(KBK): "They are wreaths placed by Nature around dying stars." 


APPENDIX 


A. A LISTING OF SOME PN CATALOGS AND DATABASES 
A.1. Nebulae 
The Strasbourg-ESO Catalogue of Galactic Planetary Nebulae. Parts I, II 


Acker, A., Marcout, J., Ochsenbein, F., Stenholm, B., Tylenda, R., Schohn, C., 1992 European 
Southern Observatory, Garching (Germany), 1992 


Planetary Nebula Image Catalog 
Balick, B.; http://faculty. washington. edu/balick/PNIC/ 


New Galactic Planetary nebulae selected by radio and multiwavelength characteristics 
Fragkou, V., Parker, Q.A., Bojicic, I.S., & Aksaker, N. 2018, MNRAS, 480, 2916 


Planetary nebulae in the UWISH2 survey 
Gledhill, T.M., Froebrich, D., Campbell-White, J., & Jones, A.M. 2018, MNRAS, 479, 3759 


An Atlas of Images of Planetary Nebulae 
Górny, S.K., Schwarz, H.E., Corradi, R.L.M. & Van Winckel, H. 1999, A&AS, 136, 145 


The coordinated radio and infrared survey for high-mass star formation— IV. A new 
radio-selected sample of compact galactic planetary nebulae 
Irabor, T., Hoare, M.G., Oudmaijer, R.D., et al.2018, MNRAS, 480, 2423 


A. Catalog of Relative Emission Line Intensities Observed in Planetary Nebulae and 
Diffuse Nebulae 
Kaler, J.B. 1976, ApJS, 31, 517 


Version 2000 of the Catalogue of Galactic Planetary Nebulae 
Kohoutek, L. 2000, A&A, 378, 843 


Gallery of Planetary Nebula Spectra 
Kwitter, K. B. & Henry, R.B.C.; tinyurl.com/63ed7txz 
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Catalog of Intensities, Analysis & Abundances in Galactic Planetary Nebulae 
Kwitter, K. B. & Henry, R.B.C.; tinyurl.com/PN-analysis 


The IAC Morphological Catalog of Northern Galactic Planetary Nebulae 
Manchado, A., Guerrero, M.A., Stanghellini, l, & Serra-Ricart, M. 1997, IAU Symp. No.180, 24 


The Macquarie/ AAO/Strasbourg Ho Planetary Nebula Catalogue: MASH 
Parker, Q.A., Acker, A., Frew, D.J., et al. 2006, MNRAS, 373, 79 


HASH: the Hong Kong/AAO/Strasbourg Ho planetary nebula database 
Parker, Q.A., Ivan S Bojici¢, I.S., & and Frew, D.J. 2016 J. Phys: Conf. Ser. 728 032008; 
http://vizier.u-strasbg.fr/vizier/MASH/ 


A catalogue of 108 extended planetary nebulae observed by GALEX 
Pradhan, A.C., Panda, S., Parthasarathy, M., Murthy, J., & Ojha, D.K. 2019, Ap&SS. 364, 181; 
https://arxiv.org/pdf/1910.05543.pdf 


The SPM Kinematic Catalogue of Planetary Nebulae & Extragalactic Planetary Neb- 
ulae 
Richer, M. G., López, J. A., Díaz-Méndez, E., et al. 2010, RMxAA, 46, 191 
http://kincatpn.astrosen.unam.mza 


First deep images catalogue of extended IPHAS PNe 
Sabin, L., Guerrero, M.A., Ramos-Larios, G., Zijlstra, A.A., Awang Iskandar, D.N.F. 2021, 
arXiv:2108.13612 


First release of the IPHAS catalogue of new extended planetary nebulae 
Sabin, L., Parker, Q.A., Corradi, R.L.M., et al. 2014, MNRAS, 443, 3388 


An evolutionary catalogue of galactic post- AGB and related objects 
Szczerba, R., Siódmiak, N., Stasińska, G., & J. Borkowski, J. 2007, A&A, 469, 799 
https://www.ncac.torun.pl/postagb 


The Fornax3D project: Planetary nebulae catalogue and independent distance mea- 
surements to Fornax cluster galaxies 
Spriggs, T.W., Sarzi, M., Galán-de-Anta, P.M., et al. 2021, A&A, in press 


An imaging and spectroscopic study of the planetary nebulae in NGC 5128 (Centaurus 
A) 
Walsh, J.R., Rejkuba, M., & Walton, N.A. 2015, A&A, 574, A109 
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A.2. Central Stars 


VizieR Online Data Catalog: Central stars of planetary nebulae in Gaia DR2 
Chornay, N & Walton, N.A. 2020, A&A, 638, 103 


One star, two star, red star, blue star: an updated planetary nebula central star 
distance catalogue from Gaia EDR3 
Chornay, N & Walton, N.A. 2021, arXiv:2102.13654v1 


All Known Binary CSPN 
DeMarco 2009 PASP, 121, 316 


Planetary Nebulae in Gaia EDR3: Central Star identification, properties and binarity 
González-Santamaría, L, Manteiga, M., Manchado, A., Ulla, A., Dafonte, C., & López Varela, 
arXiv:2109.12114v2 


Binary central stars of planetary nebulae identified with Kepler/K2 Jacoby, J., Hillwig, 
T., Jones, D., Martin, K., DeMarco, O., Kronberger, M., Horowitz, J., Crocker, A., & Dey, J., 2021, 
MNRAS, 506, 5223 


Post-CE Binary Central Stars of PNe 
D. Jones; https://www.drdjones.net/bcspn/ 


Catalogue of central stars of planetary nebulae: Expanded edition 
Weidmann, W.A., Mari, M.B., Schmidt, E.O., Gaspar, G., Miller Bertolami, M.M., Oio, A.G., 
Gutiérrez-Soto, L.A., Volpe, M.G., Gamen, R., & Mast, D. 2020, A&A, 640, A10 


B. MODEL-PREDICTED YIELDS AND SURFACE ABUNDANCES OF LOW AND 
INTERMEDIATE-MASS STARS 


A stellar element yield is the total amount of that element that is synthesized by the progenitor 
star during its lifetime and ultimately expelled into the ISM. Examples of early estimates of LIMS 
yields can be found in Renzini & Voli (1981), van den Hoek & Groenewegen (1997), Buell (1997), 
and Marigo (2001). Karakas & Lattanzio (2014, Table 2) provide a listing of 16 sources of model- 
predicted yields and surface abundances published between 2004 and 2014. For each reference, the 
stellar mass and metallicity ranges are noted as well as whether or not s-process predictions are 
included. More recent papers containing AGB yields and surface abundance information include 
Miller Bertolami (2016), Karakas & Lugaro (2016), Karakas et al. (2018), Ventura et al. (2018), 
Marini et al. (2021), and numerous results relevant to the model output provided at the FRUITY 
website (http: / /fruity.oa-teramo.inaf.it/modelli.pl). Each reference provides details and specifications 
regarding the models along with the assumptions that went into them. While we shall not review or 
critique these papers here, we urge the interested reader to delve into them in order to appreciate the 
differences in treatments of processes such as extra mixing, mass loss, third dredge-up, hot bottom 
burning, and formation of C pockets, just to mention a few. 
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Table 1. MW Planetary Nebula Abundance Surveys: Sample Sizes, Spectral Ranges, & Abundance 
Medians 


Source Number A Range (A) He/H cjuo N/H O/H Ne/H S/H Cl/H Ar/H 
MW Disk 

Barker78 32 3400-7400 0.106 ee 3.63E-05  2.63E-04 5.02E-05  2.43E-06 

Aller83--87 72 3400-7800 0.107  7.28E-04  1.11E-04  4.07E-04 8.32E-05  8.91E-06  1.51E-07  2.29E-06 
Kingsburgh94 58 3150-7400; 0.112 3.98E-04  1.74E-04  4.58E-04  1.10E-04 7.38E-06 ee 1.85E-06 
Perinotto04 © 114 Optical+UV 0.106 t 1.34E-04  4.27E-04  1.01E-04  4.96E-06 t 1.59E-06 
Stanghellinio6 © 78 Optical 0.114 Ue 1.50E-04  3.25E-04  8.20E-05 Zn Ue 8.25E-07 
Girard07 37 3700-7500 0.106 Ue 1.32E-04  4.90E-04  9.34E-05  8.13E-06  1.03E-07  2.32E-06 
Pottasch10€ 32 Op+IR+UV 0.107  4.03E-04  1.50E-04  4.15E-04  1.31E-04 7.20E-06 Ue 2.70E-06 
Maciel17 230 3600-7900 0.112 5.37E-04  1.51E-04 4.47E-04 8.91E-05 7.41E-06 tee 2.51E-06 
KH20 152 Optical+UV 0.121 7.71E-04  1.67E-04  3.95E-04  1.02E-04  4.59E-06  8.91E-08 2.16E-06 

MW Bulge 

Aller87 15 3400-7400 0.107 Ue 2.69E-04  4.27E-04  7.94E-05 7.93E-06  4.02E-07  3.16E-06 
Webster88 33 3700-7350 0.120 Ue 1.80E-04  5.45E-04 $e Mice oe je 
Ratag97 20 3500-7600 0.106 Ue 1.32E-04  4.42E-04  6.76E-05  6.31E-06  2.30E-06  3.13E-07 
Stasiriska98 © 85 Optical 0.100 Ue 1.60E-04  4.00E-04  9.60E-05 Ue Ue 
Exter04 32 3600-7400 0.112 2.47E-04  4.43E-04  8.41E-05 7.06E-06 Ue 3.43E-06 
Wang07 20 Op+IR+UV 0.115 2.99E-04 1.78E-04 5.25E-4 1.2AE-4 1.18E-05  2.14E-07 1.82E-06 
Chiappini09 167 3500-7400 0.122 ee 1.22E-04  3.37E-04 7.78E-05 6.28E-06 ee 2.00E-06 
Cavichial0+17 37 3600-7900 0.110 Ue 1.05E-04  2.29E-04  5.96E-05  3.80E-06 Ue 1.48E-06 
Solar 0.085 2.69E-04  6.76E-05  4.90E-04  8.51E-05  1.32E-05  3.16E-07 2.51E-06 


4 Source references: Aller83=Aller & Czyzak (1983), Aller87—Aller & Keyes (1987), Barker78=Barker (1978b), Cavichial0=Cavichia et 
al. (2010), Cavichial7=Cavichia et al. (2017), Chiappini09=Chiappini et al. (2009), Exter04—Exter et al. (2004), Girard07=Girard, 
Köppen, & Acker (2007), Kingsburgh94=Kingsburgh & Barlow (1994), KH20=Kwitter & Henry (2020), Maciel17—Maciel et al. 
(2017), Perinotto04=Perinotto, Morbidelli, & Scatarzi (2004), Pottasch10=Pottasch & Bernard-Salas (2010), Ratag97—Ratag et 
al. (1997), Stanghellini06=Stanghellini et al. (2006), Stasiriska98=Stasiriska et al. (1998), Wang07—Wang & Liu (2007), Web- 
ster88— Webster (1988). 


b The Aller83+87 and KH20 samples comprise a mix of PNe whose carbon abundances were measured using either optical recombination 
lines or UV collisionally excited lines. The remainder of the samples used strictly UV collisionally excited lines. 


©The abundances in this study were computed by the authors but were based either partially or entirely upon spectrophotometric 
data from other researchers. 


Table 2. MW Halo Planetary Nebula Abundance Studies 
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PN: Common Name (PNG) Source He/H cn N/H O/H Ne/H S/H CI/H Ar/H 
BoBn1 (108.4-76.1) c Hawley78 0.111 ee 3.30E-04 . T.70E-05 . 5.20E-05 muss 
Barker80 0.115 Z 2.19E-04 7.59E-05 5.25E-05 3.89E-08 
TP81 0.095  1.20E-03* 2.19E-04 7.94E-05  1.00E-04 tee Ue Ue 
Pena91 0.105 — 1.45E-03* 8.71E-05 4.79E-05 = 5.75E-05 6.31E-07 Ue 5.50E-08 


Kniazev08 0.100 — 1.58bE-03* 4.3TE-05 | 6.46E-05 | 8.13E-05 1.45E-07 . 1.38E-09 . 3.72E-08 


Milanova09 0.129 3.02E-05 . 3.09E-04 


Otsukal0 0.118 1.05E-03 1.07E-04 5.51E-05 . 9.04E-05 2.07E-07 2.47E-09 2.13E-08 


KH20 0.088 — 8.75E-04* 4.81E-05 6.69E-05 1.14E-04 7.34E-08 ee 1.49E-08 
DdDm1 (061.94-41.3) Clegg87 0.100 1.40E-05 2.50E-05 1.40E-04 2.00E-05 2.90E-06 
Dinerstein03 1.40E-05 2.10E-06 


Milanova09 0.110 1.51E-05 9.77E-05 | 1.17E-04 


Otsuka09 0.102 1.05E-05 2.31E-05 1.18E-04 . 2.84E-05 2.23E-06 4.07E-08 . 6.10E-07 


KH20 0.093 1.12E-05 2.51E-05  1.15E-04 2.09E-05 2.04E-06 — 2.63E-08 | 6.46E-07 
H4-1 (049.3--88.1) Hawley78 0.106 ee 6.30E-05 2.20E-04 5.20E-06 4.00E-06 ee ee 
TP79 0.098 | 2.45E-03* 7.41E-05  3.16E-04 6.31E-06 7.94E-07 ee «1.58E-07 
Barker80 0.107 ee 6.31E-05 2.19E-04 5.25E-06 ee 4.90E-08 
Dinerstein03 «8.00E-06 <5.00E-07 


Milanova09 0.120 1.55E-03 3.55E-05 . 2.57E-04 


Otsukal13 0.108 1.04E-03 3.85E-05 1.50E-04 2.67E-06 1.86E-07 7.57E-09 . 3.63E-08 


KH20 0.122  1.68E-03* 5.16E-05  1.46E-04 2.95E-06 1.32E-07 Ue 2.71E-08 
K648 (009.8-07.5) Hawley78 0.100 1.80E-05 4.60E-05 2.50E-06 
TP79 0.098 ee 2.46E-06 | 6.61E-05  6.17E-06 <1.66E-06 sue «3.31E-07 
Barker80 0.100 ee 1.17E-06 | 4.47E-05 2.51E-06 1.82E-08 


Adams84 0.104 5.40E-04 3.00E-06 | 4.70E-05 . 5.00E-06 
Milanova09 0.100 1.35E-03 T.41E-05 — 8.32E-05 


Otsukal5 0.104 9.41E-04 2.28E-06 | 5.39E-05 = 2.75E-05 2.53E-07 . 3.76E-09 . 4.00E-08 


KH20 0.094  1.14E-03«  2.89E-06 | 4.89E-05 . 6.63E-06 1.85E-07 e 6.06E-08 
NGC 2242 (170.3+15.8) C Garnett89 0.100 — 9.77E-05 1.206E-05 5.89E-05 . 3.63E-05 1.51E-06 
TP90 0.100 2.45E-04*  5.25E-05 1.07E-04 6.31E-05 7.76E-07 
KH20 0.125 ari et 1.08E-05 1.36E-06 1.12E-06 

SaSt 2-3 (232.0+05.7) Pereira07 E ee 8.40E-05 . 3.60E-04 ee 3.10E-06 
Otsukal4 EE 5.24E-04 3.09E-05 1.70E-04 4.79E-05 1.48E-06 e 8.51E-07 
Pagomeno18 ee ee ee 2.60E-05 7.00E-07 Z 1.66E-06 

SBS 11504-599A (135.9+55.9)  Péquignot05 0.082 9.00E-05x 3.00E-05  4.50E-06 


Jacoby06 0.076 3.09E-05 8.71E-06 1.95E-05 4.17E-06 
Sandin10 0.0706 | 7.94E-05 2.95E-05 5.50E-06 9.12E-06 


Stasiiskal0 0.089 6.92H-05 1.41E-05  6.61E-06 6.76E-06 ^ «3.16e-7 ee <3.16¢e-8 
NGC 4361 (294.1+43.6) © TP90 0.107 9.12E-05 2.45E-05  6.46E-05 3.24E-05 Sad ae 1.23E-06 
Howard97 0.100 7.76E-05  2.24E-05  L.41E-04  1L.78E-05 2.00E-06 fee 2.51E-07 

Kholtygin98 ^ -..  2.90E-05 che 5.10E-05 x I 
PRTM 1 (243.8-37.1) © Pefia90 0.107 <3.98E-05 <1.00E-04 2.24E-04 7.94E-05 2.00E-06 
Howard97 0.100 3.98E-06*  3.16E-05 1.26E-04 2.29F-05 <1.00E-05 ne 3.39E-07 
PRMG 1 (006.0-41.9) © Pena89 0.091 Ads pue 1.26E-04  3.16E-05 6.31E-07 
Howard97 0.091  1.00E-04« 1.91E-04  1.15E-04  1.70E-05 1.00E-05 ns 3.16E-07 

JaFul (002.1+01.7) Jacoby97 0.141 s 9.33E-05  3.09E-04 <1.00E-04 4.17E-06 

JaFu2 (353.5-05.0) Jacoby97 0.115  «3.16E-05  «9.12E-05 5.37E-05 6.17E-06 <5.89E-06 sus «2.88E-07 
M2-29 (004.0-03.0) Howard97 0.093 3.55E-05* 1.82E-05 5.13E-05 8.71E-06 8.51E-07 "es 4.57E-07 
Pefia91 0.129 ss 9.55E-06 —2.04E-05 5.25E-06 8.13E-07 ss 1.82E-07 
Solar Asplund09 0.085 2.69E-04 6.76E-05 4.90E-04  8.51E-05 1.32E-05  3.16E-07 2.51E-06 


a References: Adams84=Adams et al. (1984); Asplund09=Asplund et al. (2009); Barker80=Barker (1980); Clegg87=C 
(1987); Dinerstein03=Dinerstein et al. (2003); Garnett9=Garnett & Dinerstein (1989); Hawley78=Hawley & M 


egg, Peimbert, & Torres-Peimbert 
iller (1989); Howard97=Howard, 


Henry, & McCartney (1997); Jacoby97=Jacoby et al. (1997); Jacoby06—Jacoby et al. (2006); Kholtygin98=Kholtygin (1998); Kniazev08=Kniazev 


et al. (2008); KH20=Kwitter & Henry (2020); Milanova09=Milanova & Kholtygin (2009); Otsuka09=Otsuka et a. 
(2010); Otsukal13=Otsuka & Tajitsu (2013); Otsukal4=Otsuka et al. (2014); Otsukal5=Otsuka, Hyung, & Tajitsu 


. (2009); Otsukal0=Otsuka et al. 
(2015); Pagomeno18=Pagomenos 


et al. (2018); Pefia89=Pena et al. (1989); Pertia90—Peiia et al. (1990); Pefia91=Pena, Torres-Peimbert, & Ruíz (1991); Péquignot05=Péquignot & 


Tsamis (2005); Pereira07—Pereira & Miranda (2007); Sandin10—Sandin et al. (2010); Stasiriskal0=Stasiziska et al. 


(2010); TP79=Torres-Peimbert & 


Peimbert (1979); TP81=Torres-Peimbert, Rayo, & Peimbert (2081); TP90=Torres-Peimbert, Peimbert, & Peña (1990) 


b An asterisk (*) indicates that the C abundance was derived from optical line strengths; unmarked abundances were derived from UV line strengths. 


^ Halo PN classification is uncertain, based upon one or more relatively low radial velocity observations (see 83.1.3). 
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Table 4. s-Process Abundances 


PN Author? Ge/H (32) Se/H (34) Br/H (35) Kr/H (36) Rb/H (37) Cd/H (48) Te/H (52) Xe/H (54) 
1C418 Sharpee07 ee ee ee 3.70E-09 Ue ee ee 4.66E-10 
Madonnal18 Z 2.51E-10 tee 5.89E-09 ee Zr 5.13E-10 
I1C2501 Sharpee07 tee tee Z 9.62E-10 aZ tee Z 9.40E-11 
1C4191 Sharpee07 e tee 1.25E-09 Z tee 1.81E-10 

IC5117 Sterling16 zx 4.14E-09 ee 7.13E-09 9.52E-10 1.71E-10 tee 
J900 Otsuka20 zx e tee 7.11E-09 7.28E-10 tee tee 2.61E-09 
NGC2440  Sharpee07 zz tee see 9.59E-10 tee tee tee 4.70E-11 
NGC3918  Garcia-Rojas15 ee 3.31E-09 ee 6.31E-09 3.24E-10 ee ee 3.09E-10 
Sterling17 Ue Ue 8.91E-10 Ue Ue Ue 
NGC5315 Madonnal7 ee 3.98E-09 3.39E-09 3.98E-09 <7.41E-10 ee ee 2.69E-09 
Sterling17 ee 1.90E-08 ee ee 
NGC7027  Péquignot94 tee ee ee 2.74E-08 e ee ee 3.36E-09 
Sharpee07 tee ee 7.34E-10 9.69E-10 tee tee tee 9.92E-11 
Sterling16 1.81E-09 4.72E-09 ee 1.16E-08 1.14E-09 1.08E-10 
Sterling17 ee ee 2.63E-09 ee Zr ee 
Madonnal8 ee 3.80E-09 4.37E-10 1.35E-08 1.07E-09 ee 5.62E-10 
TCL Aleman19 Z Zr tee 1.26E-08 
SMP47 Sterling17 Ue Ue Ue 6.61E-10 
SMP99 Sterling17 e E Zr 3.89E-09 


Multi-Object Survey 


120 pxeP Sterling15 tee 2.43E-09 tee 6.91E-09 
2.29E-08, 2.69E-10 ee 2.82E-08, 7.59E-10 


Solar Abundances 


SOLAR Asplund09 4.47E-09 2.19E-09 3.47E-10 1.78E-09 3.31E-10 5.13E-11 1.51E-10 1.74E-10 


a References: Aleman19—Aleman et al. (2019), Asplund09=Asplund et al. (2009), García-Rojas15—García-Rojas et al. (2015), Madonnal7=Madonna et 
al. (2017), Madonnal8=Madonna et al. (2018), Otsuka20—Otsuka & Hyung (2020), Péquignot94=Péquignot & Baluteau (1994), Sharpee07=Sharpee, 
Zhang, & Williams (2007), Sterling15—Sterling et al. (2015), Sterling16—Sterling et al. (2016), Sterling17=Sterling et al. (2017). 


b Top row: Median Se/H for 68 objects and median Kr/H abundances for 39 objects. Bottom row: Abundance maximum, minimum. 
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Table 5. Iron Abundances and Depletion Values? 


Shields78 Garstang78 Perinotto99 
PN Fe/H Depletion Fe/H Depletion Fe/H Depletion 
Cn 3-1 m e d - 6.0E-07 -1.72 
Hu 2-1 an b x ny 9.0E-07  -1.54 
IC 2165 1.00E-06 -1.51 
NGC 1565  6.31E-06 -0.70 
NGC 2022 1.58E-06 -1.30 ee E 
NGC 6543 vee a vee Ue 3.0E-06 -1.04 
NGC 6720 Ue Ue 7.00E-06 -0.70 
NGC 6741 2.51E-06 -1.11 
NGC 6886 2.51E-07 -2.11 
NGC 7009 Ue Ue 2.40E-05 -0.11 Ue mE 
NGC 7027  1.00E-06 -1.51 Ue Ue 4.0E-07 -1.86 
NGC 7662 Ue Ue 4.20E-07 -1.88 
a 1X gas/H] = log( X/ H)gas — log(X/H)o, where X represents the abundance of the 


element in question. 


Table 6A. Extragalactic Planetary Nebula Abundance Studies? 
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pE—— Number of Objects Elements Spectral Region o/g*C Comments 
FORNAX [dE4, -5] 
Danziger78 i He,N,O,Ar OPT 3.21E-04 
IC 10 [IBm, 10] 
Magrini09a 2 He,N,O,Ne,S,Ar OPT 1.30E-04 
LARGE MAGELLANIC CLOUD [SB(s)m, 9] 
Osmer76 3 He,N,O OPT 1.58E-04 
Dufour77 2 He,N,O,Ne,S,Ar OPT 2.88E-04 
Maran82 1 He,C,N,O,Ne,S,Ar UV, OPT 2.14E-04 
Aller83 8 He,C,N,O,Ne,S,Ar OPT 1.79E-04 
Leisy96 16 He,C,N,O,Ne,S,Ar UV, OPT 2.76E-04 
Dopita97 10 He,C,N,O,Ne,Mg,Si,S,Ar UV, OPT 2.46E-04 
Stanghellini05 24 C UV 
Leisy06 120 He,N,O,Ne,S,Ar OPT 2.04E-04 
Bernard-Salas08 18 Ne,S IR 
Chiappini09 110 He,N,O,Ne,S,Ar 1.70E-04 
Delgado-Ingladall 3 O,Fe OPT ,NIR Fe/O only 
M31 [SA(s)b, 3] 
Jacoby86 1 He,N,O,Ne OPT 4.67 E-04 disk 
Jacoby86 2 He,N,O,Ne OPT 2.35E-04 halo 
Jacoby99 12 He,C,N,O,Ne,S, Ar OPT 2.19E-04 bulge 
Jacoby99 3 He,C,N,O,Ne,S, Ar OPT 2.88E-04 disk 
Richer99 14 O OPT 3.80E-04 bulge 
Kwitterl2 16 He,N,O,Ne,S,CLAr OPT 4.17E-04 outer disk 
Sanders12 51 O OPT 2.82E-04 disk 
Sanders12 17 O OPT 3.16E-04 halo 
Balick13 2 He,C,N,O,Ne,S,CLAr OPT 3.29E-04 outer disk 
Corradil5 9 He,N,O,Ne,S,CLAr OPT zc 
Fang13,15,18 9 He,C,N,O,Ne,S,Ar OPT 3.06E-04 northern spur 
Fang13,15,18 7 He,C,N,O,Ne,S,Ar OPT 3.74E-04 giant stream 
M32 [cE2, -6] 
Jenner79 1 He OPT 
Richer99 4 O OPT 2.49E-04 
Richer08 11 He,N,O,Ne OPT 3.09E-04 


@ Galaxy morphological types were taken from NASA/IPAC Extragalactic Database https://ned.ipac.caltech.edu/, while the T-number for each 
galaxy was obtained from the RC3 catalog (de Vaucouleurs et al. 1991). 


b References: Aller83—Aller & Czyzak (1983); Balick13—Balick et al. (2013); Bernard-Salas08—Bernard-Salas et al. (2008); Chiappini09=Chiappini 
et al. (2009); Corradil5b=Corradi et al. (2015b); Danziger78=Danziger et al. (1978); Delgado-Ingladall—Delgado-Inglada et al. (2011); 
Dopita97=Dopita et al. (1997); Dufour77=Dufour & Killen (1977); Fang 13,15,18=Fang et al. (2013, 2015, 2018); Jacoby86=Jacoby & Ford (1986); 
Jacoby99=Jacoby & Ciardullo (1999); Jenner79=Jenner, Ford, & Jacoby (1979); Kwitter12=Kwitter et al. (2012); Leisy96=Leisy & Dennefeld 
(1996); Leisy06—Leisy & Dennefeld (2006); Magrini09a=Magrini & Gonçalves (20092); Maran82=Maran et al. (1982); Osmer76=Osmer (1976); 
Richer99—Richer, Stasiriska, & McCall (1999); RicherO8=Richer & McCall (2008); Sanders12=Sanders et al. (2012); Stanghellini05—Stanghellini, 
Shaw, & Gilmore (2005) 


€ O/H* values for studies of three or more objects are medians. For studies of one or two PNe, values correspond to the single value or the pair 
average, respectively. 
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Table 6B. Extragalactic Planetary Nebula Abundance Studies 


Auihetyenr? Number of Objects Elements Spectral Region O/ H*C Comments 
M33 [SA (s)cd, 6] 
Magrini09b 91 He,N,O,Ne,S,Ar OPT 1.64E-04 gradient 
Bresolin10 16 He,N,O,Ne,Ar OPT 2.75E-04 gradient 
M81 [SA(s)ab, 2] 
StanghellinilO 19 He,N,O,Ne,S,Ar OPT 1.86E-04 gradient 
NGC 147 [E5 pec, -5] 
Gongalves07 6 He,N,O,Ne,S,Ar OPT 1.29E-04 
NGC 185 [E3 pec, -5] 
Ford83 1 He,N,O,Ne,Ar OPT 7.59E-05 
Richer08 5 He,N,=O,Ne OPT 1.35E-04 
NGC 205 [E5 pec, -5] 
Richer08 13 He,N,O,Ne OPT 1.29E-04 
Gongalves14 14 He,N,O,Ne,S,Ar OPT 1.15E-04 
NGC 300 [SA(s)d, 7] 
Stasińska13 25 He,N,O,Ne,S,Ar OPT 1.95E-04 gradient 
NGC 3109 [SB(s)m edge-on, 9] 
Peña07 r4 He,N,O,Ne,S,Ar OPT 1.26E-04 
Flores-Durán17 ib He,N,O,Ne,S,Ar OPT 1.32E-04 
NGC 4449 [IBm, 10] 
Annibalil7 4 He,N,O,Ne,S,Ar OPT 2.00E-04 
NGC 5128 [SO pec, -2] 
Walsh99 2 N,O OPT 2.37E-04 
Walsh12 40 He,N,O,Ne,S,Ar OPT 2.88E-04 
NGC 6822 [IB(s)m, 10] 
Dufour80 1 He,N,O,Ne OPT 1.26E-04 
Hernández-Martínez09 11 He,N,O,Ne,S,Ar 8.91E-05 
Garcia-Rojas16 18 He,N,O,Ne,S,Ar OPT 7.11E-05 
SEXTANS A [IBm, 10] 
Magrini05 1 He,N,O,Ne,S,Ar OPT 1.3E-04 
SEXTANS B [IB(s)m, 10] 
Magrini05 5 He,N,O,Ne,S,Ar OPT 9.6E-05 
SAGITTARIUS [dSp(t),xxx] 
Zijlstra06 4d He,C,N,O,Ne,Mg,S,CLAr,K OPT 2.21E-04 
SMALL MAGELLANIC CLOUD [SB(s) pec, 9] 
Osmer76 3 He,N,O OPT 1.00E-04 
Dufour77 1 He,N,O OPT 3.80E-05 
Aller81 7 He,N,O,Ne OPT 1.38E-04 
Maran82 2 He,C,N,O,Ne,S,Ar UV, OPT 2.15E-04 
Leisy96 15 He,C,N,O,Ne,S,Ar UV, OPT 1.29E-04 
Leisy06 37 He,C,N,O,Ne,S,Ar OPT 1.00E-04 
Bernard-Salas08 T. Ne,S IR 
Stanghellini09 11 C UV no O/H 
Shaw10 14 He,C,N,O,Ne,S,Ar IR, OPT 1.10E-04 


a Galaxy morphological types were taken from NASA/IPAC Extragalactic Database https://ned.ipac.caltech.edu/, while the T-number for each galaxy 
was obtained from the RC3 catalog (de Vaucouleurs et al. 1991). 


b References: 


Aller81=Aller et al. (1981); 


Annibalil7=Annibali, Tosi, & Romano (2017); 


Bernard-Salas08=Bernard-Salas et al. (2008); 


Bresolin10=Bresolin et al. (2010); Dufour77=Dufour & Killen (1977); Dufour80=Dufour & Talent (1980); Flores-Dirán17—Flores-Dufan, Peña, 
& Ruiz (2017); Ford83—Ford (1983); García-Rojas16—García-Rojas et al. (2016); Gongalves07—Gongalves et al. (2007); Gongalves14=Goncalves 
et al. (2014); Hernández-Martínez09—Hernández-Martínez et al. (2009); Leisy96—Leisy & Dennefeld (1996); Leisy06—Leisy & Dennefeld (2006); 
Magrini05—Magrini et al. (2005); Magrini09b=Magrini, Stanghellini, & Villaver (2009b); Maran82=Maran et al. (1982); Osmer76=Osmer (1976); 
Richer08=Richer & McCall (2008); Shaw10=Shaw et al. (2010); Stanghellini09=Stanghellini et al. (2009); Stanghellinil0—Stanghellini et al. (2010b); 
Stasiriskal3=Stasiriska et al. (2013); Walsh99— Walsh et al. (1999); Walsh12=Walsh et al. (2012); Zijlstra06—Zijlstra et al. (2006) 


€ O/H* values for studies of three or more objects are medians. For studies of one or two PNe, values correspond to the single value or the pair 
average, respectively. 


d Zijlstra et al. (2006) include BoBn 1, usually assumed to be located in the MW halo (see Table 2), in their abundance survey of PNe in the Sagittarius 
dwarf spheroidal galaxy. 


Table 7. Median Gas Phase 
Abundances? for the Galactic 
Bulge and Disk 
Dust Type He/H N/O Ar/H 
Bulge 
Featureless (F) 
C-Rich (CC) 
O-Rich (OC) 11.04 -0.80 6.29 
Dual-Dust (DC) 11.09 -0.23 6.51 
Disk 
Featureless (F) 11.07 -0.65 6.21 
C-Rich (CC) 11.03 -0.70 6.07 
O-Rich (OC) 11.00 -0.81 6.04 
Dual-Dust (DC) 11.10 -0.27 6.56 
Solar? 10.93 -0.85 6.40 
“Number abundance ratios for He/H 


and Ar/H are given 


in the format 


12+log(X/H), while for N/O, the values 
are log(N/O). These values are taken from 
García-Hernández & Górny (2014). 


b Asplund et al. (2009) 
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Table 8. Disk Metallicity Gradients? From Planetary Nebulae 


Element Milky Way M31 b M33 b NGC 3000 M81Č 
O/H -0.276.035 -0.021.021 -0.342+.144 -0.159.058 -0.803.292 
Ne/H -0.242+.058 -0.041.021 -0.324+.270 -0.154.069 -0.686+.438 
Ar/H -0.207.173 -0.041.021 -0.279+.126 -0.270.074 

S/H -1.014.146 


a Gradients are given in units of dex/R25 


b Peña & Flores-Durán (2019) 


© Stanghellini et al. (2010b) 


Table 9. Least Squares Fit Parameters For Alpha Elements 


Relation y-Intercept Slope Correlation Coefficient Objects 
Planetary Nebulae 4 
Ne/H vs. O/H -0.056+.08 0.94+.083 +0.72 116 
S/H vs. O/H +0.024+1.17 0.77+.14 +0.47 117 
CI/H vs. O/H -2.26+1.03 0.83+.12 +0.56 107 
Ar/H vs. O/H -1.45+.83 0.91+.097 +0.66 117 
H II Regions and Blue Compact Galaxies) 

Ne/H vs. O/H -0.84+.27 1.02+.03 +0.96 85 
S/H vs. O/H -1.631.22 1.01.03 + 0.97 73 
C1/H vs. O/H -3.44+.11 1.00+.01 e 9 

Ar/H vs. O/H -3.03.75 1.10.09 +0.80 82 


Least-squares fits shown with solid lines to the PN data of KH20 displayed here. 


b Least-squares fits shown with dashed lines to a sample of H II region and blue 
compact galaxy data computed by Milingo et al. (2010) in the cases of Ne, S, 
and Ar. For Cl, we used the MW disk abundance gradients reported by Esteban, 
García-Rojas, & Pérez-Mesa (2015, 83.3) to derive the fit given here. 


Table 10. Properties of AGB Models 


SourceÜ ^ Metallicity Progenitor Mass (Mo) HBB Onset (Mo)0 


MONASH 0.007 1.0-7.5 3.00 
MONASH 0.014 1.0-8.0 3.25 
MONASH 0.030 1.0-8.0 3.50 
LPCODE 0.001 1.0-3.0 2.50 
LPCODE 0.010 1.0-3.0 2.50 
LPCODE 0.020 1.0-4.0 3.00 
ATON 0.004 0.0-8.0 2.50 
ATON 0.008 1.0-8.0 2.50 
ATON 0.014 1.0-8.0 3.00 
FRUITY 0.003 1.3-6.0 2.50 
FRUITY 0.008 1.3-6.0 2.50 
FRUITY 0.014 1.3-6.0 2.50 


a Each source name refers to one or more investigators and their associated 
publications, all of whom make use of a specific modelling code. The 
source names are explained in the text here as well as in Henry et al. 
(2018). 


b These values were estimated by plotting log(C/O) versus progenitor mass 
and determining the mass maximum where log(C/O) begins to descend 
as mass continues to increase. 


TT 


ESA/Hubble & NASA. 
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Figure 2. Aitoff projection showing the Galactic distribution of all 3540 PNe currently in HASH. True, 
likely, and possible PNe are the red, green, and black symbols, respectively. Credit: Fig. 1 in Frew (2016) 
(OCambridge University Press. Reproduced with permission of the Licensor through PLSclear. 
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Figure 3. Hertzsprung-Russell diagram of a complete evolutionary track for a 2 Meosolar metallicity star 
from the main sequence to the white dwarf phase. In the cooler section of the post-AGB phase, wiggles 
in the track are caused by numerical convergence difficulties. The blue track shows a born-again evolution 
model [triggered by a very late thermal pulse — see 84.2 in Herwig (2005)] of the same mass, however, shifted 
by approximately Alog Teff = —0.2 and Alog L/Lo = —0.5 for clarity. The red and green stars mark the 
position of the central stars of planetary nebulae for which spectra are shown in Herwig (2005, Fig. 8). The 
number labels for each evolutionary phase indicates the log of the approximate duration for a 2 Mọ case. 
Larger or smaller mass cases would have smaller or larger evolutionary timescales, respectively. T'his figure 
and its caption are republished with permission of Annual Reviews, from Annual Reviews of Astronomy and 
Astrophysics, Falk Herwig, 43, 435, 2005; permission conveyed through Copyright Clearance Center, Inc. 
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Figure 4. Crossing time vs. main sequence mass for solar metallicity. Crossing time is the interval between 
the star having an effective temperature ~10*K and when it reaches its maximum effective temperature. 
Based on data in Table 3 from Miller Bertolami (2016). 
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Figure 5. Median abundance values for disk PNe in Table 1 are plotted as 12 + log X/H for each element 
displayed along the horizontal axis for the surveys listed in the legend, using colored filled circles. Surveys 
are linked to their related references in a Table 1 footnote. Survey symbols are plotted left to right, early to 
late, in the order of survey date. Solar values from Asplund et al. (2009) are shown as black stars and are 


connected by a solid line. 
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Figure 6. Median abundance values for bulge PNe in Table 1 are plotted as 12+ log X/H for each element 
displayed along the horizontal axis for the surveys listed in the legend using colored filled circles. Surveys 
are linked to their related references in a Table 1 footnote. Survey symbols are plotted left to right, early to 
late, in the order of survey date. Solar values from Asplund et al. (2009) are shown as black stars and are 
connected by a solid line. 
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Figure 7. Median abundance values for halo PNe in Table 2 are plotted as 12 + log X/H. For each element 
ratio on the horizontal axis, each symbol in the plot represents the median abundance of all studies that 
pertain to an individual halo PN. Objects are identified by colored filled circles as shown in the legend. Solar 
values from Asplund et al. (2009) are shown with black filled circles and are connected by a solid black line 
to aid comparisons with PN values. 
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Figure 8. Values of 12+log(O/H*) versus morphological T-type for 19 external galaxies and the MW listed 
in Tables 6A,B. Each circular symbol represents the measured median value(s) of PNe within each galaxy 
as a function of its morphological T-Type. Symbols for individual galaxies are color-coded, and multiple 
symbols for that color represent measurements by different authors, as detailed in Tables 6A and 6B. 
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Figure 9. The relative percentages of dust types in the Galactic bulge, disk, and Magellanic Clouds. 
Galactic data are from García-Hernández & Górny (2014, Table 5), while the MC data are provided in 
Stanghellini et al. (2012, Table 6). Sample size for each region is indicated in parentheses on the horizontal 


axis. 
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Figure 10. Starting with the upper left panel and proceeding in clockwise fashion, Ne/H, S/H, Ar/H and 
CI/H, are plotted versus O/H, where values are expressed in the form of 12+log(X/H). Observational data 
(open circles) are from KH20 error bars represent median uncertainties and are shown in the lower right of 
each panel. We also show least-squares fits in each panel, where the solid line represents the fit to the PN 
data. For Ne, S, and Ar, the dashed line shows the fit to abundances for H II regions and blue compact 
galaxies compiled by Milingo et al. (2010, Table 6) for comparison purposes. The analogous fit for Cl was 
inferred from a study of Galactic disk H II regions by Esteban, García-Rojas, & Pérez-Mesa (2015, 83.3). 
'Table 9 displays the regression parameters and uncertainties. 
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Figure 11. A comparison of observations with model predictions of final AGB surface abundances. Upper 
panels: log(C/O) vs. He/H for Galactic disk PNe, with observational data taken from KH20. Observed C 
abundances forming the C/O ratios in the upper three panels were determined strictly from UV lines (C III] 
1909À and C IV 15494; filled circles) or optical lines (C II 4267A; open circles). Models produced by the 
MONASH (left), LPCODE (center), and ATOM (right) codes are shown, with model metallicity indicated 
by the color-coding defined in the legends. Each model within a metallicity set is represented by a cross (X) 
and has a unique progenitor mass, where the mass range for each set is given in Table 10. Lower panels: 
Same as the upper panels but for log(N/O) vs. He/H. Error bars for the observations represent median 
values of uncertainties and are shown in the lower right corner of the two left panels. The black dashed lines 
in the lower three panels define the upper-right quadrant containing Type I PNe. 
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Figure 12. Like Fig. 11 but for log(C/O) vs. log(N/O), showing AGB surface abundance predictions from 
the MONASH (ul), LPCODE (ur), ATON (ll) and FRUITY (lr) codes along with observational data of disk 
PNe from KB20. 
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Figure 13. Upper: Predicted values of final AGB surface abundances of He/H vs. progenitor mass in solar 
units from Karakas & Lugaro (2016); Karakas et al. (2018) for the three metallicity values indicated in the 
legend. Lower: Same as the upper panel but for log(N/O) versus progenitor mass. 
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Figure 14. Top: Comparison with proto-solar abundances of the contribution of LIMS plus rotating massive 
stars (red filled symbols), LIMS plus non-rotating massive stars (green crosses, only up to A—96 to avoid 
confusion at higher A), and LIMS only (red open symbols). Abundance uncertainties from Lodders, Palme, 
& Gail (2009) for each isotope are indicated. The dotted horizontal lines show deviations from the proto- 
solar values of 1096 and 5096. Middle: Contribution of LIMS to the total production of s-only isotopes. 
Bottom: Production factors normalized to °°Sm. (Note: This plot is a copy of Fig. 11 in the paper titled 
Chemical evolution with rotating massive star yields I. The solar neighbourhood and the s-process elements, 
Prantzos et al. 2018, Monthly Notices of the Royal Astronomical Society, 476, 3432. We are grateful to 
Oxford University Press for permitting us to display it.) 
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Figure 15. PNLF for PNe in M31 out to ^30 kpc. The characteristic sharp cutoff at the bright end occurs 
at m5097 ~20.17. Credit: R. Ciardullo. 
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Figure 16. The PNLF for five galaxies in the Leo I group (d~10 Mpc). Note the remarkably consistent 
bright-end cutoff across the different Hubble types. Credit: Based on a figure by R. Ciardullo. 
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Figure 17. Diffraction-limited V-band (A=554 nm) and Ng-band (A=645.9 nm) composite of Lo Puppis 
showing the dust ring and bipolar outflows in primarily scattered light. Note the stellar components (sepa- 


rated by ~33 mas) visible as white blobs above the narrow waist. The image is 38.4 AU on a side. Credit: 
ESO/P. Kervella. 
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Figure 18. Abell 39: Example of a round PN. The nebular diameter is roughly 1.2 pc. Image taken in 
[O IIIJ. Credit: WIYN/NOIRLab/NSF. 
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Figure 19. IC 4634: Example of a point-symmetric PN. The major axis from outermost arc to outermost 
arc is about 0.23 pc. Blue: H£, [O III]; Cyan: V; Red: Ha, [N H]. Credit: ESA/Hubble and NASA. 
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Figure 20. Mz 3: Example of a bipolar/bilobed PN. The major axis of the bulbous main nebula is very 
roughly 0.15 pc. Purple: [O III]; Blue: Ho; Green: [N II]; Red: [S II]. Credit: R. Sahai (JPL) et al., Hubble 
Heritage Team, ESA, NASA. 
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Figure 21. NGC 6543 main nebula. The major axis is roughly 0.25 pc. Blue: X-ray; Orange and Purple: 
optical. Note that X-rays are confined in the interior of the nebula. FLIERs are visible along the major axis. 
Concentric rings surround the main nebula. Credit: X-ray: NASA/CXC/SAO; Optical: NASA/STScI. 
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Figure 22. NGC 6543 extended halo. This halo has a diameter of about 1.7 pc. Blue: [O III]; Red: [N II]. 
Credit: Nordic Optical Telescope and Romano Corradi (Isaac Newton Group of Telescopes, Spain). 
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Figure 23. Toy model of a wind-blown bubble. From the star outward: stellar wind (gray); inner shock; 
shocked stellar wind (blue); ionized swept-up gas (the visible PN; red); ionization front; neutral swept-up 
gas (tan); outer shock; unaffected AGB wind (white). 
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Figure 24. NGC 7027. The main nebula is about 0.08 pc along the major axis. Blue: F343N ([Ne V]); 
Green: 502N ([O III]); and Red: F164N ([Fe II]). Note the concentric dust rings. Credit NASA, ESA, and 
J. Kastner (RIT). 
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Figure 25. NGC 6302. The blue-green “butterfly” nebula is about 0.9 pc in diameter. Blue: F343N 
([Ne V]); Green: F128N (Pa); Red: F164N ([Fe II]). The point-symmetic arcs emitting [Fe II] are suggestive 
of fast shocks. Credit: NASA, ESA, and J. Kastner (RIT). 
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Figure 26. NGC 7293 knots seen in Hg at 2.12 um. Image size: 5.10’ x 3.45’. The central star is outside 
of this image toward the bottom. Reprinted from Matsuura et al. (2009) Figure 2. 
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Figure 27. Two PNe with prominent FLIERs. The major axis of each nebula is about 0.2 pc. Left: 
NGC 6826. Blue: [O III]; Green: [O Ij; Red: [N II] and [S II]. Credit: Bruce Balick (University of Washing- 
ton), Jason Alexander (University of Washington), Arsen Hajian (U.S. Naval Observatory), Yervant Terzian 
(Cornell University), Mario Perinotto (University of Florence, Italy), Patrizio Patriarchi (Arcetri Observa- 
tory, Italy) and NASA/ESA. Right: NGC 7009. Violet: H8; Blue: [O II], [O III]; Green: He II, He I; Red: 
[O I], Ha, [N II]. Credit: ESA/J. Walsh. 
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Figure 28. Fg 1: The S-shaped jets appear in point symmetry across the nebula. The main nebula is about 
1.5 pc in diameter. Blue: [O II]; Green: [O III]; Red: Ha, [N II]. Credit: ESO/H. Boffin. 
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Figure 29. Sh2-188. Lines indicate the central star. This is one of the largest known PNe, with a diameter 
of about 2.8 pc. Original image credit: T.A. Rector/University of Alaska Anchorage, H. Schweiker/WIYN 
and NOIRLab/NSF/AURA. 
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Figure 30. Abell 21. Orange: Ho; Blue: [O III]. Lines indicate the central star. This is a large PN, 
with a diameter of roughly 1.6 pc. Original image credit: H. Schwelker/NOIRLab/NSF/AURA and T.A. 
Rector/University of Alaska Anchorage, and NOIRLab/NSF/AURA. 


